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1 Theoretical basics 



1.1 Introduction 

The radiative corrections (RC) for observables related to the Z resonance have been 
described in great detail in the 1989 CERN Yellow Report 'Z Physics at LEP V 
About 1.5 X lO'' Z decays have been recorded and analysed during the years of operation 
of the four LEP experiments — from autumn of 1989 to the end of 1994. In order to match 
the actual experimental precision, a completely revised analysis of radiative corrections at 
the Z resonance is needed. The aim of this contribution to the present Report is threefold. 
We will: 

(i) introduce a common language for presentation of the results emerging from different 
approaches; 

(ii) update the predictions of Z resonance observables within the Minimal Standard 
Model (MSM); 

(iii) estimate the intrinsic theoretical uncertainties of these predictions, which are mainly 
caused by the neglect of higher order contributions. 

The results, which are presented in this Report, are based on several different ap- 
proaches and on a comparison of their numerical predictions. 
The findings of the Report are based on the following computer codes: 

BHM 0, LEPTOP [§, TOPAZO §, WOH [§, and ZFITTER §. 
The design of some of these codes can be traced back to the 1989 Workshop on 'Z-Physics 
at LEP 1'; others have been developed later. All of them contain a so-called electroweak 
library, which allows the calculation of virtual higher order electroweak and QCD correc- 
tions to selected quantities for example, a weak mixing angle or a partial or total Z width. 
Strictly speaking, these are 'pseudo-observables' and not directly measurable in an exper- 
iment. Some packages (codes) contain, in addition, virtual and real photonic corrections 
(bremsstrahlung) with simple kinematical and geometrical cuts similar to those used in 
experiments, thus allowing for the calculation of (idealized) measurable cross-sections and 
asymmetries, which we will call 'realistic observables'. 

Apart from a comparison of the numerical results of different codes, we attempted a 
simulation of the theoretical uncertainties of each code. To do so we had to reach an 
agreement on the principles on which these simulations are based. We hope that as a 
result of this project the prospects of the quantitative tests of the MSM at LEP and SLC, 
and their sensitivity to potential New Physics beyond the MSM, will become clearer. 

The Report is organized as follows. In the first section we define the exact framework 
for the analysis of the electroweak data in the MSM, beginning with 'Input parameters' 
and ending with 'Basic notions', which are needed for discussing the numerical results 
collected in the subsequent sections. These deal with pseudo-observables and realistic 
observables in the sense introduced above. The first three sections represent a homogenous 
presentation of the material related to or obtained with the above mentioned codes. In the 
fourth and fifth sections we have collected items particular to the different codes, such as 
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explicit analytical formulae and descriptions of the design and of the use of the packages. 
All the authors of the report agree on the following conclusions from this study: 

• The differences between results of different codes are small compared to existing 
experimental uncertainties. Thus improvement of experimental accuracy at LEP 1 
and SLC is welcome even at the present level of theoretical accuracy. 

• At present the most promising are measurements of gylg^ various P- and C- 
violating asymmetries and polarizations. 

• The real bottleneck for improved theoretical accuracy in 5^^/(7^ is presented by 
the uncertainty of the input parameter a = a{M^). The improved accuracy of 
this important parameter calls for new accurate measurements of the cross-section 
e+e~ — > hadrons at low energies (Novosibirsk and Bejing accelerators, etc.) 

• The estimates of theoretical uncertainties are highly subjective and their values 
partly reflect the internal philosophy of the actual implementation of radiative cor- 
rections in a given code. 

• In many cases the one-loop approximation in the electroweak gauge coupling is ad- 
equate enough at the present level of experimental accuracy. At the same time, 
however, it should be stressed that a complete evahiation of the sub-leading correc- 
tions, 0{G'j^M'^m^) would greatly reduce the uncertainty that we observe, one way 
or the other, for all observables. 

• In case the next generation of experiments at LEP 1 and SLC improves accuracy 
considerably (a problem not only of statistics but mainly of systematics) the full 
program of two- loop electroweak calculations should be carried out. 

1.2 Input parameters 

Within the MSM, any measured quantity can be calculated in terms of a small set of 
input parameters. Once all possible relations between the parameters of the Lagrangian 
in the MSM are exploited, we may choose a set of them which, in the electroweak part, 
consists of one interaction constant — say. the fine structure constant a — of the masses 
of all particles and of the CKM fermionic mixing angles. The latter ones are of little 
importance for Z resonance observables, which are dominated by flavour-diagonal neutral 
current interactions. 

In practical calculations one prefers to make use of the most precisely measured pa- 
rameters. Three of them characterize the gauge sector of the MSM: 

a = a{0) = 1/137.0359895(61), 
Gf, = 1.16639(2) X 10-^ GeV-^ 

= 91. 1887 ± 0.0044 GeV, (1) 

where G^, the four-fermion /x decay coupling constant, is defined through the muon life- 
time by 
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Light fermion masses contribute via the electromagnetic couphng constant a, running 
from very small momenta up to M^, the typical scale of the Z resonance, and yielding 



a 



1 - Aa ' 

with Aa consisting of leptonic and hadronic contributions 

Aa = Aai + Aat ■ 
The leptonic contribution is known explicitly: 



a r 5 mf ^ ( ^ r.^f\,l3i + l 



(3) 



(4) 



(5) 



where 



ft 



1-4 



M2 ■ 



(6) 



With lepton masses taken from the last edition of the Review of Particle properties 0, 
from (151) we get 



Aai = 0.0314129 . 

In the small mass approximation, (^) reduces to the well-known expression 



(7) 



Aai = — yfln 



m? 3 



0.0314177 . 



The hadronic contribution is being calculated by a dispersion relation, 

ds 



Aah = -yile 
37r 



Aral ■^{M'i — ie — S 



-R{s) 



(9) 



from experimental data on the cross-section ratio R{s) = o"e+e- (s) / ctq^s), where cXe+e- (s) 



to 



r 



hadrons) and ctq = 4/37r« /s. In 1990, the value Aah has been updated 



Aah = 0.0282(9) , (10) 

leading to a-\M^) = 128.87 ± 0.12 § Q. 

By using a instead of a, the electromagnetic coupling constant at the correct scale 
for Z physics is chosen and one automatically avoids the problem of the light quark mass 

Very recently, three new analyses have been published. The first, @, leads to Aan = 0.02666 ± 
0.00075 and Qf-^M^) = 129.08±0.10. The second, leads to Aa^ = 0.02732±0.00042and a~i(M^) = 
128.99 ±0.06. The third, leads to Aa,, = 0.0280 ± 0.0007 and a~^(iV4) = 128.899 ± 0.090. 
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singularities. So the relevant input parameter is a, and not a in — spite of the extremely 
high accuracy of the latter. It should be stressed that the quantity a used by us is not 
only numerically, but also in principle, different from the quantity a{Mz) = 1/127.9(1) 
(see, for example Ref. 0, page 1304). The latter is defined in the modified minimal 
subtraction scheme (MS) in terms of a bare electric charge, while the former is a physical 
quantity expressed in terms of charged fermion masses. 

The treatment of the light fermion masses (all but the top quark mass) in the elec- 
troweak corrections is not the same in different codes. First we should mention that the 
exact expression (|^) for the photonic vacuum polarization is not the only place where 
fermionic power corrections, m^/M^, may appear]^. There are fermionic contributions 
to the self-energy insertions of the heavy gauge bosons, which also contain power mass 
corrections. Moreover additional power corrections appear, due to the presence of axial 
couplings. The latter may be exphcitly calculated for leptons; for quarks, they are not 
taken into account by the dispersion integral (P). In principle, there are also finite mass 
corrections in three- and four-point functions. In some codes these finite mass terms are 
neglected completely, some other codes retain them in the two- and three-point functions. 
In the latter case, for the quarks some effective masses are used, which are selected in 
order to reproduce with a sufficiently high accuracy the hadronic vacuum polarization 
contribution Aoh 0. We emphasize that these finite mass terms have no numerical rele- 
vance. They are being retained just in order to investigate the sensitivity of corrections to 
finite light fermion masses. They also help sometimes to compare results of independent 
calculations. Such a comparison has been done once for Ar, and an agreement of up to 



12 digits (computer precision) was found pA 



The T lepton mass has been retained in all codes in the phase-space corrections and 
in the Born-level matrix element. The masses of the heavy quarks b and c also give 
rise to phase-space corrections, but are treated together with QCD final state interaction 
corrections (see the QCD part of this Report). We take the following input values for rrir, 
nib and rric [0: 



rur = 1.7771 ±0.0005 GeV, 
rrib = 4.7 ±0.3 GeV, 

rric = 1.55 ±0.35 GeV, (11) 

where m^, rric are the pole masses converted in the actual calculation of radiative correc- 
tions to MS masses using QCD perturbation theory. The meaning of these MS mass 
definitions and their relation to the pole masses are explained in the QCD Part of this 
Report by Chetyrkin, Kiihn and Kwiatkowski. The QCD part of the MSM is character- 
ized by these running quark masses, and by the strong coupling constant ^^(M^), which 
will be assumed to be [|1^ and [p!6[| : 



^^(M^) = = 0.125 ±0.007 . (12) 



^One should emphasize that the dispersion integral (Q) automatically includes all these power correc- 
tions for the hadronic part of the photonic vacuum polarization. 
^Details on effective quark masses be found in |12[ or p3[. 
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Note, that from pure QCD observables at the Z peak one obtains 



ag(M^) = = 0.123 ±0.006 . (13) 

Finally, there are two yet unknown input parameters left: the top quark mass and the 
Higgs boson mass, 

mt, M„, (14) 

although the situation has considerably improved with the recent CDF indication [1^ of 
evidence for the t quark with rrit = 174 ± 17 GeV. In the following, by rrit we will always 
imply the t-quark pole mass. 

1.3 Codes to calculate electroweak observables 

The radiative corrections for the measured physical observables must be included into 
the theoretical predictions for them, and the consistency of the theory is verified by 
comparison with the data. This particular step may assume different aspects related to 
the actual implementation of the comparisons, but in the end it usually takes the form of 
some constraint on rrit and and to a lesser extent on and a. It is our main goal to 
present and discuss the most accurate theoretical predictions for the measured quantities 
and to introduce a reliable estimate of the associated theoretical uncertainties. 

The results of this Report have been obtained by a critical comparison and examination 

of: 

- the variations in the results following from the different, although not antithetic, 
formulations of the various groups; 

- the internal estimates of the uncertainties induced by the still missing higher order 
corrections. 

The groups participating in this project can be identified by the names of the codes 
that they have assembled, namely: 

BHM 

Burgers, Hollik, Martinez, Teubert; 

LEPTOP g] - ITEP Moscow group 
Novikov, Okun, Rozanov, Vysotsky; 

TOPAZO 0] - Torino-Pavia group 

Montagna, Nicrosini, Passarino, Piccinini, Pittau; 
WOH H 

Beenakker, Burgers, Hollik; 
ZFITTER - Dubna-Zeuthen group 

Bardin, Bilenky, Chizhov, Olchevsky, S.Riemann, T.Riemann, Sachwitz, 
Sazonov, Sedykh, Sheer. 
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All these codes may be used to calculate pseudo-observables. Three of them, BHM, TOPAZO, 
and ZFITTER, may calculate also realistic observables. One code, TOPAZO, is additionally 
able to calculate the full Bhabha cross-section, including complete s and t channel ex- 
change contributions. Concerning the electroweak corrections, four codes are based on 
completely independent theoretical computational schemes, while BHM/WOH use basically 
the same framework. The treatment of QED corrections in the three QED dressers is 
based on completely independent theoretical methods. As a consequence, it is extremely 
difficult to describe them in a common language. On the other hand, the treatment of 
QCD corrections is to a large extent common in all five codes. It is based on papers 
included in Part II of this volume and on references quoted therein. 

In this Report, we will attempt to present the collected material homogeneously as for 
long as possible. Because of this, all complicated description of what is entering in the 
various, alternative theoretical formulations, has been shifted to the last two sections of 
the Report. 

1.4 The language of effective couplings 

A natural and familiar but approximate language for the basic ingredients of the physical 
observables of the Z resonance is that of effective couplings. How it will translate into the 
particular realizations and schemes of the various codes can be studied in References [0- 
1^ and in the references quoted in these. The more theoretically oriented reader should, 
however, be aware that basic differences do indeed exist and that the language of effective 
couplings is not universally realized in all the approaches. The basic notions of effective 
couplings, which are common to all approaches, can be easily introduced. Nevertheless, 
this deserves to be seen in several logical steps. 

1.4.1 The Born Approximation for e"'"e~ ff 

The matrix element of the process 




Figure 1: Feynman graph for the reaction (^) in the Born approximation 
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s 
1 
s 



QeQfla ® 7" + X (^^:^.^7a ® 7° - ^:^f7a ® 7^75 
algllalb ® 7" + 9ld7a75 ® 7^75 



(16) 



(17) 



with X being the propagator ratio 

X = 



s- Ml + isVjM^ ■ 
And in the Minimal Standard Model (MSM) 



r(3) 
'/ 



r(3) 



2Qf sin^ e 



(19) 



where (^f , gl are vector and axial vector couplings of the Z-boson, are the electric 
charges of fermions in units of position charge, Jj- the projection of weak isospin, and 
9 the weak mixing angle in the Born approximation. Its value depends on the choice 
of computational scheme. In ([161) and ([T7D a short notation for bilinear combinations of 
spinors u and v is used: 



(20) 



In the matrix element (IB) and (|T^ the contributions from 7 and Z exchange diagrams 
are unambiguously separated, and the Z exchange contribution is presented in a factorized 
form (ITf 



1.4.2 Electroweak non-photonic corrections 

The higher order electroweak non-photonic corrections are indicated symbolically by the 
blobs in Fig. They include all possible self-energy and Zff vertex insertions, together 
with non-photonic boxes insertions (the WW and ZZ boxes in the one-loop approxima- 
tion), and lead to a slightly more complicated structure in the matrix element, which is 
valid as long as we neglect the external fermion masses: 




+ 




+ 




Figure 2: Feynman graphs for higher-order corrections to the reaction (ITH) 
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s 



T:1{s, t)7. ® 7" - (5, t)7a ® 7^75 

-^llis, thai, ® 7° + -^at (S, thai, ® 7"75l | • 



(21) 



The products of interaction constants in {\i7\) are replaced by the running a{s) — this 
time including the imaginary part — and by four running electroweak form factors, T^ - . 
In the MSM, the explicit expressions for the form factors result from an order-by-order 
calculation and are certain explicit functions of the input parameters (Q), ([TTp, ( [T^ ) 
and dll). 

Several comments should be given on the structure of the matrix element (|2T|) after 
inclusion of higher order electroweak corrections: 

• This structure is unique, but implies the introduction of complex valued form factors 
which depend on the two Mandelstam variables s, t; the dependence on t is due to 
the weak box diagrams. 

• The separation into insertions for the 7 and Z exchanges is lost. In the proposed 
realization, the first term in (|2l|) contains the running QED coupling where only 
fermionic insertions are retained. This ensures a gauge invariant separation. 

• The weak boxes are present in (pT]) as non-resonating (~ s — M^) insertions to the 
electroweak form factors J^tj ■ At the Z resonance, the one- loop weak WW and ZZ 
box terms are small, with relative contribution < 10~^. If we neglect them, the t 
dependence is turned off. The t-dependence would also spoil factorization of the 
form factors into products of effective vector and axial vector couplings. 

• Full factorization is re-established by neglecting, in addition, the other non-resonating 
loop contributions, such as the bosonic insertions to the photon propagator, and 
photon-fermion vertex corrections. All the neglected terms are of the order 0{ajf-). 
The resulting effective vector and axial -vector couplings are complex valued and 
dependent on s. The factorization is the result of a variety of approximations, valid 
at the Z resonance to the accuracy needed, and indispensible in order to relate the 
pseudo-observables to actually measured quantities. 

In the complete codes — BHM , TOPAZO , ZFITTER — it is possible to control the nu- 
merical influence of all of these approximations. 



1.4.3 Z pole approximation 

After the above-mentioned series of approximations, we arrive at the so called Z boson 
pole approximation, which is actually equivalent to the setting s = in the form factors. 
In the Z pole approximation, the one-loop diagrams contributing to the blobs in Fig. ^ 
can be visualized in the diagrams of Fig. 
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(a) 



(b) 



+ 




+ 




+ 




+ 




Figure 3: The one-loop Feynman graphs at the Z boson pole 



In Fig. ^ (a) all particles (Z, if, W ^ g, /, v) contribute through their weak neutral currents, 
while in the Z7 mixing (b) only charged ones contribute. This, of course, is strictly true 
only in the unitary gauge, while in the renormalizable gauge, Faddeev-Popov and Higgs- 
Kibble ghosts will also appear. The vertex diagrams (c) and (d) contain Z, W ^ as vertical 
lines. The vertex diagrams (e) and (f) contain trilinear gauge boson interactions. The 
vertex diagrams contain, in principle, also Higgs boson exchanges — as vertical lines in 
diagrams (c) and (d), and as Z^H virtual states in diagrams (e) and (d). They give, 
however, a negligible contribution because of small Hjf Yukava couplings 0. 

The set of the diagrams of Fig. ^ has to be complemented by the contributions from 
the bare 'unphysical' parameters (masses and couplings) in order to get the physical 
amplitude of Fig. |^. This is depicted in Fig. ^ as an example, the cross representing the 
contribution from the bare mass M°. 

The W boson self-energy enters the analysis only through the W mass (see Fig. 
The W self-energy contains all types of particles entering through their charged weak 
currents. Explicit expressions for the decay width of the W are not discussed in this 
Report, (see the existing literature [T^). 

If the mass of the W boson were known with the same accuracy as that of the Z, 
the conceptual picture of the electroweak corrections would be as simple as that of QED. 

^In some renornialization schemes one should add also the self-energy insertions for external fermion 
lines. In some other schemes, they vanish, when the field renormalization counterterms contributions are 
added. 
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z z 



+ 




a b 
Figure 4: The Z-boson self-energy 



W W 




a b 
Figure 5: The W boson self-energy 



Every physical observable would be expressed in terms of a, M^, M^, the masses of all 
fermions and the mass of the Higgs boson. Unfortunately, is not known with an 
accuracy comparable to that of M^. Therefore, instead of we have to choose as 
the most accurately measured dimensional observable. 



1.4.4 Effective couplings for Z decay. 

Formulae (0) and (|21| ) are very similar in their structure, a property which makes the 
language of effective couplings so convenient. These formulae refer to a general // — >■ /'/' 
annihilation process. After the simplifications described in 1.4.2-1.4.3, the factorization 
allows us to write for the form factors 

Tt^is = Ml,t = 0) = gt{Ml)gf {Ml) . (22) 

This implies that the Z part of the amplitude (0) is obtained from that of (0) by 
substituting 

gli - QliiMl) . (23) 

Simultaneously, we obtain for the matrix element of the decay process Z ff, Fig. ^, 
as follows: 



^ Zff 



M/7o 



(24) 



where Uf, Vf, are wave-functions of the fermion /, antifermion / and Z boson. 
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Figure 6: Feynman graph for higher-order corections to the Z ^ ff decay 

In the one-loop approximation, the blob in Fig. ^ can be visualized as the diagrams of 
Fig. 1^. Note that Fig. |^ (a) contributes both to initial and final state virtual corrections. 
The effective decay couplings Ql^, are nothing else but the reduced electroweak form 
factors, which are now constant because of the kinematical constraint s = M^, which is 
equivalent to the Z pole approximation. They will be the basic objects in the following 
presentation. We will use for them the simplified notations, 

= T^eQUMl) , (25) 

where the capital letters V, A are introduced in order to distinguish effective couplings 
Qy^, dressed by higher order interactions, from their Born-like analogs g^^a- 




Figure 7: Feynman graphs of the Z decay in the one-loop approximation 

In the next two subsections we will define the quantities to be computed and compared: 
the pseudo-observables, and the realistic observables. 
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1.5 Pseudo-observables 



The pseudo-variables are related to measured cross-sections and asymmetries by some de- 
convolution or unfolding procedure. The concept itself of pseudo-observability is rather 
difficult to define. One way to introduce it is to say that the experiments measure some 
primordial (basically cross-sections and thereby asymmetries also) quantities which are 
then reduced to secondary quantities under some set of specific assumptions. Within these 
assumptions, the secondary quantities, the pseudo-observables, also deserve the label of 
observability. Just to give an example, we quote the de-convoluted forward-backward 
asymmetry, where, typically, only the Z exchange is included and initial and final state 
QED corrections, plus eventual final state QCD corrections, are assumed to be subtracted 
from the experimental data. We have analyzed 25 of such pseudo-observables in details, 
namely: 



mass of the W 

hadronic peak cross-section ah 

partial leptonic and hadronic widths F,^, Fg, F^, F^, F^, F,^, Fc, F^ 

the total width F^ 

the total hadronic width F/j 

the total invisible width Finv 

ratios Ri,Ri),Rc 

asymmetries and polarization ^4^^^ , A^^ , A'^^^ , A'^^ ,P'^,P^ 

effective sine sin^ ^^g**, sin^ 9^^ 



The effective sine are defined by 



4\Qf\sm'eis^l-J, (26) 

with Qf being the electric charge of the fcrmion / in units of the positron charge. By 
definition, the total and partial widths of the Z boson include final state QED and QCD 
radiation. The explicit formulae for partial widths will be presented in subsection 1.9. 
Moreover, we have defined 



Fh = r„ + rd + r, + r, + r5, (27) 

Tinv — — Fe — F^ — Ft- — F/i , (28) 

Ri - ^, (29) 

i e 

Rb,c = ^, (30) 

i h 



^With lepton universality there exists only one leptonic effective sine, sin^ ^^.g"* = sin^ ^^g. 
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= 12-^- (31) 

z z 

The quantity at is the de-convoluted hadronic peak cross-section, which by definition 
includes only the Z exchange. To this end we would like to emphasize that in our calcu- 
lations we indeed assumed that 

Tinv — , (32) 

then the total Z width becomes 

= 3r, + Te + + + . (33) 

Unlike the widths, asymmetries and polarizations do not contain, by definition, QED and 
QCD corrections; furthermore, they will only refer to pure Z exchange: they are nothing 
but simple combination s of the effective Z couplings, introduced above, 

= ^^ 

P..{r) = (34) 

where we define 

and (7^ , g^^ defined by Eq. (|2^) are the effective neutral current vector and axial- vector 
couplings of the Z to a fermion pair //. Ppg(r) is the r polarization forward-backward 
asymmetry. One should realize that as a consequence of the adopted definition of pseudo- 
observables, there exist even more relations among them than are indicated above; it is, 
for example, = — P^, if in addition lepton universality is assumed, which is granted 
automatically in the MSM. For this reason we will present numerical and graphical results 
only for a selected subset of them. 

1.6 Realistic observables 

Realistic observables are the cross-sections cr-^(s) and asymmetries [s) of the reactions 



e+e-^(7,Z)^//(n7), (36) 

for a given interval of s = around the Z resonance, including real and virtual photonic 
corrections ('dressed by QED'). We will present results for / = /x and / = had, both in a 
fully extrapolated set-up or with simple kinematical cuts. For / = e both the s channel 
exchange cross-section and the complete Bhabha reaction will be treated. In total, we 
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considered 14 cross-sections and asymmetries for different combinations of flavour / and 
cuts; they will be defined in section 3. The kinematical cuts that will be imposed are of 
two types: 



and: 



^min ^ ^— ^ ^ ^min 



^min ^ ^— ^ ^ ^min 



• ^acoU < ^ 

• E± > Eth ■ 

Where s' is the invariant mass of two final-state fermions, 6- refers to the outgoing scat- 
tering angle of the outgoing fermion, 6'acoii is the acoUinearity angle between the outgoing 
fermions and E± the outgoing fermion (anti-fermion) energies. By s channel Bhabha 
scattering we mean that the s — t and t — t exchange interferences have been subtracted 
from the observable under consideration. 

We would like to mention in passing that the QED dressers have options which allow a 
model independent interpretation of realistic observables. Since this is beyond our present 
scope, we refer the reader for details to the existing literature [p0[]-[pl|. 



1.7 Calculational schemes 

Before entering into a detailed study of the numerical results it is important to underline 
how an estimate of the theoretical uncertainty emerges from the many sets of numbers 
obtained with the five codes. First of all, one may distinguish between intrinsic and 
parametric uncertainties. The latter are normally associated with a variation of the input 
parameters according to the precision with which they are known. Typically, we have 
\Aa-\Ml)\ = 0.12, |Amfe| = 0.3 GeV, {Amd = 0.35 GeV etc. These uncertainties will 
eventually shrink when more accurate measurements become available. In this Report we 
are mainly devoted to a discussion of the intrinsic uncertainties associated with missing 
non-leading higher-order corrections, although some results on parametric uncertainties 
will be also given. An essential ingredient of all calculations for radiative corrections to 
pseudo-observable is the choice of the renormalization scheme. There are many renormal- 
ization schemes in the literature: 



the on-shell schemes in various realizations |^-|^, ||^ 



the scheme [^-^\ 
the * scheme 



the MS scheme p3|-|3l 
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One cannot simulate the shift of a given quantity due to a change in the renormahzation 
scheme with one code alone. Thus the corresponding theoretical band in that quantity 
will be obtained from the differences in the prediction of the codes, which use different 
renormahzation schemes. On top of that we should also take into account the possibility 
of having different implementations of the full radiative corrections within one code, — 
within one well specified renormahzation scheme. Typically we are dealing here with the 
practical implementation of resummation procedures, of the exact definition of leading 
versus non-leading higher order corrections and so on. Many of these implementations are 
equally plausible and differ by non-leading higher order contributions, which, however, 
may become relevant in view of the achieved or projected experimental precision. This 
sort of intrinsic theoretical uncertainty can very well be estimated by staying within each 
single code. However, since there are no reasons to expect that these will be the same in 
different codes, only the full collection of different sources will, in the end, give a rehable 
information on how accurate an observable may be considered from a theoretical point of 
view. 

Of course, only a complete two-loop electrowcak calculation, combined with some 
re-summations of the leading corrections (renormahzation group improvements), would 
ultimately solve the problem of the theoretical accuracy. 

1.8 Main features of different approaches 

In this subsection we will discuss common features of and main differences between the 
electroweak libraries of the five codes. 

Common features: 

• All five codes use as input parameters the most accurately known electroweak pa- 
rameters G^,M^ and a, (the codes, with the exception of LEPTOP, also use a), in 
order to calculate the less precisely measured pseudo-observables. 

• They use the same expressions for final state QED and QCD corrections (radiation 
factors) . 

• All codes include essentially the same internal gluon corrections of the order of aag 
in the W and Z self-energy quark loops. 

• All codes include leading two-loop corrections of the order of Cj^m^; all gluonic 
corrections of the order of asGf^m^; leading gluonic corrections aa^ in the vector 
boson self-energies. 

Main differences: 

• Each code uses a different renormazation framework. 

• Some codes define an electroweak Born approximation, others give no physical em- 
phasis to a Born approximation and employ only the notion of an Improved Born 
Approximation (IB A) which includes the leading electroweak loop corrections. 
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• Certain codes include higher-order electroweak corrections, which are documented 
in the hterature but are numerically irrelevant, such as the irreducible two-loop 
photon vacuum polarization and the Higgs corrections proportional to a^m^ . Some 
codes use, on top of these higher order corrections, resummation of one-loop terms. 

• They differ by the choice of the definition of the weak mixing angle. 



1.8.1 BHM/WOH, ZFITTER 

The three codes BHM 0, WOH ^ and ZFITTER rely on different realizations of the 
on-mass-shell renormalization scheme. They perform the renormalization procedure with 
systematic use of the counterterm method for the basic parameters a, M^, M^. The BHM 
and WOH codes lead back to one approach and are thus not completely independent of 
each other. 



The description of the renormalization schemes may be found in p8| for BHM/WOH and 

in 



24 for ZFITTER. The weak mixing angle formally appears as 



cl = cos'e^ = {M^/M,y, 

si = sm'e^ = l-cl. (37) 



This corresponds to the definition proposed by Sirlin The masses in ( pTj) are the 

physical W and Z boson masses. Thus, is not an independent quantity and is used 
mainly for internal bookkeeping. For fixed values of rrit and each observable, including 
the Fermi constant G^, is expressed in terms of a, and with the corresponding 
quantum corrections being taken into account. Because of lack of precision in the W 
mass, M^, it is not taken as an experimentally measured input quantity. It is instead 
replaced by the more accurate value of G^. 

The quantum corrections for fi decay after removing the Fermi-model like QED 
corrections are contained in the parameter related to by the equation^ 

Since this equation contains a instead of a{M^) a large fraction of pc is of purely electro- 
magnetic origin via Aa. The parameter pc in (^) takes into account the W propagator, 
vertex and box corrections due to one-loop diagrams and presently available higher order 
terms. The one-loop corrections were first calculated in Refs. |^5|-|^. The correction 
due to the photon self-energy is QED renormalization group improved: a geometrical 
progression brings Aa into the denominator of pc- Hence, pc is usually written in the 
form 

Pc = , (39) 

1 — Ar 



^The naive Born approximation for muon decay of Eq. (38), pc = 1, was used in the 70's for predictions 
of the masses of the W and Z bosons with a crude accuracy by using the value of s,^ from the ratio of 
neutral and charged currents in neutrino interactions. 
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where Ar contains all the one-loop corrections to the muon-decay with the inclusion and 
the proper arrangement of the higher order terms. (For more details see the explanations 
presented in subsection 1.10.2.) 

For given values of rrit and , Eq. (|38D fixes 9^ and, hence, by the experimental 
value of Gf^. In practice, Eq. (|38D is solved with respect to iteratively, since Ar is a 
complicated function of M . This equation for iteration of reads 



Ana 



\ ^ V2G,Ml[l-Ar{M^)] 



(40) 



Therefore, appears as an rrit, M^, dependent prediction, which can be compared 
with the experimental values from UA2 and CDF [p^ . 

After fixing M^^ in this way all the other observables are expressed in terms of 
Gfj,, a, M^,mt, and d^Q. As a consequence of this procedure, purely leptonic Feynman 
graphs (e.g. a vertex correction in the Z ^ fifi decay) also turn out to be implicitly rrit 
dependent through the M^. 

In conclusion, Eq. (^) uses the amplitude of /i decay with the inclusion of weak 
corrections in order to establish the interdependence between M^, rrit and and the 
best measured parameters a,G^,M^. Thus, appears as a prediction as well as an 
intermediate parameter for the calculation of Z-boson observables. 

The flowchart of the BHM/WOH, ZFITTER approach is shown in Fig. 8. One should 
have in mind that in spite of a common presentation of the basics of these codes there 
are certain differences in the realization of the on-mass-shell renormalization schemes 
between BHM/WOH on one side and ZFITTER on the other — for example, in the use of 
different gauges and in the different treatment of field renormalization. 



'' In leptonic processes the dependence on cts is of higher order; it comes from gluonic exchanges 
between virtual quarks in W, Z, Zj and 7 self-energy diagrams with quark loops. 
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FLOWCHART OF ZFITTER/BHM/WOH 

Select minimal set of parameters in the MSM Lagranian: 
ao, ^wo' ^zo' ^ffo' "^/o (including rrito)] note that a^^, 
and VEV 77 are not among these. 



Define renormalization Z-factors for each bare parameter and each field 
(Z-matrices for Z— 7 and fermion mixing — for ZFITTER only). 
Fix Z-factors on mass shell. Use dimensional regularization (1/e, /i). 



Lagrangian now depends only on physical fields, couplings and 
masses, and on counterterms (Z-factors). 



Expand Z-factors; Zj= l+a/j, where a = a{0) and //s are functions 
of physical input M^, M^, M^^, nif and 1/e and fi. 



Calculate one-loop electroweak amplitudes with graphs, including 

loops and coimterterms; 1/f and // drop out. 



Improve one-loop results by RG-techniques and by proper resummation 
of the higher-order e.w. terms. Define improved Born approximation. 



Select experimental inputs: a{0), M^, (r^). 



Cet from G^= {it/ a/2) (a/s^ cj^ M^) pc, where pc depends 
on mu M^, a(0), M^, M, and 4 = 1- /M^. 



Calculate Z decay observables, with rrit and M^j free, 
in terms of G^, a(0), M^. 



Introduce gluonic corrections into quark loops and QED + QCD 
final state interactions in terms of a, as{M^), rrib (M^), rrit. 



Compare the results with electroweak experimental data, 
exhibit M^, rrit, M^j, and 0.3 (M^) dependence. 

Figure 8: BHM/WOH ZFITTER flowchart 
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1.8.2 LEPTOP 



The authors of LEPTOP avoid the use of counterterms in the formulation of the theoret- 
ical framework. They do not use resummation of one-loop electroweak corrections, thus 
limiting themselves to a simple one-loop approximation, completed by selected leading 
electroweak two-loop corrections. In contrast to all the other codes, LEPTOP uses a Born 
approximation. According to LEPTOP, the weak mixing angle 9 is defined by (s = sin 9, 
c = cos 9) : 



Ga = , 41) 

sin^2^ = 4vr^_ 

This gives s2=0.23117(33), c=0.87683(19). Such a 9, which by definition does not depend 
on rrit and , is used to determine the so-called a-Born approximations for electroweak 
observables. 

The bare gauge couplings q;°, and the bare mass M° are expressed in terms of 
a, Gn and M^, 1/e and where ji is the 't Hooft's scale parameter and 2e — A — D in 
dimensional regularization scheme. The W mass, M^, is treated on an equal footing with 
the other observables. 

The a-Born approximation automatically includes purely electromagnetic corrections. 
In terms of s and c, the expressions for the hadron-free pseudo-observables are very simple 
in that approximation 

(M^/MJ^ = c , (43) 
{9{r = if , (44) 
{9l/9if = l-MQf\s'. (45) 

The a-Born approximation with the due account of the final state QED and QCD radia- 
tion factors gives simple expressions for the observables in the hadronic decays of Z boson 
as well. 

The electroweak corrections for all observables are calculated in the framework of 
LEPTOP in terms of a^, M°, m^, m^^. and (1/e, fi) in one-loop approximation. In this 
approximation m° and M° can be replaced by the physical masses rrit and M^. 

By expressing a^, a^, M° in terms of a, G^, and (1/e, n) one derives formulae 
in which terms (1/e, //) cancel out. Thus introduction of counterterms is avoided. 

The explicit analytical expressions for the LEPTOP electroweak corrections in terms of 
G^, a, c, s, t = [mt/M^Y' and h = (M^/M^)^ are given in section 4.1. 

The flowchart for the LEPTOP approach is shown in Fig. 9. 
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FLOWCHART OF LEPTOP 



Select the three most accurate observables: 
G^,M,,a{M,) = a 



Define angle 9 ( s = sin 6*, c = cos 9) 
in terms of G^, and a: 

G, = {7,/V2)a/s\^Ml 



Define the Born approximation for other electroweak 
observables in terms of G^, and 9. 



Introduce bare couplings {a ,a^,a^), masses (M^,M^, 
(including m^)) and VEV i] in the framework of MSM. 



Express a^, a^, Af^ in terms of G^, and a in the one-loop 
approximation, using dimensional regularization {l/e,/j,) 



Express one-loop corrections to all other electroweak observables 
in terms of a^, a^, M^, rrit, and hence in terms of G^, 

M^, a, rrit, M^. Check cancellation of 1/e and ji. 



Introduce gluonic corrections into quark loops and QED and QCD 
final state interactions for hadronic decays 

(in terms of a,ds{M^),mi,{M^),'mt). 



Compare the Born results and Born + one loop results with 
experimental data on Z-decays and M^. 



Make a global fit for three parameters rrit, M^^, ds{M^). 



Predict the central values of all electroweak observables 
and their uncertainties. 



Figure 9: LEPTOP flowchart 
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1.8.3 TOPAZO 



The TOPAZO code uses the MS (modified minimal subtraction) scheme for all types of 
interactions, including the electroweak ones, as introduced in Its approach is quite 
different both from those of BHM/WOH, ZFITTER, and LEPTOP. Its main steps are presented 
in the following items: 

• The MSM Lagrangian is assumed. Excluding fermion masses and mixing angles the 
MSM is a three-parameter theory: i.e. g^, M° and sin^o- At the level of defining a 
renormalization scheme, two essentially different approaches may be distinguished: 
to prescribe precisely what a parameter of the Lagrangian is or to prescribe precisely 
what a counterterm is. It is a matter of convention, since only their combined effect 
appears in the confrontation with the data, and once the latter is chosen it would 
then seem natural to follow the consensus with respect to QCD. The bare parameters 
are fixed by considering three data points — a(0), G^, and M^. These quantities 
are computed up to one-loop diagrams and fitting equations are written: 

p° = /-,(«,G^,M„A) , ^ = l,2,3, (46) 

where p° are the bare parameters and A = —Ijin — 4) + 7^; — Invr. Everything is 
worked out explicitly in the 't Hooft-Feynman gauge. 

• The fitting equations are solved perturbatively (order- by-order renormalization), but 
gauge invariant (fermionic) higher order leading terms are always re-summed. 

Strictly speaking, one usually solves an implicit equation /(a. A) = 0, where a is 
some parameter and A is a coupling constant, by expanding a around oq, the solution 
of /(do, 0) = 0. However, one could just as well expand a around a, the solution 
of some other equation h{a, A) = 0, as long as a is gauge invariant. Thus in our 
implementation a(0) will evolve into a(M^) and the lowest-order approximation for 
= sin^ 9o is 



1 
2 

1 + 



4 7ra(M^) 



(47) 



z 

2v^7r2 



X {UV finite combination of two-point functions}, (48) 



where (p^) ^ — 1 is determined by an ultraviolet finite combination of two-point 
functions. The detailed definition of is left for the section 4.2 — Eqs. ( p27| ) 
and (p^-p3D. 

Eq.(|47|) is an algebraic solution of the following relation^ 



G,, 



(49) 



^For analogous relations for BHM/WDH, ZFITTER see Eq. (p§) and for LEPTOP - Eq.fil 
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If no resummation at all is performed in TOPAZO then = 1 and, in this limit, 
one recovers the weak mixing angle as defined by LEPTOP. The parameter is 
strictly connected to the Z wave function renormalization and usually all bosonic 
contributions are expanded up to first order, in this case the p parameter being 
instead denoted by p^, Eq. (|227|) , but the possibility of a resummation for a certain 



gauge invariant subset in the MS framework is built in. It should be noticed that 
p^ or its variant p^ represent the natural extension of the Veltman's p-parameter, 
as defined at low energy [Q, to the scale M^. Its asymptotic behavior, for large 



mt, is exactly the same as dictated by Ap — see Eq. (q^)- That is why all the 
higher-order leading [m^ corrections will simply be added to p^Ap^)- 



By a proper redefinition of the bare coupling [^, corresponding to a gauge field 
re-diagonalization in the quadratic part of the Lagrangian, in the ^ = 1 (t'Hooft- 
Feynman) gauge the following properties are fulfilled: the sum of all Zjj vertices, 
X]t,{vertex}, is ultraviolet finite and the Z — 7 transition satisfies ^^^(O) = 0. In the 
renormalization procedure of TOPAZO it is observed that, unlike QED, no one-to-one 
correspondence exists between the bare parameters and experimental data points. 
Therefore no attempt is made to give an all-orders relation as for , but rather 
it is observed that the mixing angle defines the distribution of the vector current 
between Z and 7. It is an accident of the minimal Higgs system (or more generally 
of representations where the Higgs scalars only occur in doublets) that the vector 
boson masses are not both free and beyond lowest order different definitions of the 
weak mixing angle receive different radiative corrections. Also, the W mass has 
no special role in TOPAZO, being computed like any other quantity, once the bare 
parameters are substituted through their explicit expressions. Typically, 



sin^6'o = + {C(a)UV finite bosonic corrections} 

+ {0{q) A-dependent terms} , (50) 

and in computing for instance Z ^ ff the latter cancel their A dependence against 
the Z — 7 transition, leading to an ultraviolet finite result. According to the chosen 
strategy, all bosonic terms are expanded or, alternatively, a gauge invariant sub-set 
of the bosonic corrections is re-summed in (in the MS environment). In the 
second case a (usually small) remainder is left, which compensates against vertex 
corrections. The corrected Z propagator, leading to mass renormalization and to a 
Z wave-function renormalization factor, is ultraviolet finite by inspection. 

The processes Z ff, e+e~ ff are computed up to one-loop diagrams with 
lowest-order (partially re-summed) parameters plus tree diagrams with one-loop 
corrected parameters. Ultraviolet finiteness is checked both analytically and nu- 
merically (independence of A). Roughly speaking, the central objects in TOPAZO 
are the S-matrix elements for a given process, so particular care is devoted to the 
proper treatment of the wave-function renormalization factors (see Ref. |55||). 
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1.9 Phenomenology of Z boson decays 



The decay width of Z — > // is described by the following equation used by all five codes: 

= r(Z ^ //) = ^NlV,{{glfRl + {g{fR{\ , (51) 

where gl^ and (7^ are effective electroweak couplings defined in section 1.4.4, Nl = 1 or 3 
for leptons or quarks (/ = /, g). The factors i?^ and i?-^ describe the final state QED and 
QCD interactions and take into account the mass mj. The radiation factors are universal, 
to a large extent, and will be described below in this subsection. 

The standard width Fq in Eq . (|5TD is known with very high accuracy: 



G 

To = ^ J" = 82.945(12) MeV 
24v 27r 



(52) 



For the decay into leptons // the radiation factors i?'^^ are especially simple. For 
charged leptons: 



4Fn 



a 



(53) 



The mass term in Eq. (|53D is negligible for I = e, /i and is barely visible only for I 
T {ml/Ml = 3.8 X 10-^). 
For the neutrino decay: 



F, = 8((7^)2Fo, 



9.. = 9, 



(54) 



For the decays into quarks qq the radiation factors are given by 



1 + 61 + ^val - 12.76706a^ + 12^^a,5,^ , 



(55) 



1 + 51 + 



5„ - 2/p) J(2) 



+ 



-12.76706 - 2JFx(3) 



2 s^am ' 

mt 



(56) 



where 



3 1 

/ 44 2 s 
y675 135 J m'f 



5^ = 1.40923 + I — - — In ■ 2 , 2 ' 



26 



5lr 




II 54 m? ' 



(57) 



and 



as{s) 



TT 



71 



X(3)( 



37 7 

-— + / + — a; + 0.0132 
LZ oi 



5651 8 23 



x] 



216 



36 



2 ^ 67, 23, 
vr^ + C 3 + — + — 
^ 18 12 



(5^ 



where I = \iax, rrit = mt(poie) (and s = for Z decay). 

Note that we have automatically absorbed the finite mass terms into the QCD correc- 
tion factors and that in the large logarithms ln(m^/s) have been absorbed through 
the use of running parameters 



m 



m{m ) exp 



a,(m2) fj{x) 



mem 



4 



1 + ^as(m) + Kal{m) 



(59) 



where m = rripoie and Kf, ^ 12.4, ~ 13.3. The various codes may differ in the construc- 
tion of g-^ and g-^ but, a general consensus has been reached for final state radiation. In 
this way we obtain the hadronic and total Z width as 



^ r{z^ff) + 4roNiR 

f=udcsb 

2 

.3 



-0.41318 (T.V'^ 

f=efj.T 



a 



s ' 



(60) 



Actually, the hadronic width and the partial qq widths deserve some partial comment, 
while the complete theoretical framework is described in The singlet QCD con- 



tribution which is simple and unambiguous for the hadronic width becomes ambiguous, 
starting at (9(q;^), for individual qq channels. In fact, some sort of agreement has recently 
been reached on these matters but we want to summarize the roots of the problem. From 
a pragmatic point of view there is a hierarchical description where 
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V, = V[Z ^qq{g) + q'q'] (61) 

for all g', such that rriqi < ruq. However, we could have a democratic description where 
the final states qq + q'q' are assigned for | to Tq and for the other | to Vqi. The two 
descriptions agree — fortunately for the leading terms. In general, one could also decide 
that such final states should not be assigned to any specific channels in such a way that 

^ r, ^ r, = ^ r(z ^ gg) + ^ r(z ^ gggV) + ■■■ (62) 

<? 9 q,q' 

In particular the 0{a^^) contribution to Vy cannot be assigned to any specific fiavour. In 
our approach the T^^'^^ corrections terms are included according to isospin, i.e. propor- 
tional to I^^ , and therefore cancel in the hadronic width while summing over the first four 
flavours m, c?, c, s. However, it must be pointed out that no general consensus has been 
reached on this particular point, so that the prediction for F^, F^, Fc may differ from one 
code to the other, although the total hadronic width will remain the same. 

1.10 Electroweak corrections: Basic notions 

Above, we introduced a universal language describing the general features of radiative 
corrections, the language of effective couplings. At a secondary level, each computational 
scheme, being a particular representation of the same general concept, may use different 
building blocks, which, however, are so deeply related that one could even attempt the 
construction of some sort of 'Rosetta stone'. Since one of our main motivations was to 
build up a theoretical framework that includes some estimate of its own uncertainty, we 
must necessarily spend some time in discussing some of the specific building blocks. 

Indeed, even though there is a high degree of universality in the realization of the 
theoretical uncertainty, in practice the way in which this realization becomes effective is 
strictly related to the actual implementation of higher-order radiative corrections into the 
various codes. In order to discuss the numerous effects involved in the calculations we 
have to introduce additional notions. 

1.10.1 Compcirison of notations of different codes 

As explained earlier, all electroweak loop effects in Z boson decays are concealed in the 
effective couphngs g^^ and g"^. The treatment of these couplings, unlike that of the radia- 
tion factors -R^^, differs from one code to another (having several different realizations). 
In order to concentrate on these realizations let us for a moment forget about radiative 
factors . Then 

V,A 

^f-^NlV,[{gif + {g{f]. (63) 

The five codes deal with three main realizations. 

In BHM/WOH , ZFITTER the following notations are introduced: 
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(65) 



where 



Kf = 1 + SUf 

SO that Eq. (|63D may be rewritten in the form 



(66) 



= ToNipj[A{lj - 2QfslKff + 1] . (67) 

Here the electroweak corrections affect 5p/, 5k f and s^^. They will be discussed in the 
next subsection. 
In TOPAZO, 



= 4 To iV/p, \(lf -2Qff + 5gl)\ {if + 5g 



(68) 



where is defined in (^), in ( |227| ) and 6g^,6g^ — see ( pT]) — contain (part of) the 
bosonic corrections as well as the vertex corrections. 
In LEPTOP, 



9 



f 



r(3) 



1 + 



9l 

9{ 



R ' 



(69) 



(70) 



and Va,r are simple functions of {rrit/M^Y and (M^/M^)^, as described in section 4.1. 



1.10.2 Basic notions of different codes 

Here we will concentrate on a specific realization of the language of effective couplings in 
different codes. In particular it is important to clarify the notion of leading, non-leading 
and eventually of remainder terms in BHM/WOH, TOPAZO and ZFITTER. 

We start our presentation by considering two well known objects, Ar (which is used 
by BHM , TOPAZO , ZFITTER) and the partial Z width. There are many notions which are 
common to both these objects. 

The quantity Ar is nothing but the effective coupling of the /i-decay, when it is being 
implemented with higher order corrections, see the couple of Eqs. (^) and (0). In this 
case there is only one effective coupling, as the p decay is a purely weak process mediated 
by a charged current. For the partial Z width, the situation is more involved, since we 
have a decay which is mediated by the neutral current. That is why it is described by two 
effective electroweak couplings. The introduction of Ar will be very useful in clarifying the 
notion of leading and of remainder contributions to radiative corrections. However, from 
a general point of view this is only one of the many specific realizations of the effective 
coupling in the p decay. 
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As an example of the leading-remainder splitting we consider and Ar 
and 5pf (0), and Kf and 5k f (|65| , |6^ ) — the electroweak corrections to Ar and to F/. 
Next we subdivide pc as introduced in Eq. ( P^D and p/ and ft/ of (0, into a leading 
term, A,, and remainder. 



Arem, tcrms as follows: 



P/ 




(71) 



1 — 5pf 1 — Ap — Ap/ 



(72) 



^2 



Kf = l + = (l + Aft/,rcm) (^1 + ^Ap^j . (73) 



From Eqs. ([71|-[73|) one may notice that each coupling contains some universal, flavour- 
independent piece (Aa, Ar^, Ap, Ap-^), and some flavour- dependent remainder. The 
former comprise the leading terms (Aa, Ap) they are re-summedto all orders in accordance 
with the renormalization group equation for Aa and with the proper inclusion of leading 
irreducible terms for Ap The latter are normally small, since all potentially large 
contributions have already been subtracted and shifted to the leading terms. For this 
reason, they are placed freely either into the denominators, ( |7T| , (|72D, or as a factor of the 



leading contribution, (73); in doing so we follow Ref. ^^l- case however, namely 

Eq. (^), there exist arguments in favour of keeping the remainder in the denominator. 
In Ref. it was proved that the right-hand side of Eq. (|7lD actually means that, on the 
two-loop level at least, all fermionic mass singularities are located exclusively in Aa. This 
may be argued by a partial expansion of a simplified case (Ap = 0) up to the two-loop 
order terms, 

- — X X ~ T. — ir^ 1 + 1 — ^ ' (74) 



Aa - Arrcm (1 - Aa) V 1 - Aa 
from which one can see that the scale of Arrem is actually a{M^] 



One should emphasize that since there is no Aa in Eqs. ([72D, [73|) and since the non- 



leading terms of the order A Aj-em are unknown, no arguments in favour of putting Aj. 



as it is done in (|7^, |73D can be presented. The size of these uncontrolled terms should be 
treated as an intrinsic theoretical error. 

All the remainder terms have a similar structure: 



Ar,em = Ar^i°°P'" + Ar2i°°P'""= + ^Ap^ - Aa , (75) 



s 



Ap/,rcm = Ap"°°P'" + Ap2i°°P'"-=-Ap, (76) 



Akj,,^^ = Aft}'°°P'" + Afi;21oop,aae _ -^^Ap^ . (77) 
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They contain all the terms known at present: the complete one-loop 0{a) corrections 
(two-, three-, four-point functions) and complete two-loop 0{aas) insertions to two-point 
functions, from which the leading 0{a) and 0{aas) terms are subtracted: 



(X) 



3a 

w w 



M2 




(7J 



where braces in ([T8|) and below mean that this expression describes simultaneously both 
quantities Ap and Ap^, which appeared in Eqs. (|75D -(|77D. 

The term X in ( fTHf ) is a next-to-leading order term, whose proper treatment is rather 
important: 



X = ne 



M2 



M2 



Hoop 



MS 



Ap" 



(79) 



In Ref. it was argued that in the on-mass-shell renormalization scheme this term 
should be re-summed together with Ap", since in the MS scheme it appears to be auto- 
matically incorporated. Similar arguments in favour of such a resummation were presented 
in Ref. . In X the UV divergences are removed according to the MS renormalization 
scheme with p = M^. The separation of X is not unique; it makes the resummation 
dependent on the renormalization procedure. 

The leading contribution Ap is built out of the same terms as ([THD. But they are 
normalized by rather than by a/s^M^, as is required by the resummation proposed 
in Ref. ||: 



Ap 



{X) 



where 



Ap" + Ap" + Ap""= + Ap""= + (X) 



NcXt 



1 + xtAp 



(2) 



Ml. 



C2 



TT 



+ (X) , (80) 
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W~W FR 1 



(81) 
(82) 



and Xt is the Veltman heavy top factor ||5T|. The coefficients c\ and C2 describe the first- 
and second-order QCD corrections for the leading xt contribution to Ap, calculated in 



Refs. [p^ and |53]. Correspondingly: 



Cl 
C2 




-TT^ (2.155165- 0.180981 n/) . 
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The function Ap(^)(m^/M^) describes the leading second order Xt (two-loop electroweak) 
correction to Ap, calculated first in the = approximation in Ref. [|54| and later in 
Ref. for an arbitrary relation between and rrit. 

The partial decay width Z ^ bb contains an additional rrit dependence due to vertex 
diagrams (see Fig. 10). 




Figure 10: Top quark exchange diagrams which contribute to F;, (unitary gauge) 

As a consequence, the effective couplings pf, and Kf, contain additional leading terms of 
the order 0{G^m^). The complete one-loop approximation for the Z ^ bb partial width 
was calculated in Ref. |Q . We first redefine remainder terms by an additional subtraction 
of the leading one-loop term originating from these diagrams: 



where 



,rom 

Apb 



Apfe^rem " 2Apfe 
A/«;,,rcm + Apb , 



a 



87rs2 M2 
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(86) 



Following papers |55| and |5^, the two-loop order QCD and electroweak leading terms 
in the Zbb vertex are implemented by an additional re-definition of effective couplings pb 
and Kb. 



Pb 

Kb 



P6(i + n) 
Kb 

i + n ' 



^7) 



where Tb is given by 



Tb = -2xt 



l-^a,(m,) + x*r(2) 



(89) 



A compact analytic representation for the two-loop functions p*-^^ and r*^^) -^^^g q\^q given 
in Ref 0. 

We have just discussed what is known in the literature about the treatment of the cor- 
rections of the order 0{G^j_mf). The Zbb vertex also contains a logarithmically enhanced 
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term, 0[aln{mf /M'^ )], whose contribution is comparable to the leading one. Recently, 
QCD corrections were also calculated for this term, (see contribution by Kwiatkowski and 
Steinhauser in the QCD Part of this Report). This correction, however, can be nearly 
completely absorbed into the final-state QCD corrections. What remains is approximately 
one hundred times less than the QCD correction for the leading term as given in (|89|) . 

What has been presented so far in this subsection, is to an extent inherent in all 
five codes. However, in order to reach a better understanding of our set of theoretical 
predictions and of the strategy adopted to extract their related intrinsic errors we have 
to devote more time to discussing other realizations of the effective coupling language. In 
TOPAZO — see, Eq.(B3) — the partial Z-widths reads: 



4roiv// 



lf-2Q,s^ + 6gi)\(lf+6gf:' 



(90) 



where s^, defined by (0) and ( |229| ) (see also Ref. [0), and is to be put in partial cor- 
respondence with the s'^Kf of Eq. (|65|) . Moreover, 6g-^,6g-^ contain bosonic corrections 
as well as vertex corrections. Again, the basic point under examination is the leading- 
remainder splitting. The main idea beyond this realization is to write a system of equa- 
tions that connect the MS parameters of the theory in terms of a,G^ and M^. These 
equations contain the effects of radiative corrections up to a certain order in perturbation 
theory and must be solved consistently and respecting gauge invariance. As a result of 
this procedure we end up in Eq. (|90|) with a defined according to Eq. ( |227| ). This 
parameter, which properly re-sums the gauge invariant fermionic corrections, contains all 
isospin breaking terms and includes, beyond 0{a), all the presently known higher-order 
terms — the 0{a^) |58l,C(aaJ ^ and O^aal) H corrections. 



The bare weak mixing angle, according to the previous strategy, will always be ex- 
panded around s, Eq.( p29|) , where from Eq. ( p32|) and (|233|) is used. These two 



quantities (p^ and p^) may differ whenever a resummation of bosonic corrections is per- 



formed for the weak mixing angle. Once the divergent terms are treated in the MS 
framework, (the leading term) will receive a correction (the remainder) As^. At this 
point, when the Z wave function renormalization factor 11^ is included, we usually ex- 
pand all the remaining corrections up to first order — we decompose II^ into its leading, 
fermionic part and a remainder AII^ (see Eq. ( |231[ )). This AII^ is always expanded. 
After inclusion of vertices and fermion wave-function renormalization factors we end up 
with 
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[2Ff - '^VfAUz 
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'2F[ - \lfAVL 
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- 2QfAs^ 



(91) 



where is again defined as (|4l|) and F^^ are the flavour- dependent vertex corrections. 
Additional versions of this realization will be discussed in the framework of the theoretical 
options. Here we need only mention some of the problems connected to the resummation 
of bosonic contributions in this realization (as well as in others). It is strictly related to 
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the question of gauge invariance. As it is well known, the vertices and the bosonic parts of 
the vector boson self-energies are not separately gauge invariant. On the one hand there 
is no rigorous procedure for re-summing the vertices, and on the other hand any attempt 
to isolate the gauge variant parts and to throw them away is not unambiguously defined. 
Even when we identify the universal contributions from vertices and boxes (for instance 
by working in the i?^ gauge), to be combined with the vector boson self-energies in order 
to get the one-loop gauge invariant dressed propagators, we still have the freedom to shift 
from one part (to be re-summed) to the other some process- independent, ^-independent, 
UV finite function of with the proper asymptotic behavior. Strict enforcement of 
gauge invariance is one of the roots of Eq. ([9T|), in that AII^ is not absorbed into the 

parameter and that particular care must be devoted to the proper definition of the 
remainder, As^, when resummation of irreducible terms is performed. This fact accounts 
for one of the many structural differences with other realizations. 

In LEPTOP, due to the proper choice of the a-Born approximation, the electroweak 
radiative corrections turned out to be small — for around 175 GeV. The smallness 
of corrections is a result of the mutual cancellation of different equally important terms. 
Therefore, functions V/^ representing loop corrections are not subdivided into leading 
and remainder terms and no resummation is performed. 
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2 Options, theoretical uncertainties 



During the completion of this work it has become increasingly evident that there is a need 
to quantify the effect of our partial lack of knowledge of the missing higher-order terms in 
radiative corrections. Thus we have introduced the notion of option, which refers to a set 
of possible and plausible alternative implementations of the full machinery of radiative 
corrections within a given renormalization scheme. In order to explain the treatment 
of higher-order terms and the interplay between pure electroweak and QCD corrections 
we must once more remember that the effective coupling description of the Z width has 
several different realizations. 

Quite independent of specific details, all the realizations single out two main compo- 
nents in each observable: 

= 0b + ^0 . (92) 

The term Ob, giving the bulk of the answer, is often called the Born approximation (the a- 
Born approximation in LEPTOP, or improved Born approximation in BHM/WOH, TOPAZO and 
ZFITTER), or the leading contribution to O. The term AO represents a small perturba- 
tive correction, often called remainder or non-leading contribution. Different realizations 
usually have different ways of performing this splitting so that, while they agree at the 
0{a), there are differences which start at 0{a'^). 

Independent of the particular realization of g^jQ^, these effective couplings are com- 
plex valued functions, due to the imaginary parts of the diagrams. This, however, will 
have some relevance only for realistic distributions while for pseudo-observables they are 
taken to be strictly real. Given that the universal language has in practice many al- 
ternative realizations and that the actual implementation of the options is very much 
code-dependent, we have not created a common set-up for the electroweak options. From 
this simple fact follows the need to discuss at length the physical motivations behind our 
options. Although we have tried, as much as possible, to illustrate them from a general 
perspective, it is also evident that some space has been left to analyze certain specific 
issues. A description of the implementation of the options in the various codes, with the 
relative flags, will be given in section ^. 

2.1 Factorization of QCD Corrections 

No matter which realization we are using, one problem naturally emerges when final- 
state QCD corrections are switched on. This problem is connected with the folding of the 
non-corrected widths using the QCD factors Ry^- One can take the point of view that 
non-universal and flavour-dependent couplings should also be folded or, to the contrary, 
that only the universal effective couplings should multiply the QCD radiation terms. This 
option merely reflects our ignorance of the mixed 0{aag) corrections, with the noticeable 
exception of the bb partial width where the leading rrit part of these corrections is actually 



known |5^. If we consider the bb partial width we can write 



{glf + {If 
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.{Z hh) 



At this point naive factorization would imply 



1 -4xi 



{9lY + {9lY 



(93) 
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while the computed FTJR term gives 
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(94) 



(95) 



so that the correct QCD coefficient in front of the heavy top factor turns out to be 
—2.290 tts/vr instead of a^/n. Of course, an approximate but factorized expression is still 
possible and can be written as 
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(96) 



where the uncertainty has now moved to order o?^Gyjn'i. Actually, we have not yet 
specified the scale of in the previous equations but these corrections have been im- 
plemented such that the universal QCD factor is computed as 1 + a^iM^^jn^ while the 
specific FTJR term is computed with agirrif). As a consequence the 1 + ^^(M^)/?? factor 
is not included for the 6— quark asymmetries and for the effective 6— quark mixing angle. 
In general, however, the 0{aa^}j corrections for the Zff vertex is presently not known 
and we must accept an intrinsic uncertainty associated with the two procedures, i.e. fac- 
torization or non-factorization of electroweak and QCD corrections. Once again, if Y\ is 
the (improved) Born qq partial width and 



EW 



n (1 



+ A« 

' EW 



then 
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(97) 



(9^ 



is roughly assumed as the corresponding uncertainty. This type of uncertainty can be 
illustrated very well by asking how correct it is to shrink an electroweak blob to a point 
before allowing for QCD radiation. 



2.2 Genuinely Weak Uncertainties 

In this section we briefly discuss the main ingredients which enter the pure weak cor- 
rections to the pseudo-observables — the resummation of the one-particle irreducible 
vector boson self-energies, the scale in vertex corrections and the linearization of the cor- 
responding S-matrix elements. We have avoided in this chapter any intensive usage of 
the lengthy formulae introduced in the first part of this Report and summarized in some 
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detail in section |^. Indeed, there is a simple set of problems in the implementation of 
radiative corrections which is inherent to perturbation theory and does not depend on any 
specific approach. Only formulae belonging to different realizations represent the tech- 
nical transcription of these implementations. Real progress is usually achieved whenever 
a new term becomes available, otherwise we are left with heuristic arguments or with 
ingenious attempts to improve upon perturbation theory. Different paths along this road 
represent our present degree of inaccuracy and to understand these differences already 
gives a hint on how to proceed. The structural and logical essence of these differences can 
be explained without a massive use of equations. 



2.2.1 Leading— Remainder splitting 

Before we move to discuss the next source of uncertainty we must recall again that, 
generally speaking, the effective couplings g-^, g-^ contain a leading and usually re-summed 
part and a non-leading (remainder) one, quite independent of the specific realization. For 
instance, in one of the realizations the building blocks are pf and Kf, while in another they 
are p^, and 6g-^, 6g-^, and in a third they are s^, Va_ and Vr. The way in which the non- 
leading terms can be treated and the exact form of the leading-remainder splitting give 
rise to several possible options in the actual implementation of radiative corrections that 
in turn become another source of theoretical uncertainty. For instance, for all the objects 
Ar, Apf and Akj, we can introduce the decomposition into leading and remainder. Since 
we know how to proceed with all objects in the leading approximation the only ambiguity 
is due to the treatment of the remainders. Clearly, after the splitting Ar = Ar^ + Arrcm 
there are in principle several possible ways of handling the remainder^: 



Ar 



1 - Ar^ - Arrcm 



l-Ar, 



1-Ar, 



1 + Arrc: 

1-Ar, 



-Arv 



1 _|_ ^'•rc 
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l-Ar^ ' "'rem 



(99) 

Actually, these options differ among themselves but the difference can be related to the 
choice of the scale in the remainder term. A complete evaluation of the sub-leading 
(9(G'^M^m^) would greatly reduce the associated uncertainty. At the moment these sub- 
leading terms are simulated just by varying the scale in the remainder. Typical choices are 
described in section |^ and may vary from realization to realization. A different approach 
to the problem H] will be illustrated in section |2.2.5. 



^In case of Ar, the last two expansions are not valid indeed, see a discussion around Eq. 
present them here nevertheless for the sake of illustration. 



Q. We 
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2.2.2 Scale in vertex corrections 



Another possible option, which is realized by some of the codes (but not all), has to do 
with the scale of a in the non-leading corrections, in particular the vertex corrections. To 
make this point clear we use one of the realizations of the effective couplings and analyse 
the Z ^bb decay in some details. The essential ingredient will be Apf, which we write as 



Spbirrit) 



Apd + Spb{mt) 

1 + + + O («ln m^) + Sp^ 



JVL 
b 



(100) 



where 6pb is a correction specific to the bb channel. The question naturally arises as to 
what to use for the scale in S and in the sub-leading logarithmic term — G^, a{M^) or 
q;(0). Obviously the same kind of option will be present in the vertex corrections for light 
fermions where the expansion parameter is formally a / (47r sin^ 6 cos^ 6) . Different ways to 
interpret a in these expressions give rise to different results, from Q;(0)/(47r sin^ ^ cos^ 6*) 
to G^M2/(2v^7r2). The difference between possible identifications of coupling constants 
in the 0{a) corrections represents, of course, effects of 0{a'^). The fact that the neutral 
current amplitude is automatically expressed in terms of Gfj,M'^ is a possible heuristic 
argument to adopt the same strategy in the evaluation of the presently known 0{a) 
corrections, but in order to be on the safe side, the differences should be considered as a 
theoretical uncertainty, at least according to many authors. For instance, we know from 
a specific calculation that for inclusive quantities the final state QED radiation is 
controlled by a(M^) and not by a(0). The option of variation of the scale of non-leading 
corrections has been implemented by the majority of the codes. 



2.2.3 Linearization 

Another example to be discussed is the following. In almost any realization we have 
the possibility of using an expanded versus a non- expanded option — to linearize our 
expressions. To give a partial illustration of this possibility we use the realization of 
Ref. IMI where we have 



sin' 




+ \ Sg[^ + 1 (4 — 1) Sg'- expanded 



non-expanded 



On the same footing, the 5g^^,5gf^ terms in each partial width will be expanded up to 
first order in the expanded option. More generally, the difference between the two options 
in the evaluation of O^, where O is given by Eq. (p2D, is equal to (AO)^, a two-loop 
reducible but non-leading contribution, and from this point of view we clearly do not need 
a specific exemplification. By comparing the two options we obtain a rough estimate of 
the importance of the missing non-leading two-loop effects. It goes almost without saying 
that when this option is implemented on top of the leading-remainder splitting and of the 
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choice versus a(0) we can have rather different behaviours in the global theoretical 
error according to the size of the remainder. Thus whenever we need to quantify the 
effect then the remainder must be considered in detail and the effect becomes realization- 
dependent. This has some relevance for the de-convoluted leptonic forward-backward 
asymmetry, where the leading term is small by itself due to accidental cancellations. To 
summarize we can say that the theoretical error on is very sensitive to the sum of 

•1 FB 

various factors: 

• expansion versus non- expansion 

• definition of the leading part of sin^ 

• scale of a for the non-leading terms. 

A further exemplification of what we call option for radiative corrections, in part 
connected to a selection of the scale in the coupling but generalizable to all parameters 
appearing in the remainder, can be illustrated as follows. Suppose that a given quantity 
A, function of the parameter a, is given in perturbation theory by the following expansion: 



A 



a + g\a^ + /i(a)J + g'^ 

= -a + gh{a) + 0{g^) 
a = a/(l — ga) , 



/2(a) +0{g' 



(101) 



and that only the fi term is actually known. It could be decided that d is the effective 
expansion parameter (or that in the full expression we change variable a ^ d), which of 
course in the truncated expansion introduces the option. 



= a + ghid), (102) 

so that AA = g'^a^f[{a) would be our estimate of the associated theoretical uncertainty. 
Sometimes it is only a rough estimate, since there is no guarantee that the irreducible 
terms (7^/2 are essentially small in size. 

2.2.4 Resummation 

Another source of theoretical uncertainty is connected with the treatment of the physical 
Higgs contribution. As we know for large Higgs masses (M^ ^ M^) there is a correction 
term in Ap which is only logarithmic, in contrast to the heavy fermion case, the so-called 



Veltman's screening effect [^. With respect to this correction Ap^ we can: 



expand it to first order in a as is sometimes done for all bosonic corrections; 

re-sum the leading part of it, Ap^, for relatively large values of M^, e.g. > 
M^exp5/12; 
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• re-sum in p the whole physical Higgs contribution. This requires a further comment, 
since this term is not UV finite by itself and therefore the resummation procedure 
must be understood strictly in the MS scheme with a scale /i = M^. 

As an additional and rather general comment, which somehow collects many of the previ- 
ous considerations, we stress once more that there are different ways of implementing the 
resummation of the vector boson self-energies. These choices, which in turn are deeply 
related to the proper definition of remainder, differ from code to code — at least in their 
default settings. We have already illustrated the splitting Ar^ — Arrom and simply add a 
few additional considerations. Resummation is very often the main recipe for separating 
a small remainder from the bulk of the effect. 



One choice consists in a resummation which includes the square of the Z — 7 mixing- 
term ^ with the option of strictly keeping in the resummation only the one-loop 
irreducible terms. The default thus corresponds to an additional mass counter- 
term which enters the field and coupling renormalization constants and modifies 
the quantity Ar. The resummation of the modified Ar leads automatically to the 
factorization property which takes into account the proper summation of all the 
leading higher-order reducible terms. 

Even more generally we can distinguish among complete expansion of the one-loop 
self-energies, partial resummation of fermionic self-energies or partial inclusion of 
bosonic self-energies in the resummation. There are two considerations to be made 
at this point. Sometimes accidental cancellations occur among the fermionic and 
the bosonic sectors, which would suggest a similar treatment for both; however, the 
bosonic sector is not gauge invariant by itself. Thus any resummation of bosonic 
parts must properly identify some numerically relevant but gauge invariant sub-set. 



For completeness we recall that one of these identifications gives in the ^ = 1 gauge |60 



S.Ap') = S,AP%=i-^e^'-^{p'-Ml)l{p'), (103) 

Sill u 



where I{p ) = —Bo{p ,M'^,M )/(167r ), Bq is a scalar formfactor and S denotes a 
possible identification of the gauge invariant part of S. 

To spend an additional word on electroweak uncertainties and to understand the im- 
plications of some of the procedures used as a possible estimator of the theoretical error 
we consider a fictitious quantity X defined as X = 1 — 4 sin^^^g. Each code will define 
some effective Born approximation to sin^ 6'^^ and here we distinguish between three ba- 
sic possibilities: complete expansion of the corrections (E), a resummation which only 
includes the fermionic self-energies (FR), or a global resummation (R). We refer here to 
the complete set of formulae given in section § and simply quote the adopted strategies: 
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• BHM/WOH adopted as the default of the resummation of the entire set of self-energies, 
including the Z — 7 mixing term; 

• LEPTOP, which expands all those contributions not re-absorbed in the running of a 
into a; 

• TOPAZO's default, which re-sums in the MS framework the term of Eq. (|233|) , 
while an option is left in which only its fermionic content is re-summed; 

• ZFITTER's default, which also re-sums in the MS framework the X term of Eq. (|82[) , 
allowing asr options the resummations of the leading terms only. 

Let us define Xq as the effective Born approximation of X and define a non-leading 
(remainder) part, X = Xq + [a{M^) /tt]Xi + 0{a^). The three different procedures, E 
or PR or R, will find Xq = 0.07528, 0.08672, 0.07184 (for = 175 GeV, M^ = 300 GeV 
and = 0.125). In fact, there is no unique way for global resummation (R), so that 
the last number can vary a little (say, from 0.07184 to 0.07416). Using the corresponding 
value for sin^ we get Xi = —1.344, —5.976, —0.016. Whenever we compute and use 
the square of the remainder to estimate part of the uncertainty (there are other options 
around) we are bound to see rather different behaviors, depending on the adopted leading- 
remainder splitting. It should be clear that the remainder itself is subject to several 
possible options, from the scale of the coupling to the choice of those terms which are to 
be considered small and perturbative. These options themselves contribute, sometimes 
sizably, to the final uncertainty, quite independently from linearization (expansion) and 
in the end every adopted procedure is somehow equivalent to a proper choice of the scale. 
Thus a preferred set of options is equivalent, in some sense, to an ideal optimization of 
the perturbative expansion. 

2.2.5 Estimate of the missing terms in higher orders 

While some of the realizations of the effective coupling description make use of a (partial) 
resummation of higher-order terms, thus trying to improve upon ordinary perturbation 
theory, we have other realizations where the mixing angle is defined only in terms of 
G^, and of the running constant a{M^) while everything else is strictly expanded. 
Thus, a Born approximation is defined in terms of an s^, the definition of which was given 
in Eq. ( |4^ ) and the realization is constructed in terms of the one-loop approximation 
with respect to the genuinely electroweak interactions. As has became clear from the 
previous discussion, there is no common set of implemented procedures for describing 
the theoretical uncertainties but only some rather general guidelines. Certainly, when we 
move to a concrete implementation, it happens that every code has its own internal ways 
of estimating the missing higher-order effects. Codes which do not foresee, for theoretical 
reasons, the possibility of expanding the remainders with respect to the leading terms 
or of playing with the choice of the scale in the remainders, introduce another set of 
options based on a different estimate of the not-yet-calculated diagrams or terms in a 
given diagram. The basic and alternative idea is that each of the Born relations (P3|-^) 
will receive a genuinely electroweak correction proportional to aVi, as given by Eq. (|180|) . 
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The main point is related to the recent observation |]HT]] that the sub-leading two-loop 
corrections to the vector boson self-energies, of order O(G^m^M^), could be numerically 
close to the leading ones, of C(G^m^). Thus the latter can be considered as an estimate of 
the uncertainties in the Vi. At the same level the 0{aa'l,ml) corrections, to be discussed 
in the next section, can also be used as an estimate of the uncertainty. In order to have 
the correct asymptotic behavior of the uncertainties for mt ^ M^, it is assumed that 
these universal corrections are multiplied by a factor 2/t = 2M^/m^. Indeed, the 0{t'^) 
is completely under control, although the sub- leading 0{t) is not and thereby 1/t follows 
with the usual safety factor of 2. Thus if the leading higher-loop corrections calculated 
in Refs. ||55| and are denoted respectively by 



6Vf , 6vfi , (104) 
then the corresponding estimates of the missing terms are assumed to be 

2 

AV = (2/t) 6V , with t = ■ (105) 

z 

There is also a specific correction to the Z bb vertex dependent on M'^/mf, which 
is presently computed up tp OiG'^^mf) term |5^, while the sub-leading terms are still 
unknown. If we denote it by 5(j/^ — see the second term in Eq. ( |22U| )) — then an upper 
bound on the 'Higgs theoretical uncertainty' can be estimated as 

Ar, = -^ro5/. (106) 

TT 

For rrit = 175 GeV and = 300 GeV this amounts to 0.02 MeV, which is much smaller 
than other uncertainties and could therefore be neglected. 



2.3 QCD Corrections on Electroweak Loops 

The effect of QCD corrections is not confined to the final-state radiation or the 0{aas) 
vertex corrections to Z ^ 66 but it will infiuence all vector boson self-energies through 
virtual gluon exchanges within quark-loop insertions. All codes include these two-loop 
diagrams [0] by first decomposing the WW, ZZ, Z7, and 77 self-energies into two basic 
building blocks, n^(mj,mj) and Il^{mi,mj), which are given by the expansion 

n„ = iv.(n<!i + ^^:^f n('i + ...) . (io7) 

with Cp = [N^ ~ 1) / C^Nc). For instance for each isodoublet. 



7-7 



i=l 
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z - z e^Y. [^l^v (p^ mi) + a- (: 



P ,'mi,mi 



i=l 



(108) 



In the limit where we neglect light quark masses, four different cases have to be considered: 
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n^{mt,mt), U^imt,mt), U^Jmt,0) , n^,^(0,0). (109) 

An important issue is related to the renormalization scale to choose for a^. The default 
that we have adopted is to select ^ = rrit for contributions from the t — b doublet while 
/i = is assumed for light quark contributions. A practical difference emerges in the 
various implementations of where sometimes the full expression is used, while in 
other cases a Taylor expanded (in q'^/mf) version has been used. A step forward has been 
made with the evaluation of the (9(aa^) corrections to Ap, the AFMT term Given 
the usual definition of Ap as 



P 



1 - Ap ' 

Ap = Ar.x,(l + 5^" + 5"^°) , (110) 
we have the QCD contribution to Ap, up to three loops and in the heavy top limit, as 



NcXtas [62 + ttsd^ j , (111) 



with tts = a^/n. As stated, this correction has been computed in the heavy top limit and 
therefore only the leading part of ^3 is available. Actually, this new calculation makes 
the QCD corrections to Ap much more stable with respect to the renormalization scale, 
since we now have 



6'^''° (p) ^ -0.910 {mt) - 1.069 (m*) + 2.609 {rrit) In (^^^ + 0(aj) , (112) 

where ag{mt) is evaluated with five flavours and the n/ that appears in the final AFMT 
result is interpreted as the total number of flavours contributing in the Feynman diagrams 
— Hf = 6. Recently the same result has been cast into a slightly different form by use 
of the notion that the corrections to Ap in terms of rrit are almost entirely contained in 
iv? {mt) / nil\ rh{mt) being the running mass evaluated at the pole mass. The correspond- 
ing differences in S'^'^^ amount to ^ 5% of the total QCD correction. Concerning the 



treatment of the AFMT term it should be noted that there is, at present, some disagree- 
ment on the value for nf — for instance, LEPTOP uses nj = 5 on the basis of decoupling 
of heavy fermions in vector theories, leading to a 2% change in 6 

A strictly related topic concerns the inclusion and the magnitude of the ti threshold 
effects on Ap and therefore on all the electroweak parameters. These effects have recently 

-"^^A recent and independent evaluation of the QCD corrections to the p-parameter has been presented 
by K.G. Chetyrkin, J.H. Kiihn and M. Steinhauser showing disagreement with the original AFMT 
result. Meanwhile, a revised version of the AFMT calculation has appeared in hep-ph 16.02.1995 showing 
agreement with the CKS result. For understanding the effect we have added a Note in proof. 
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been estimated by both dispersion relation and perturbative methods |Q however an 
uncertainty remains as to their magnitude (see Ref. for a detailed discussion). Most 
of the QCD corrections to Ap beyond the leading-order QCD terms can be discussed 
and evaluated by absorbing them into the O^a^;) term computed with an adjusted scale 



. QCD 

Ap^ = N,xt 



3 V 3 / TT 



Apr . (113) 



where Ap^^° takes into account the non- leading top effect with the scale set to /i = m^, 



as well as the light quark effects with /i = M^. For instance the result of Ref. |^ can be 
summarized by saying that it corresponds to a very high accuracy, to a scale, 

e = 0.3211°:^^° . (114) 

Incidentally, the original AFMT formulation corresponds in this language to a scale ^ = 
0.444. We note that the absorption of the non- leading QCD terms into a rescaling of 
has only been performed, so far, at the level of the leading ^2 '"'^ term, even if we have at 
our disposal the full O^aa^) correction factor. 

A comment is in order here. Our original idea was to incorporate the tt threshold 
effects through an opportune rescaling factor ^. However, the intrinsic theoretical error 
on ^ deriving from the threshold analysis is very difficult to assess, since non-relativistic 
approximations also play a certain role. In view of the present situation, the majority of 
the codes has agreed with a specific strategy: 

• The default is represented by computing the pseudo-observables to include the com- 
plete three-loop AFMT term at a scale of ^ = 1, Eq. ( p.l2| ), which is equivalent to 
the use of Eq. ( [1731) with ^ = 0.444. 



In order to estimate the size of the non-leading QCD effects, the 5^'^^ correction 



factor has been implemented according to the formulation of Ref. |5^, with a scale 



which gives the maximum variation with respect to the AFMT term — ^ = 0.248 
— and the difference between this and the AFMT calculation is used as an estimate 
of the corresponding uncertainty. It is certainly not the ideal solution, but we 
cannot rely on other analyses, not incorporating the (9(q;q;^) term. In this context a 
further subtlety arises, since the effective scale for the AFMT term (^ = 0.444) gives 
a correction outside the range required by the present Green function analysis |65 



In fact, the result of Ref. and the corresponding error estimate, Eq. ( |1 14| ) , also 
leave the AFMT expansion at the edge or even a bit outside the error range. 

2.4 Parametric Uncertainties 

The parametric uncertainties are those related to the input parameters, a(M^), 
and the masses, M^,mf, etc. Among them the largest uncertainty comes from a{M^): 
a~^{M^) = 128.87 ± 0.12 (see, however, the new results |P-|11|), while the relative un- 



certainty in is an order of magnitude smaller and that in is 50 times smaller. For 
the b quark mass we used 4.7 ± 0.3 GeV as a sort of conservative average. We have tried 
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to compare the effects of variations in tlie input parameters in a way wliicli, liopefuUy, 
will also give useful informations in the future, providing some sort of evolution of the 
uncertainties as a function of the errors in the input parameters. Assuming independence 
from the actual central values of a(M^) and of we computed the derivatives of the 
pseudo-observables with respect to a{M^),mf,. Assuming, also, that the errors are small 
enough and that the dependence is therefore approximately linear, this result will allow 
us to give the uncertainties, even when the input parameters or the errors on them change 
with time. Actually, a linear approximation is good enough for the derivative with respect 
to a^^ but not with respect to rrib- In order not to have problems with the latter near a 
local extremum we have defined a maximum derivative given by 



Vf = sign(/)- 



Ax 

6h = I max(/) - /I , 

Si = |/-min(/)|, (115) 

where / is the corresponding central value. In order to avoid lengthy tables we have 
computed this set of derivatives for our pseudo-observables at the standard reference 
point, where rrit = 175 GeV, = 300 GeV and as{M^) = 0.125. The results are shown 
in Tables and | for BHM, LEPTOP, TOPAZO and ZFITTER. 



2.5 Structure of the Comparisons 

Having introduced the main ingredients of the calculations we are now ready to explain 
in more detail the structure of the comparison. As already pointed out, we focus on 25 
pseudo-observables and vary mt,M^ and ag{M^) in a given range of values. It is worth 
mentioning that we could introduce at this point the notion of an adapted set-up, in 
the sense that prior to the introduction of options we have made an effort to show that 
different codes running under as similar as possible conditions give very close answers. 

In fact, each code will produce a set of results according to some well specified preferred 
set-up relative to the various options briefly discussed. As far as the external building 
blocks are concerned — QCD and QED correction factors — some effort has been made 
in order to reach a common default. 



The QED final state radiation for the partial widths is computed at the scale pi[] . 



The FTJR correction |5j has been split into a universal QCD factor computed at a 
scale /i = and an internal correction factor computed at the scale ii = rrit. The 
former is not included for the 6-quark asymmetries (effective mixing angle sin^ ^efr)- 



The AFMT three- loop effect |^ with n/ = 6 is included in the default with a scale 
fi = mt. 

The 0{aas) vector boson self-energies (fl^^^) are included (almost always the full 
expression) with two scales: n = rrit for the t — h doublet and /i = for the light 
quarks. 
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• The introduction of an effective scale is not part of tlie preferred set-up but 
rather is included in the uncertainty. 

This procedure for estimating the size of the non-leading QCD effects has been imple- 
mented in all codes but LEPTOP, which uses the (2/i) x AFMT term as an estimator of the 
uncertainty. 

In this way the preferred set-up of a code refers to a specific choice of the elec- 
troweak options, everything of course embedded in a given choice of the renormalization 
scheme. The result of this procedure is given by five sets of predictions for the 25 pseudo- 
observables as functions of rrit, and ag{M^). On top of our predictions, each group 
has adapted the various codes to run under all its options (typically up to 2^ 2^). It 
must be clearly realized that no common set of options has been created and that the 
options of one code have been designed independent of the options of all other codes. 
In the end, the 2" electroweak options are folded with 2 QCD options, the inclusion of 
AFMT being the default versus a rescaling 0.248 mj in the 0{aa^) term represent- 

ing the uncertainty. LEPTOP, however, is using a different procedure (see above). Note, 
that electroweak-QCD factorization simulates 0{aag) non-controlled terms, therefore it 
can also be considered a QCD option. For each pseudo-observable O we have collected 
Oadapt and 0+,0-, i.e. 

0+ — max Oi 

i 

C>_ = minOi (116) 

i 

where the index i is running over the options. The differences 0+ — Oadapt and Oadapt ~ 
0_ calculated by a given code are our internal estimates of the theoretical uncertainty 
associated with O. while the different results for Oadapt as obtained by the various codes 
may be considered as giving our estimation of the scheme dependence. Here internal must 
be understood as the estimate that one particular code can produce by varying internally 
its options on the implementation of radiative corrections. The corresponding error bars 
on the theoretical predictions are in some cases very asymmetric, merely reffecting the 
specific ideas or, even better, personal taste, beyond the choice of preferred set-up. Clearly 
the way in which the different realizations have been built into the codes is very indicative 
of the original strategy. To give an example, we may note that the reason why some code 
does not include the one-loop corrected axial-vector coupling of the Z to fermions in the 
definition of the p parameter is because its preferred set-up is the expanded solution. In 
some cases, the theoretical uncertainty internally estimated by one code could turn out to 
be large with respect to those of other codes. Basically, we do not attribute any particular 
relevance to this fact, as the global indication should include some sort of average among 
the codes. In the general discussion of the results we have tried as much as possible to 
trace the roots of the phenomenon, whenever it appears, and any further consideration 
should be left to the potential users of our analysis. Each point in the error bars has the 
same content of probability and the width of the theoretical bands associated with each 
pseudo-observable should be taken as an indication of the relevance of that quantity for 
the analysis of the LEP data. Sometimes the theoretical error, which cannot be reduced 
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unless we come up with a full two-loop calculation, becomes of the same order as the 
present or projected experimental error. Even before entering into a full discussion of the 
results we may anticipate the predictions for sin^^^g by choosing some reference point, 
such as rrit = 175 GeV, = 300 GeV and ag{M^) = 0.125. The prediction is 



eflf 
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The corresponding width for the theoretical band is therefore 0.00011, 0.00016, 0.00008, 
0.00010 or 0.00018 for BHM, LEPTOP, TOPAZO, WOH and ZFITTER. This can be compared 
with what we obtain by combining the LEP results on all the asymmetries: sin^ = 
0.2321 ±0.0004 



In order to present our results we have focused on the W mass, on the primary set 
of pseudo-observables chosen by the LEP collaborations for fitting — r^,Ri,A'- , and on 
the b- and c-quark related charge asymmetries and ratios of partial widths. 

2.6 Experimental data and theoretical predictions 

To set the scene for our discussion we must first introduce, in Table ^, the relevant 
experimental data. Clearly a detailed discussion of our results should take into account 
the whole range of values for the input parameters but many of the relevant conclusions 
can already be drawn by considering Tables where we have reported 



TTI and a 



sample of 11 quantities as measured by the LEP collaborations. For comparison we have 
fixed a reference point, rrit = 175 GeV, Af^ = 300 GeV and = 0.125, and reported the 
predictions from BHM/LEPTOP/TOPAZO/WOH/ZFITTER, including the estimated theoretical 
errors and the average. The content of the tables should not be confused with a fit but 
should only be taken as a first introduction to the theoretical results. Our reference point 
is a mere consequence of the range indicated by CDF for the top quark mass and of 
the most recent prediction for a^(M^) from LEP. As far as is concerned, we must 
admit a high degree of arbitrariness. As we have detailzed, the differences among the 
(theoretical) central values for each quantity are basically (even though not totally) a 
measure of the effect induced by a variation in the renormalization scheme. As can be 
seen from Tables ^-|^ the ratios of the maximal half-differences among the five codes over 
the experimental errors are: 



Ae(M^) = 2.5 X 10 



-2 
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Ae(re) 


= 6.7x10' 


-2 




= 3.9 X 10" 


-2 
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-1 
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FB ' 


= 5.6 X 10" 


-2 


Ac(A^ ) 


= 7.8 X 10" 


-2 




= 2.5 X 10" 


-2 



Another piece of information is obtained by looking at the theoretical uncertainties as 
estimated internally by each code. A very conservative attitude would be to report the 
half-difference between the global maxima and minima among the five codes. By consid- 
ering again the ratio with the experimental errors we obtain 



A,(M,) 


= 7.2 X 10" 


-2 
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= 1.6 X 10" 


-1 
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By comparing the various Ac and we can obtain a rough evaluation of the global theo- 
retical error associated with the most relevant quantities in the analysis of the LEP data. 
As far as the differences among codes in their preferred(adapted) set-up is concerned, 
we can safely conclude that the ratios of their predictions to the experimental errors are 
usually less than 0.15. However, the most important message to be derived from this 
simple exercise is that very often the theoretical uncertainty can be larger than what is to 
be expected on the basis of a simple comparison of the results from different calculations. 

A wider sample of results is shown in Figures 11-23, where, again, we have fixed some 
reference point — = 300 GeV and a^{M^) = 0.125 and where 100 GeV < rrit < 
250 GeV. In every figure the corresponding experimental points(data), as they are given 
in ref. (|ll5|]), are shown at = 178 GeV. The full collection of our results refers instead 
to = 60, 300, 1000 GeV and to ^^(M^) = 0.118, 0.125, 0.132[3. 

Here we discuss the main features of the comparison and for each quantity we indulge in 
presenting an estimate of the global uncertainty by roughly considering the half-difference 

^^A preliminary version of the comparison among our results has been presented in Bq] 
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between the maximum and minimum among all predictions. Admittedly this is not a 
very rigorous procedure and therefore it should be treated with due caution. It should 
be stressed again that the Ac, factors presented above show the half-difference of 
the predictions over the experimental errors. In the following discussion we will mainly 
analyze two quantities: 

dci dg — the half-differences, either among central predictions or between the maximum 
and the minimum among all predictions. 

It should be noted that we will not address in this section the question of how the transfor- 
mation from primordial distributions to the secondary quantities will affect their precision. 
In the following we have examined some of the pseudo-observables at the standard refer- 
ence point and tried to present the state of the art for their theoretical predictions without 
including all sorts of parametric uncertainties, but rather we have limited our discussion 
to the genuinely theoretical ones. Below we discuss 13 pseudo-observables. 

• 

There is a certain spread in the predictions for increasing values of rrit substantially 
independent of and increasing for large M^, with the formation of two clus- 
ters, represented by BHM/WOH and ZFITTER on one side and LEPTOP/TOPAZO on the 
other. The maximum difference in central values is seen for rrit = 250 GeV where 
dc reaches a 10.5(12) MeV for = 300(1000) GeV, made even more significant by 
the addional fact that the bands essentially do not overlap. The situation improves 
for central values of m^, where we see at most a half-difference of ~ 4.5 MeV even if 
the clustering is already evident. 

No pattern of any particular relevance is observed. Almost independent of we 
notice that for intermediate rrit there is a substantial agreement for all M^, while for 
high rut the agreement is better at large values of , and for small rut it improves 
for low values of M^. All the error bars tend to become wider for increasing rrit 
and M^, with the possible exception of TOPAZO/ZFITTER. The largest half-difference 
among central values is for low and high rut where dc can reach ^ 0.03 MeV. A 
safe estimate of the overall theoretical error at rrit = 175 GeV is of about 0.03 MeV 
(0.030 <dg< 0.031 for the full range of and M^). 

• 

Due to the final state QCD corrections is much more sensible to variations in a^. 
However, we have verified that there is no substantial variation among the codes as 
a function of a.,, a sign that final state QCD corrections are well under control. For 
instance we find 1.4 MeV < dg < 1.5 MeV for 0.118 < < 0.132 at rut = 175 GeV 
and = 300 GeV 0- For low rrit the result of the comparison does not show 
any particular pattern, while for high rrit the agreement improves with high M^, 
showing again some sort of correlation between the two variables. For intermediate 

""^^This is not a trivial consequence of the fact that all five codes use the same radiation functions, since 
their implementation is usually different. 
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rrit, instead we find smaller differences around central values of M^. In general for 
rrit > 150 -7- 175 GeV the various error bands, while growing, have a tendency to 
overlap. The maximum deviation among codes is for high rrit and low M^, where dc 
can reach approximately 0.6 MeV. For rrit = 175 GeV = 0.125 instead we get, as 
the global estimate of the uncertainty, 1.3, 1.4 and 1.4 MeV for = 60,300 and 
1000 GeV (with tiny variations in respect of a^). 

Ri 

This quantity has a role of its own since quite often it is used for extracting a precise 
determination of a^^M^). Indeed, up to some degree of accuracy, the two variables 
are related by Aa., ~ nAR/R so that a difference of 0.01 in Ri is equivalent to an 
error of 0.002 in a^. There is some common trend in all our results for Ri, namely 
the BHM predictions always stay a little higher, while the other codes tend to cluster, 
apart from the TOPAZO results, which, for very low rrit and high tend to converge 
towards those of BHM. For rut > 175 GeV the error bands tend to overlap so that 
each code has a central prediction within the other error bands, again apart from 
the BHM point, which sometimes is fully contained only within the WOH predictions 
and lies at the upper boundary of the LEPTOP/ZFITTER ones. Error bars are often 
very asymmetric, expecially for WOH/ZFITTER. The maximal deviation is for high 
mt,Mj^, where dc may reach 0.006, whereas the global estimate of the uncertainty 
for mt = 175 GeV gives 0.0085,0.0090 and 0.0095 for a^{M^) = 0.118,0.125 and 
0.132, independent of M^. These values correspond to an error of 0.002-^-0.003 in the 
determination of a^^M^). We have also analyzed in more detail the dependence 
of the ratio Ri for rrit = 175 GeV and = 300 GeV, including the case when 
QCD is switched off. Indeed, a determination of from Ri is usually achieved by 
writing Ri = Ri^a^ = 0)(1 + (5*^*^°) and by using the most updated formulation of 
QCD corrections (see for instance Ref. |6^). In this way, some relevance should 



also be attributed to a comparison of various predictions for the ratio Ri, unfolded 
of QCD correction terms. For = the WOH prediction is lower than the others, 
which cluster around 19.946 (WOH is —0.78 x 10"^ below the average). When QCD 
corrections become active, BHM remains higher (0.5 -h 0.6 x 10~^ above average) while 
the other codes form a cluster. 

Rc 

The ratios of the bb and cc partial widths to the total hadronic width share some 
common features. The experimental errors are 0.0020 and 0.0098 and the two quan- 
tities are —38% correlated. Our central predictions for Rb are all within 2 x 10~^ 
and even the inclusion of the theoretical uncertainty only gives dg ^ 6 x 10~^. 
The overall theoretical error at rrit = 175 GeV is 1.7 X 10"^. For Rc the global 
uncertainty is well contained within 6.0 x 10~^. Rb shows some clustering which 
becomes more and more evident for large rut, with BHM giving the higher prediction 
and LEPTOP/TOPAZO/WOH/ZFITTER forming a lower cluster. The behavior of dg as a 
function of is practically flat. 



sm 



2 /ilept 



The reported value of the leptonic effective weak mixing angle is the average of all 
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forward-backward and polarization asymmetries from LEP. Therefore the analysis 
relies on the hypothesis that the peak forward-backward asymmetry can simply be 
connected with the remaining asymmetries through the use of the same sin^ O^q. For 
all values of and the agreement is less satisfactory for low rrit where ZFITTER 
remains on the higher side, BHM/WOH form a lower cluster and LEPTOP/TOPAZO are 
somewhere in between. The general trend is to have a convergence of all codes 
for high mt- The maximal half-difference among central values is as for low rrit, 
where for all M^,ag we find dc ~ 9.0 x 10~^. For — 175 GeV the overall 
uncertainty is estimated to be ci^ ~ 1.1 -i- 1.4 x 10^*^ over the whole range — a^. 
For = 300 Gev and = 0.125 we find dg = 1.3,1.1,2.1 x 10"^ for mt = 
100, 175, 250 GeV. In view of the supposed relevance of this parameter and of its 
projected experimental error we have to admit that the status of the theoretical 
predictions is far from satisfactory but totally related to the unknown higher order 
effects. To give an example, we could say that the knowledge of the sub-leading 
0{G'j^M'^mf) corrections to Ar would greatly improve the situation — for instance 
for ZFITTER, which dominates the error. 

• sin'' els 

From many points of view the situation is very similar to that described for sin^ ^'g. 
Let us remember that sin^ O^q differs from sin^ O^q because of flavour-dependent cor- 
rections, which are mj-dependent and not negligible. There is some kind of crossed 
behavior in our predictions, with an agreement substantially better for intermedi- 
ate values of mt and deteriorating at the boundaries. For low mt the comparison 
is similar to that for sin^ 91^: a higher prediction from ZFITTER, a central clus- 
ter LEPTOP/TOPAZO and a lower one BHM/WOH. For high mt the highest prediction 
is from BHM, with a lower cluster of the other codes. The global estimate of the 
uncertainty for mt = 175 Gev and = 0.125 is of 1.15(1.05,1.45) x 10"^ for 
— 60(300, 1000) GeV, with an uncertainty ±0.1 x 10~^ due to a variation in a^. 

FB 

In presenting results for the leptonic forward-backward (peak) asymmetry, as well 
as for any other leptonic asymmetries, we follow the indication of the experimen- 
talists who keep A^^^ until the end in their standard model fits, without necessarily 
identifying it directly with sin^6'gg. We first start by analyzing the comparison 
among the central values. Here the result of the comparison is rather good and 
essentially we always start at low mt with two clusters — a higher one containing 
BHM/WOH and a lower one containing TOPAZO/ZFITTER and with LEPTOP somewhere 
in between. There is a fast convergence of all results for increasing mt, expecially for 
high M^. Typically for mt = 175 GeV and = 300 GcV we find a half-difference 
of about 8.4(8.9,9.5) x 10"^ for = 0.118(0.125,0.132). When we come to the 
inclusion of the theoretical uncertainty it is immediately evident that the previous 
comment still applies and the global error for the standard reference point becomes 
1.9 ^ 2.0 X 10^^ at a, = 0.118, 2.0 ^ 2.2 x 10"^ for a, = 0.125 and 2.0 ^ 2.3 x 10~^ 
for ctg = 0.132 , where the variation with is illustrated. We will come back to 
A''^^ and to the related uncertainty while discussing additional theoretical options 
which have not been included in the working set. 
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Also for the left-right (peak) asymmetry there is some general behavior in our pre- 
dictions. At low rrit we start with maximal disagreement, two clusters with BHM/WOH 
in the higher one and LEPTOP/TOPAZO in the central one, to reach convergence of 
results for high m^. For low rrit there is also a considerable spreading of the error 
bands. Given the fact that the SLD measurement seems to require a much higher 
value of rrit, we have considered the overall uncertainty for = 250 GeV with the 
result of 1.82(1.59, 1.74) x lO'^ for M„ = 60(300, 1000) GeV for = 0.118. These 
values become 1.84(1.61, 1.74) x 10"^ and 1.86(1.63, 1.77) x 10"^ for = 0.125 and 
0.132. 

FB ' FB 

Practically everything we state for the general behavior of A^^^ can be repeated for 
the b- and c- quark charge asymmetries. As for A^^, the agreement is worst at 
low rrit and there is a general convergence at high when rrit is growing. For 
low and large rrit TOPAZO has the tendency to give a higher prediction. The 
global uncertainty for A^ at rrit = 175 GeV is 6.2 -^ 7.4 x 10~^. There is no general 
agreement among codes on the definition of ^^g. In particular LEPTOP and TOPAZO 
include mass effects into the definition of this pseudo-observable. For the b quark, 
the effect of its non-vanishing mass leads to: 

where v is the b quark velocity: 



4m2(M2) ^ ^ 



where rhf,{M'^) is the running 6-mass in M5'-sheme defined by Eg . (|59|) below. How- 
ever this is a very tiny effect. An estimate of it, based on TOPAZO results, gives an ab- 
solute deviation between the massless and massive definitions of 1.6, 1.0, 2.0x 10~^ at 
rrit = 100, 175, 250 GeV for = 300 GeV and = 0.125. Similar conclusion can 
be drawn for A'^^. The global uncertainty at rrit = 175 GeV is dg = 4.5 -i- 6.0 x 10~^. 

For the hadronic cross-section the theoretical error dg is always contained in a range 
of 0.0075 -r- 0.0095 nb, with very little dependence on a^,. Indeed, for our standard 
reference point {mt = 175 GeV and = 300 GeV) the largest half-difference 
among central values is 0.0065 nb, 0.0065 nb and 0.0070 nb for = 0.118,0.125 
and 0.132. Our estimate of the overall theoretical error at rut = 175 GeV is of 
about 0.009 nb and the largest half-difference is seen at large rrit and where it 
reaches 0.0095 nb. Over the whole range of the parameters the prediction of WOH 
remains higher with a tendency for ZFITTER, though only at very low values of rut, 
to slightly converge towards WOH. In any case, the WOH error bands are sufficiently 
large to include the other predictions or to overlap with the other bands. 
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To summarize our results we have also presented in Table § the largest half-differences 
between central values or between maximal and minimal predictions among codes in the 
range 150 GeV < < 200 GeV, 60 GeV < M„ < 1 TeV and 0.118 < a, < 0.125. In 
Tables ^ abd ^ we have fixed mt = 175 Gev and illustrated the largest half-differences 
among central values or those between the maximum of the maxima and the minimum 
of the minima among the five codes in two situations: = 0.125 fixed and 60 GeV < 
< ITeV, or = 300 GeV and 0.118 < < 0.125. In this way the separate 
contributions from or to the theoretical errors are indicatively given. A final but 
partial indication of our results can be provided by computing some of the quantities 
which usually enter the SM fits (at = 175 GeV) and by collecting all available sources 
of error 



BHM = 2497.4+°:|^(t/i.)+^;^(M^,aJ ±6.7Aa-^ ±0.8 Am^MeV 

LEPTOP = 2497.2 ± l.l(t/i.)+|°(M^,a5) ±6.8 Aa^i ±0.8 AmfeMeV 

TOPAZO = 2497.4l°j(t/i.)+^;2(M^,aJ ±6.8A«"^ ±0.8AmbMeV 

WOH = 2497 Atli{th.)tli{M„,a,)MeV 

ZFITTER = 2497.5l|]|(t/i.)+^:?(M^,aJ ±6.8A«-^ ±0.8AmbMeV (121) 



BHM = 41.436^|^:|][{^(t/i.) ±0.042(M^,aJ ±0.013 Aa~^ ±0.007 Am^nb 

LEPTOP = 41.439+°:°°^(t/i.) ±0.040(M^,aJ ±0.012 A«-^ ±0.008 Am^nb 

TOPAZO = 41.437l[5:[![|^(t/i.)lJ]:[^^J(M^,aJ ±0.012 Aa"^ ±0.007 Am^nb 

WOH cTo^ = 41.449l°:°?°(t/..)l°:°^J(M,,a,)nb 

ZFITTER = 41.44i;°:°°°(t/i.)+°:°^J(M^, ± 0.013 Aa"^ ± 0.007 Am,, nb (122) 
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where we have allowed, as usual, 60 GeV < < 1 TeV and 0.118 < a^, < 0.132. 
Whenever available, the parametric uncertainties have been inserted. It is assumed that 
Arrtf, is given in GeV. 

2.7 More on Theoretical Uncertainties 

There are some options on electroweak radiative corrections which, although implemented 
in the codes, have not been used as working options for producing our comparisons with 
the experimental data. The main argument for this exclusion is related to their tendency 
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to produce results deviating sensibly from the average. As a matter of fact, this is the 
visible consequence of their inclusion, but quite often we have theoretical reasons against 
them. This short section is devoted to a summary of those effects and of their eventual 
influence on the comparisons. At the end one should not forget that the design and 
the implementation of options is indeed something very peculiar to a given realization. 
The following short considerations give another illustration of a very important fact: 
theoretical uncertainties should be treated with due caution, realizing that they contain, 
in any case, a large degree of arbitrariness. 

2.7.1 TOPAZO 

With TOPAZO, the most striking effect is related to an expansion of the non-leading terms 
to 0{a), once bosonic sclf-cncrgics were not re-summed. Roughly speaking, we can assume 
that a certain quantity X is given by (see also 2.2.4) 



where Xq, by construction, will include all the re-summed contributions. As soon as we 
allow the two options. 



and interpret the resulting difference as a theoretical uncertainty, then, under some cir- 
cumstances, we end up with errors considerably larger than those presented above. The 
main problem is represented here by a clash between accidental cancellations and gauge 
invariance (even the notion of cancellation between fermionic and bosonic sector is gauge- 
dependent). As it is well known, the bosonic self-energies are gauge- dependent and, more- 
over, the fermionic ones tend to dominate away from the intermediate rrit region. That is 
how TOPAZO allows, among other options, for a strict resummation of the fermionic self- 
energies alone. Of course, one could just avoid resummation altogether, but the option of 
expanding versus non-expanding the remainders, with respect to the leading terms, still 
applies and the theoretical uncertainties would become, in this case, sensibly nit depen- 
dent. Finally, we stress once more that the identification of a gauge invariant part of the 
bosonic self-energies is not at all a unique procedure and therefore some degree of arbi- 
trariness is always hidden in a global resummation. We admit that, this option allows for a 
nice construction of a small remainder, in a situation, however, where we are working with 
one-loop contributions and higher order reducible ones. What is left out of our analysis 
is, in any case, related to the two-loop irreducible terms, about which nothing is presently 
known. To summarize we could say that in one case it is approximately the 'square' of 
the one-loop bosonic corrections that we use to estimate the theoretical error, while in the 
other we can decide to re-sum part of it into the leading term and make the remainder 
small, even though we are still missing information about not-yet-computed higher orders 
and their effects, which could make the smallness of the remainder inadequate. 



X = Xo + -X, + 0{a^) , 



(124) 
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The final reason why TOPAZO has excluded this option in presenting the pseudo- 
observable Tables is therefore totally related to the abnormal (as compared to the other 
codes) size of the errors for some of the quantities, noticeably A' . Here we say again 
that the expansion option is bounded to produce larger errors whenever the remainders 
are not (one way or another) kept small and in some codes this effect is not seen, sim- 
ply because nothing equivalent to the expansion has been implemented. To give a quick 



idea of the effect we present in Table |Tl| the shift in the central values and in the error 
bands for some of the pseudo-observables. Clearly the largest effect is seen for the lep- 
tonic forward-backward asymmetry where the theoretical error becomes comparable in 
size to the experimental counterpart. Needless to say, when this option is activated, a 
large fraction of the uncertainties become TOPAZO dominated. Additional consequences 
will be introduced and discussed in the next chapter. 

2.7.2 ZFITTER 

About ZFITTER, one can also mention several peculiar moments, related to the specific 
design of its options. We begin with by making clear that nothing resembling what is 
described in the previous section was observed. To a large extent this is due to the fact 
that the coefficient Xi (see discussion in 2.2.4) happened to be small in the framework 
of the ZFITTER renormalization scheme. Having accepted this statement, however, one 
should not conclude that it has particular advantages over the other schemes. The size of 
the coefficient is simply a numerical accident, without any deep physical meaning. 

However, ZFITTER does contain additional options, which eventually were not left 
among its working options. An example is given by the array of expansions (pOf ). All four 
expansions have been implemented, and it was noted that the third and fourth expansions 
enlarge the theoretical uncertainty for some observables (M^ and sin^Olg), roughly by 
factors of two and four correspondingly. As it was pointed out in section 1.10.2, the third 



and the fourth expansions contradict to the conclusion of paper |^ about fermionic mass 
singularities. It not surprising that this was the argument in favour of excluding these 
options from the working set. 

Another interesting example is related to the leading-remainder splitting problem. 
Three variants of the resummation of the leading terms in Ap, see Eq. (|80|), were imple- 
mented as different sub-options: 

1) only the first term, i.e. Ap", is re-summed; 

2) all but X terms of (|80|) — the content of square brackets is re-summed; 

3) the whole expression (ROl) is re-summed (ZFITTER default). 



There is a noticeable increase in the uncertainty when we include the second option with 
respect to a situation where only the first and the third are retained among the working 
options. One should emphasize, however, that this increase in the error is not dramatic: 
for example for sin^ it amounts to a nearly rrit independent uncertainty, of the order 
5 X 10~^. Examining the reasons for this enhancement, it was revealed that for the second 
option the remainder terms are about 5-10 times bigger compared to those for the first and 
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the third options. Moreover, they are only several times smaller than the leading terms. 
On the basis of this observation, the second option was termed a pathological option 
and initially excluded from the working set, since its effect contradicted to an accepted 
strategy. In playing more with ZFITTER options a striking property was observed. Some 
of the ZFITTER options do not possess the additivity property. As an example, we mention 
that the 'scale of remainder option' alone produces nearly the same uncertainty as when 
it is applied 'in conjunction' with the following two other options: 

- the one dealing with expansion (pOf), first two rows; 

- the one dealing with the three above-mentioned variants of resummation for Ap: 
l)-3). 

As a consequence, the combined effect of all three options leads practically to the same 
uncertainty, and the latter is independent from the actual number of options included — 
two or three Ap resummation sub-options are left in the working set. This eventually 
lead to a decision to retain all 1) -r- 3) among working options. 
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3 Realistic distributions 



To give predictions for pseudo-observables has not been our only task and we have also 
devoted a noticeable effort in order to present the most updated analysis for realistic 
observables. This process requires as fundamental ingredients a QED dresser and therefore 
the comparisons have been restricted to BHM, TOPAZO and ZFITTER, since they allow for 
the treatment of QED diagrams involving the emission of real photons with results which 
are dependent on energies and experimental cuts. Roughly speaking we can distinguish 
between s-channel processes and Bhabha scattering, e^e~ e^e~ . For Bhabha scattering 
the commonly accepted procedure is the so-called t-channel subtraction, where the s — t 



and t — t contributions are subtracted from the data by using the code ALIBABA |^8[. In 
the procedure the 'most reasonable' values of rrit and of are used leading to additional 
sources of errors in the analysis. Given the present situation we have also performed a 
comparison between ALIBABA and TOPAZO inspite of the fact that the electroweak library 
of ALIBABA has not been constantly updated. Thus in the comparison will emerge an 
intrinsic difference due to the improved electroweak and QCD formulation of TOPAZO. 
Our comparisons can be divided according to the following scheme: 

• Fully extrapolated set-up for muonic and hadronic channel with the possible inclu- 
sion of a cut on the invariant mass of the final fermion pair, the so-called s' cut 
(BHM, TOPAZO, ZFITTER); 

• e+e" ^ (e+e-, s-channel for 40° < 6. < 140°, = 10°, 25°, and ^th(/i^) = 
20GeV (^th(e^) = IGeV), (TOPAZO, ZFITTER); 

• e+e- ^ e+e- for 40° < 9^ < 140°, = 10°, 25°, and ^th(e^) = 1 GeV, 
(ALIBABA, TOPAZO). 

Moreover, we have fixed the following set of values for the c.m.energy: 88.45, 89.45, 
90.20,91.1887,91.30, 91.95, 93.00, 93.70 GeV. All the resuhs refer to a given reference 
point, rrit = 175 GeV, = 300 GeV and ag(M^) = 0.125. A first comment con- 
cerns the s' cut, which should not be confused with a cut on the invariant mass of the 
event after initial state radiation (sometimes also used in the experiments). For realistic 
observables we have avoided all references to any specific set of effective formulas and to 



their realizations, the interested reader having available the existing literature ||68|1-||69 
Weak radiative corrections, depending on the assumptions of the electroweak theory, have 
already been discussed from the point of view of the options that arise in their practical 
implementation. Here we only mention that we have fully propagated those theoretical 
uncertainties from the pseudo-observables to the realistic ones with the result that — for 
the fully convoluted cross-sections and forward-backward asymmetries — the final results 
include a theoretical error bar. In the following we present a short discussion of the main 
ingredients entering the calculation of realistic observables and critically compare some 
of the results obtained with BHM, TOPAZO and ZFITTER. 
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3.1 De-Convoluted Distributions 

To illustrate in more detail the construction of realistic observables (RO) we start from 
the concept of de-convoluted quantities. For a given process we construct (T^(cTg), the 
forward (backward) kernel cross- sect ion, including electroweak corrections, and eventually 
the comprehensive of a cut on the angular acceptance. QCD corrections are included 
while all QED corrections are left out. After a first comparison at this level we proceeded 
by introducing QED final-state radiation (FSR), initial-state leptonic and hadronic pair 
production (PP), initial-final-state QED interference (INT) and, finally, the kernel distri- 
butions folded with initial-state QED radiation (ISR). One should emphasize, however, 
that the INT contribution was simply added — at 0{a) — and was not folded with the 
ISR. As far as the ISR is concerned we did not find it opportune to fully review the various 
treatments and implementations but have tried as much as possible to illustrate the origin 
of possible discrepancies whenever they arise. For instance, the differences that we find 
in the results of the various codes are dominated by pure weak (and QCD) corrections 
around the region of the peak and by QED radiation along the tails, where, however, 
the experimental error is considerably larger. Among the de-convoluted quantities, the 
most relevant are those computed at s = M^, which have an obvious counterpart in the 
pseudo-observables, that we have already computed, A'' and a'^, which hereafter will be 
characterized by an index 0, as A''*^ and ctq. There is, however, a noticeable difference 
between the two sets, represented by the interference of the Z — 7 s-channel diagrams 
and by the presence of imaginary parts in the formfactors, the latter being particularly 
relevant for the leptonic forward-backward asymmetry. 

As far as the propagation of electroweak uncertainties from pseudo-observables to 
realistic ones is concerned we notice that all three codes see an enhancement of the 
theoretical errors. For instance for the standard reference point we denote with ctq the 
PO hadronic (peak) cross-section and with (T^(M^) the realistic one and get: 

BHM cr^ 
TOPAZO CT^' 
ZFITTER 

the enhancement factor being 1.9, 2.3 and 1.4 for the three codes respectively. The reason 
behind this increase of the induced theoretical error is that for ctq, as defined in Eq. (|3lD , 
we first compute Te,Th and and then construct the combination FeF/i/F^ without 
any further expansion in a, while for (j^{M'^) we include the expansion in the option set 
(linearization of the cross-section), thus enlarging the error. 

3.2 Final-State Radiation 

A substantial difference exists between fully extrapolated RO and RO in the presence of 
cuts. This is illustrated bt the fact that without cuts we can simply use the well known 
correction factor 1 + for each partial channel and there is therefore no ambiguity in 

FSR. A possible source of discrepancy can instead be introduced when cuts are present, 
due to a different treatment of final-state higher-order QED effects. This can lead to 



= 41.436l°:°03nb a\Ml) = 30.366^°;°^^ nb 
= UA37iro>^ <y\Ml) = 30.375^°:°^^ nb 
= 41.44lino5nb a^Ml) = 30.373t°:°°^ nb , (125) 
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differences wliich in general depend on tlie experimental cuts required and may grow for 
particularly severe cuts. It has already been shown that two possible prescriptions, — 
completely factorized final-state QED correction versus factorized leading-terms and non- 
leading contributions summed up — lead to differences in the cross-section for Bhabha 
scattering of the order of 0.5% far from the peak, whereas the asymmetry is substantially 
left unchanged. To be more specific, the final-state QED corrections amount to a leading 
correction term: 



cut V 



2-Q)\n 



So 

s 



In 



1 



(126) 



where Sq represents a cut in the reduced invariant mass. Renormalization group argu- 
ments suggest that such a leading term should be exponentiated, which is of no practical 
importance at low thresholds but could give sizeable effects at high thresholds. By defining 



-^CUt,!!-^) 



the leading term resummation can be implemented as follows: 



exp 



l + F^^,(s)-F^,,(s)Fl,{s 



(127) 



:i28) 



where spurious terms Fl^^{s) x -F(.ut,r('S) are confined at least at 0{a^). Indeed, several 
prescriptions for treating final-state correction are possible, all equivalent at 0{a). One 
reasonable recipe could be to define the leading term in a different way. For instance, in 
the presence of an acoUinearity cut, the infrared logarithm could be defined as 



I = ln(l - x) 



(129) 



where x is given by 



X = max(so/s, yr) , 



Ut being 



1 - sin(C/2) 



(130) 



(131) 



1 + sin(C/2) ' 

and ( the maximum acoUinearity allowed. Another possibility would be to exponentiate 
the full 0{a) contribution F^i ^^ even though there is no guarantee that the experimental 
cut-dependent terms do exponentiate; in this case spurious terms appear already at O(a^). 
Alternatively, one could choose to factorize only a leading 0{a) term and simply add the 
0{a) correction due to the acoUinearity cut (this simulates the choice in Ref. 



3.3 Initial-State Pair Production and QED Interference 

Next we come to the inclusion of initial-state pair production in the realistic distribu- 
tions. A fermionic pair of four-momentum radiated from the e"*" or e~ line gives a 



correction |]70| : 
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^rmpair ^pair ' ^pair 



Urn? 7(v^-A)2 dq^ds' 

^ /. . 33X7 • (132) 



dz Jim? dq^ds' 

Also for this term there are different treatments — we can exponentiate the pair produc- 
tion according to the YFS formahsm |^ or the same pairs can be included at 0{a^). In 



the end, however, we found a reasonable agreement among the results of the three codes 
relative to the specific inclusion of pair production. The main features of this correction 
term are as follows. 



TOPAZO/ZFITTER employ the KKKS formulation |T0[ without exponentiating the 
soft+virtual part, which is added linearly to the cross-section. BHM, instead em- 
ploys the YFS formalism [|71|. The independence of the results from the soft-hard 



separator has also been successfully investigated. 

• r-pairs are not included. 

• The lower limit of integration Zmin, adopted for leptonic pair production, is 0.25 and 
the soft-hard separator A has been fixed in the region where we see a plateau of 
stability. 

Initial-final QED interference has been introduced in the calculations, including the 
effects of hard photons. For a fully extrapolated set-up this means that all photons, 
up to the maximum available energy, are taken into account, while s' cuts or energy 
thresholds and acoUinearity cuts will restrict the available phase space. In order to proceed 
step-by-step we have introduced the procedure of comparing our results in the sequence 
ap^g^{NN,NY,YN,YY), where the first argument in parenthesis denotes inclusion or 
exclusion of PP, and the second refers to the interference. In this way the relative influence 
of QED corrections has been checked and kept under control. To illustrate the trend, we 



present in Tables |T2| and |13| a comparison for the hadronic cross-section and the muonic 
forward-backward asymmetry. In particular, to continuously keep under control our 
comparisons, we have introduced and analyzed the ratios 

, c/O (C-TOPAZO, conf) 

r C-TOPAZO, conf, O) = — ^— '——^ , 133 

^ ' ' 7 c/0 (C-TOPAZO, NN) ^ ^ 

where O = a^, cr^, A^^, conf = YN/NY/YY and C = BHM.ZFITTER. Moreover, dO denotes 
the relative variation for cross-sections and the absolute deviation for the asymmetry. 
Whenever pair production or QED interference have similar effects, the corresponding 
ratio assumes values of around 1. When the ratio goes to zero it is a signal that the 
apparent agreement does not reflect a similar agreement in the NN quantities and is 
therefore due to accidental compensations. Finally, if this ratio grows in modulus it gives 
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an indication that the agreement at the NN level is not respected when PP or INT are 
introduced. 

To illustrate this fact, consider Table |T^, which refers to . We find, for instance, 
that the YY cross-sections for ZFITTER-TOPAZO agree in five digits at = 89.45 GeV, 
90.20 GeV and 91.30 GeV, giving 10.042 nb, 17.992 nb and 30.514nb respectively. More- 
over, a closer look reveals that 



= 89.45 GeV 10.067 - 10.068 nb 
v/s = 90.20 GeV 18.039 - 18.040 nb 
= 91.30 GeV 30.590 - 30.590 nb 

(134) 



for ZFITTER/TOPAZO in the NN configuration. 



3.4 Imaginary Parts of the Formf actors 

The presence of imaginary parts in the weak formfactors introduces additional possibilities 
for the implementation of radiative corrections with respect to pseudo-observables. At 
the level of the kernel cross-sections and remembering the quite general subdivision into 
leading parts and remainders introduced in section 2.2.1, we can basically select two 
possible options: 

• Imaginary parts confined in the remainders, 

• Imaginary parts inserted into the leading terms. 

To illustrate the effect of imaginary parts we select / = yU and consider the asymmetric 
part of the angular distribution for e"'"e~ A*'*'/^"- If we denote with x-) ^"^^ Xz 
corrected 7 — 7 and Z — Z propagators, this asymmetric part is proportional to 



oc TZe 



{9*^)^X7X1 + [^e (^^,^*)]^ Ix^l^ + boxes 

= T^egllmxjlmxz +2neg^lmg^nexjlmxz 

+ {TZegj' {Ueg^f + 0(7^exJ + O [Ueg^ x Xm^) + O (im^) ,(135) 

where TZeXz is suppressed around the Z peak and TZeg^ is usually small. 



3.5 Initial-State QED Uncertainties 

An estimate of the theoretical uncertainty due to QED radiation can also be derived, 
so that our results may contain two sources of theoretical error, ±A{EW) ± A{QED). 
For instance such an estimate can be performed by using the following algorithm. Since 
the main source of QED theoretical uncertainty is the treatment of final-state radiation, 
namely the use of a completely factorized formula versus a leading-log factorized plus a 
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non-log additive one, whenever the theoretical error is required, we can run over the two 
possible options and return the corresponding uncertainty. Moreover, when the large- 
angle Bhabha scattering is considered, we must realize that the main approximation 
adopted by TOPAZO is to treat the t and s + t contributions to the cross-section at the 
leading logarithmic level. Thus the size of the convoluted t and s + t terms is computed 
as the difference between the full Bhabha prediction and the pure s-channcl one and the 
theoretical error is estimated by assuming 1% of the difference. In particular, this means 
that the QED theoretical error depends on the detailed experimental set-up, growing when 
tightening the experimental cuts and, for the Bhabha case, when enlarging the angular 
acceptance to smaller angles, since this increases the contribution of the non-s terms to 
the Bhabha cross-section. In the case of Bhabha scattering, the calorimctric measurement 
problem arises if a cross-section not inclusive of the energy of the outgoing fermion is 
considered. For thresholds of the order of 1 GeV the effect is of order 0.01 0.02%, but 
for higher energy thresholds (£'th) the contribution can grow considerably. 



3.6 Comparisons 

Let us now consider the question of de-convoluting the peak quantities in order to ex- 
tract pscudo-observables like cTq and A''". By definition the de-convoluted asymmetry 
does not include final-state QED radiation, while the de-convoluted cross-sections in- 
clude both QED and QCD final-state corrections. For ^/s = GcV the three codes 
predict a hadronic cross-section of 30.366 nb, 30.375 nb and 30.373 nb respectively. The 
corresponding de-convoluted quantities — no QED corrections — are 41.400 nb, 41.409 nb 
and 41.402 nb which, in turn, means that for BHM the effect of extracting QED corrections 
amounts to 11.034 nb, for ZFITTER 11.029 nb, and for TOPAZO, also 11.034 nb with a 0.045% 
difference. If we introduce 

a'^^Via'^] (l + 5conv) , (136) 
where T> denotes de-convolution, we get 

(5conv(B,T,Z) = -0.2665 , -0.2665 , -0.2664 . (137) 
From our previous tables we also find that for the same choice of input parameters 

(7^(B,T,Z) = 41.436 , 41.437, 41.441 nb. (138) 

The corresponding differences as compared to the dc-convolutcd obscrvables, which are 
0.036 nb, 0.028 nb and 0.039 nb, give an estimate of about 0.011 nb for the uncertainty 
in the effect of the imaginary parts from the weak formfactors. Coming now to the 
asymmetry we find: 



^pg(B,T,Z) = -0.00082 , 
at y/s — M^, becoming after de-convolution, 

VA^^{B,T,Z) = 0.0169 , 



-0.00125 , -0.00109 (139) 
0.0166 , 0.0166 . (140) 
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Thus the effect of de-convolution is 0.0177, 0.0179 and 0.0177 in BHM , TOPAZO and ZFITTER, 
with a BHM/ZFITTER-TOPAZO difference of 2.0 x 10"^. It must be noted that the quantity 
usually reported in the literature is the de-convoluted peak asymmetry with a pure Z 
exchange. In this case our predictions become 



r'Apg(B,T,Z) = 0.01544(0.01544) , 0.01536(0.01536) , 0.01528(0.01531) , (141) 

where in parenthesis we have included the corresponding prediction for Therefore 
this additional and conventional filtering of the asymmetry brings the BHM-TOPAZO dif- 
ference to 8.5 X 10"^ and the ZFITTER-TOPAZO difference to -8.0 x 10"^. The various 
combinations of de-convoluted asymmetries are presented in Table |T^. The complete set 
of de-convoluted quantities is presented in Table |15|, where, around the peak, we see the 
largest difference among codes of 0.04% and of 0.02% for cr'^, . For the muonic cross- 
section we find at a/s = 1.4785 nb, 1.4790 nb and 1.4794 nb, while the corresponding 
de-convoluted quantities are 2.0015 nb, 2.0019 nb and 2.0022 nb respectively. Thus the 
effect of de-convolution in is -0.2613,-0.2611 and -0.2611 for BHM, TOPAZO and 
ZFITTER. An interesting question, which arises in this contest, is related to the possibil- 
ity of performing a real and significant test of the QED corrections by comparing the 
predictions of different codes. We could define 

iO') = - 1 (142) 

QED ^ ' ^ ' 

for a given observable O, with Oq being the de-convoluted observable as predicted by 
the i-th code. For a given code this 5^^^ represents the effect of the convolution over an 
electroweak-corrected observable. It must however be noted that the comparison of 5^^^ 
of different codes does not give an unambiguous information on QED corrections since we 
start already from slightly different kernels. Let us quantify this statement. Let A^^^^^ be 
the absolute QED correction for some observable O as computed by the z-th code: 

0^ = 01 + A;^^ . (143) 



Then 



A 



'^qed(O0 = -f^ , (144) 
so that the difference between 5^^^ of different codes can be written as 

A* A^' 

- 5qed(0^) = ^ - ^ • (145) 
By adding and subtracting A^^^^^/Og one obtains 

5 {0')-5 fO^) = A* — — A A' - ) (U6) 

"qEdI*-^ J ^QEdK'-^ J ^QED ninJ QED ^QED^' ■ K^^"^) 
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Now it is clear that the difFerence between 5^^,^ of different codes depends also on the 
differences in the pure weak and QCD libraries of the codes. In particular, for the leptonic 
asymmetries in the Z peak region, which are small, the first term in the r.h.s. of the last 
equation becomes very large, so that 5qed(0*) ~ '^qed(^'') weak- and QCD- dominated, 
whereas it becomes QED-dominated far from the peak only. 

Coming back to our original strategy we have made a constant effort to understand 
the systematics inherent in the extraction of pseudo-observables from the realistic dis- 
tributions that the experiments should consider as a theoretical uncertainty. The main 
ingredient contributing to the systematics is of course the QED radiation, inclusive of 
initial-final interference and of initial-state pair-production. But in the extraction, some 
relevance must also be attributed to the imaginary parts of the formfactors and to the 
Z — 7 interference. The latter may have some influence, since not all the codes have 
the same splitting of the 7// vertices among the formfactors. We have already devoted 
some detailed discussion to the differences among the convolutions around the peak. In 
Table |T6|, we give all the remaining results, corresponding to the full set of the energy 
points. With due caution as to the correct interpretation of our comparison we observe 
a somehow larger difference in the convolution around the tails with respect to the peak 
region, notably 0.23% at = 88.45 GeV and 0.59% at = 93.70 GeV. Since these are 
differences in b^omi and not in the total prediction, even 0.6% away from the peak is quite 
reasonable. 

Our results are presented in Figures 24 - 37 where we have reported and A^^ for 
four different set-ups: fully extrapolated, s' > 0.5 s and 40° < 9- < 140°,i?th > 20 Gev, 
6'acoii < 10°, 25°. Also reported are for two set-ups, fully extrapolated and s' > 0.01s 
and the s channel a"^, A'^ where, however, i?th > 1 GeV. 

In all these figures we have also shown the deviations (relative for cross-sections and 
absolute for the asymmetry) BHM-TOPAZO (when available) and ZFITTER-TOPAZO with 
the corresponding theoretical error bars. The choice of TOPAZO as a reference point is 
purely technical and avoids the necessity of introducing an average among the codes. It 
emerges from these comparisons that for a fully extrapolated set-up, the agreement for 
the muonic cross-section for energies below the peak and around it is quite reasonable, 
always taking the intrinsic error as a reference. Given the reasonable agreement at the 
level of de-convoluted cross-sections and once we have observed that the de-convolution 
is satisfactory for hadrons, we come to the conclusion that for muons the low-g^ region, 
where mass effects may become relevant, gives the dominant difference in cr'^. Indeed, a 
comparison for s' > 0.5 s shows a much better agreement, giving 1.4391 nb, 1.4396 nb and 
1.4397 nb for the three codes (instead of 1.4785, 1.4790, 1.4794). Since the corresponding 
de-convoluted quantities are 1.9599 nb for BHM and 1.9605 nb for TOPAZO we conclude that 
de-convolutions amount to —0.2657 for both codes, while for fully extrapolated they are 
—0.2613 versus —0.2611 with a 0.08% of relative difference. Around the peak, the energy- 
dependence of the observables as predicted by BHM looks different — a difference which 
should probably be related to the way in which BHM implements the effective coupling 
language. For the hadronic cross-section we observe a consistent agreement among the 
three codes. Finally, for the muonic asymmetry the agreement at the peak is again quite 
reasonable, but the energy dependence again looks different, with BHM/TOPAZO agreeing 
below the peak and TOPAZO/ZFITTER agreeing above it. To summarize the status of 
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the comparison for extrapolated set-up or for s'-cut we have from Figures 24 — 29 the 
following. 

• a^' 

Around the peak the maximum deviation is 9 x 10^'^ nb corresponding to 0.06%, 
the BHM predictions are always lower, those of ZFITTER higher, and TOPAZO is in- 
between. On the high energy side the BHM behavior with energy tends to differ. For 
an s'-cut of 0.5 s the maximum deviation around the peak corresponds to 0.04% and 
on the low-cncrgy side of the resonance we observe a rather remarkable agreement 
among all codes. 

• a'' 

For the hadronic cross-section we have found an impressive agreement around the 
peak, with 0.03% of maximum relative deviation, which is only slightly worse around 
the tails — 0.08% and 0.06%. The comparison remains substantially unchanged for 
a low s' cut. 

FB 

Here the agreement is quite reasonable at the peak, with an absolute deviation 
of 4.3 X 10^'^ between BHM and TOPAZO. Not completely satisfactory is the energy 
dependence, which registers a substantial disagreement with ZFITTER giving higher 
predictions below the peak and BHM after it. At least on the low-energy side the 
situation improves if an s' = 0.5 s cut is imposed. 

The agreement between TOPAZO and ZFITTER remains rather remarkable even when the 
geometrical acceptance is constrained, and as well, final-state energies and the acoUinear- 
ity angle are bounded with or without QED initial-final interference. For instance, 
for the s-channel leptonic cross-sections at ^/s — and 40° < 9- < 140°, £'th(A*) > 
20 GeV, Eth(e) > 1 GeV, we find 





6'acoll 


< 10° 


0.9802 


- 0.9801 nb 




^acoU 


< 25° 


0.9905 


- 0.9900 nb 


^^ 


^acoU 


< 10° 


0.9886 


- 0.9884 nb 




^acoU 


< 25° 


1.0012 


- 1.0011 nb , 



(147) 



where the first entry is ZFITTER and the second TOPAZO. More generally, and remembering 
that we use E'th > 20 GeV for muons and a lower cut of 1 GeV for electrons, we find from 
Figures 30 - 37: 

for the muonic cross-section there is very good agreement for all energies and 6'acoii < 
10° — agreement which slightly deteriorates for ^acoii < 25°. At the peak we have 
a relative difference of 0.02%, 0.05% respectively for two above mentioned ^acoU- 
(0.03% for a fully extrapolated set-up). For s-channci electrons the agreement is 
everywhere of the same quality: in particular at the peak we find a 0.02%, 0.01% of 
relative TOPAZO/ZFITTER deviation — for 6'acoii < 10, 25°. Thus for cr' our agreement 
is not altered by introducing cuts. 



65 



The two leptonic forward-backward asymmetries agree at the peak at the level 
of 0.3,4.3 X 10"^ for muons (^acou < 10,25°) and 2.9,0.1 x 10"^ for electrons. For 
6'acoii < 10° the agreement tends to deteriorate at larger energies, reaching 9.6 x 10~^ 
for muons and 8.7 x 10~^ for electrons at ^/s = 93.70 GeV, while remaining always 
very good for Oacoii < 25°. Globally for the TOPAZO-ZFITTER comparison shows 
for the fully extrapolated set-up a larger difference on the low-energy side of the 
resonance (1.0 x 10"'^ at ^/s = 88.45 GeV), a maximal agreement over the whole 
range of energies for Oacoii < 25°, and a larger difference on the high-energy side for 
^acoii < 10° (9.6 X 10-^ at = 93.70 GeV). 

We illustrate the electroweak theoretical error by considering at ^/s = . In the four 
different set-ups considered — fully extrapolated, s'-cut of 0.5 or 40° < 9- < 140°, Eth > 
20 GeV, and ^acou < 10° or < 25° — we find 0.095%, 0.097%, 0.066% and 0.068%, respec- 
tively, in the relative deviation between the maximum and minimum predictions among 
the codes. 

Finally, we illustrate the effect of different treatments of final state QED radiation in 
the presence of severe kinematical cuts. By adopting different strategies TOPAZO predicts 
for ct'^ at = and 40° < e_ < 140°, > 20 GeV, ^acou < 10°: 

= 0.9801 - 0.9817 nb , (148) 
while for a loose cut of Eth > 1 GeV we obtain 



= 0.9892 - 0.9893 nb . (149) 

As already discussed, the 0.16% difference at large Eth reduces to a mere 0.01% at low 
Eth- 

Coming now to the full Bhabha cross-section and forward-backward asymmetry, we 
have considered in the following a comparison between the results of TOPAZO and of 
ALIBABA. As already explained, the comparison is not fully consistent as it stands now, be- 
cause the electroweak and QCD libraries of ALIBABA are not up-to-date. Nevetheless, we 
show it because it gives an impression of the state-of-the-art. For this particular compar- 
ison the input parameters are slightly different from our default — namely rrit = 174 GeV 
and = 0.124 have been used. The ALIBABA-TOPAZO comparison for is shown in 
Tables [T^ for 40° < 6'_ < 140°, Eth > 1 GeVand6'acoii < 10°. For TOPAZO we have shown 
three sets of numbers, all including QED interference, with the first two showing the ef- 
fect of a different treatment of QED final-state radiation — I is the TOPAZO default, while 
II is the ZFITTER-like default, and III shows the effect of initial state pair production, 
which is, however, not included in ALIBABA. It should be mentioned at this point that 
the insertion of pair production is strictly valid in TOPAZO only for s-channel processes 
and that it is approximate for full Bhabha. The cross-section is shown for ALIBABA with 
its numerical error, while for TOPAZO we give first the electroweak theoretical error and 
than the numerical one. Expressing due caution in comparing the two codes we observe 
a relative difference, of 0.05% at the low-energy side of the resonance, which becomes 
0.003 -r- 0.004% around it, with a visible deterioration at high energies where we reach 
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0.86%. As already discussed ||^, this is mainly due to a different treatment of higher-order 
QED final-state corrections. 

A more detailed comparison between ALIBABA and TOPAZO is shown in Tables WE 



25 



We have shown the predictions for the full Bhabha cross-section and the forward- 

we 



backward asymmetry for two different values of 6'acoii in Tables p!^-|2(]|. In Tables |2T| - p2 



give the relative deviation between the central values of ALIBABA-TOPAZO for full Bhabha 
or for the s-channel alone. Next in Table 23 we give the difference between full Bhabha 
and s-channel results both in ALIBABA and in TOPAZO with the relative contribution, i.e. 
6 = a /a{s) — 1 and 5{A) / 5{T). In Table ^ we show the s + t, s and t forward-backward 
asymmetries. Finally in Table ^ we give a rough estimate of the theoretical uncertainty 
by considering the difference between maximal and minimal predictions from the two 
codes. For ALIBABA this takes into account the numerical error alone while for TOPAZO 
we have added linearly the electroweak and the numerical uncertainty. With due caution 
in the interpretation we extract a 0.1 -i- 0.2% before and around the Z resonance which 
becomes as large as 0.9 1.0% at higher energies. 
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4 Basic Formulae for Electroweak Radiative Correc- 
tions 

In this section we give a more detailed description of the realizations of the effective 
couphngs and g^. 

4.1 BHM/WOH basics 

4.1.1 Self-energies, propagators, and Ar 

The radiative corrections to the photon-Z propagator system (considering only the trans- 
verse parts ~ Qfj^iy) can be obtained by inversion of the matrix 

with the renormahzed self energies specified below, yielding 



where (s = /c^) 



1 



^""^ ^""^ [s + S"^ (s)] [s - M2 + S^^ (s)] - [S^^ (s)]2 ■ ^^^^^ 

The building blocks of Eq. (|150|) are the renormalized self-energies S, which are decom- 
posed into unrenormalized ones S and counter terms, as follows: 



t'\e) = S""(A;2) - 5Ml + 5Zl - M]) , 
S-'^(fc^) = s'^-(fc^)-5<+5Zr(A:^-<), 

f^i^) = - 6z;' e + {6z;' - 5zi') mi . (153) 

In the last line the abbreviations (i = 1, 2) 

5zT = -^^{szt-szD 
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and 



W W ' z ^ w w 

were used. 

The self-energies S*-' in ( |153|) are the sum of the electroweak one- loop diagrams 0, 
completed in the quark loops by the 0{aas) two-loop QCD-electroweak contributions. 

The mass counter terms for W and Z follow from the on-shell conditions for the W 
and Z propagators: 



w 



SMt 



7^e (M^) ^ 



(154) 



The other renormalization constants in ( |15^j| ) are given by the following set of equa- 
tions: 



5z; = -n\o) ^ -f^(o) 



6zl = -n"(o) 



Ml 



5z; = -n^(o) - 2 ^- ^ ^ 



5z^ = -n'(o) 



s^c^ Ml 
SwCw Ml 




< = -n^(o)-^^5^ + :3^(^-^|. (155) 

WW z 

The photon vacuum polarization 11^ (0) contains as its light fermion part the quantity 
Aa discussed in section 1.2. The last two constants 5Z^ are not independent but are 
linear combinations of 5Zj and 5Z^ . They are given here for completeness. 

By the presence of the (S^^)^ term on the r.h.s. the equations (|155|) are non-linear 
equations. It is, however, straightforward to solve them for the renormalization constants 
in terms of the unrenormalized quantities. 

The higher-order irreducible contributions to the p-parameter, as far as available, are 
built in by means of substituting 

1#-1#-1#-1# + Ap("0> (156) 

z w z w 
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for the r.h.s of (|155D, with 



Ap 



(HO) 



3xt l + xtAp^^\0 



rQCD 

^(3) 



(157) 



a 



Xt 



WW z 



87r2v^' 



M2 



(158) 



comprising the two-loop electroweak and three-loop QCD contributions. For the functions 
Ap(2) and see section 1.10.2 and 2.3. 

Around the Z peak, the Z propagator has the form 



i + n''(M2; 



e _ /Vf 2 I ^ _s_Y 



(159) 



with 



n (m: 



ds 



S (s) 



s + S^Vs) 



(160) 



The Z width is calculated from the effective coupling constants given below in Eq. 
( |175| ) together with the general formulae of section 1.10.1. 



The vector boson masses , are correlated by the Fermi constant 

vra 1 na 1 

G,j 



v^M2 s2 1 - Ar V2M2c2 s2 1 - Ar ' 

^ w w ^ z w w 

The quantity Ar(a, M^, M^,mt, M^) has the following representation: 



(161) 



Ar 




S""(0)-5M, 



M2 

w 

7-4^2 \ 

^^logcM. (162) 
w / 

The last term is the sum of the box contributions and renormalized vertex corrections 
to the muon decay amplitude after removing the Fermi-model-like virtual photonic correc- 
tions. Due to the inclusion of the higher-order reducible and irreducible terms the way of 
writing Ar in the denominator of ( |161| ) automatically takes account of the proper resum- 
mation of the leading Aa and Ap terms and some of the sub-leading terms as discussed 
in sub- section 1.10.2. 
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4.1.2 Vertex corrections 



The vertex corrections can be summarized in terms of s-dependent vector and axial vector 
form factors if the masses m/ of the external fermions are small compared to M^, both 
for the electromagnetic and the weak NC vertex. In our terminology, 'vertex corrections' 
denote the renormalized 'y{Z)ff three-point functions in one-loop order, together with 
the finite wave function renormalizations for external fermions. 

In contrast to the propagator corrections, the vertex corrections are not universal 
and depend on the fermion species. For this reason we have to list them separately for 
u, e, u, andd type fermions. In addition, the h quark is exceptional due to the virtual top 
contributions in the vertex. 

Our terminology is as follows. 

Fy^j^ and -Fy'^ denote the IR finite weak (without the virtual photon diagrams) form fac- 
tors for the Zff and 7// vertex which, together with the lowest order terms, yield the 
dressed vertices: 



fZff 



7m 
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+ Flf{s)\] 



riff 



^if ( 



(163) 



The lowest order coupling reads: 

= if - 2Qfsl , aj = if , (164) 

and the weak form factors in (|163|) are explicitly given by the following set of formulae. 

Neutral current vertex: 

neutrinos: 



^ A 

a 



(165) 



W W 



A2(s,M,) + ^'- 



4s2 

w 



-A2(.,M^) + ^A3(.,M^; 



charged fermions: 



a 

a 



Air 



Vf{vj + 3a}) 

w w 

af{3vj + aj) 

w w 



A2(5,MJ + F[ 



A2(.,MJ + F[ 



(166) 
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with 



1 3c^ 

-^A,{s,Mj-^As{s,Mj, 

w w 

1 - -s^ 3c^ 
-f^A2is, + -f As{s, MJ , 

w w 



3d 



2s2 

w 



(167) 



Electromagnetic vertex: 



with 



a 

4:71 

a 
An 



w w 



(168) 



Gl 



4s2 

w 

1 



12s2, 



A3(s,M^) , 



A2(s,M^) + ^A3(s,M^) , 



65^, 



A2(s,M^ 



4s2 



A3(s,M^ 



(169) 



The functions A2, A3 have the form, 



A2(s,M) 



A3(s,M) 



with 



-- - 2w - {2w + 3) log(-w;) 



Li2 1 + 



+2(1 + w;) 



5 2w {2w + 1) 

6 ~ T ^ 3 



TT 



Vl — 4w log(x) + -w{w + 2) log 



(170) 



VI - - 1 

W = , X = , : . 

s + ze VI - 4w + 1 



The functions and Gf^ cannot be used for b quarks. The full expressions for Fl, G 
can be found, for example, in the second of Ref. 0. 

^•^ Note that the normahzation is different: 

F^(this report) = 2s i^^(Ref. [|). 
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4.1.3 e+e 



ff amplitudes 



Around the Z resonance, the amphtude for e^e — ff can be cast into a form close to 
the lowest order amplitude: 

/l(e+e- -> ff) = + + {hox) , (171) 

where denotes the dressed photon, Az the dressed Z exchange amplitude, and (6ox) 
summarizes the terms from the massive box diagrams, which, however, can be neglected 
around the Z . 

The dressed photon exchange amplitude can be written in the following way: 



A^ 



QeQ 



f 



[(1 + F^ll, - FTW] ® f(l + Fy)r - Fl^Yl^ 



l + W{s) s 

iT is the subtracted vacuum polarization n^(s) = n^(s) — n^(0) with 



(172) 



The vertex form factors Fy^^j^ from Eq. ( |168|) are evaluated for s = M^. 
The Z exchange amplitude without the box diagrams factorizes as follows: 
A, = V2G,Ml{p,pff'^ 



X 



^75 



■(173) 



The weak corrections appear in terms of fermion- dependent form factors p/ and /t/ in the 
coupling constants and in the width in the denominator. 



4.1.4 Effective neutral current couplings 



The factorized amplitude ( |173| ) allows us to define NC vertices at the Z resonance with 
effective coupling constants g^^^, synonymously with the use of Pf, ^f-. 



jNC 



V2G^Mlpf 



lfj(3) ori„„2 ,_\^, r{3). 



The effective couplings read as follows: 



1/2 



1 - Ar 



l + n^(M2) 



1/2 



9i 



+ f!^ 



1 - Ar 



l + n''(M, 



(174) 



(175) 



The building blocks are: 
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• Ar from Eq. ([TB2D 

• Il''{Ml) from Eq. (|160D 

• U"iMl) = t''{Ml)/Ml from Eq. (HI) 

• Fy^j^ from Eq. ( [166D for s = M^. 

For given mt,M^ the values of respectively are chosen such that Eq. ( 161 ) is 

fulfilled. 

For the b quark couplings the next-order leading corrections ~ Cj^mf, ~ agG^mf are 
taken into account by performing in the one-loop expression F| the following substitution: 



^ IQttsI.cI, Ml 



with 



An 



(176) 



(177) 



The function r(2) is taken from Refs. 



], and the QCD correction term from 



Ref. 57 



The alternative form factors p and k can then be obtained via 

7/ 



^ = l-A\Qf\KfS 



pf 



(178) 



Due to the imaginary parts of the self energies and vertices, the form factors and the 
effective couplings, respectively, are complex quantities. The effective mixing angles are 
calculated from the real parts according to 



sin^ eL 



efT 



4|Q/ 



1 



neg{: 



(179) 



4.2 LEPTOP basics 

4.2.1 Electroweak loops for hadron-free observables: functions Vi 

For hadron-free observables we write the result of one-loop electroweak calculations in the 
form suggested in ref. |7^ : 



c + 



3c 



327rs2(c2 - s2 

1 3 

2 ~ 647rs2c2 



-aVra{t, h) 



R 



il9\ 



+ 



47r(c2 — s2 



1 



2 647rs2c2 



aV,^(t, h) 



(180) 
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where t = mf/M'^, h = M^/M^ and functions Vi{t, h) are normalized by the condition, 

Vi{t,h)c^t, (181) 



at t > 1. 

Each function ^ is a sum of five terms 



Viit, h)=t + T,{t) + H,{h) +C, + 5Vi{t, h) . 



(182) 



The functions t + Tj(t) are due to (t, 6) doublet contribution to self-energies of the 
vector bosons, Hi{h) is due to W^, Z and H loops, the constants Cj include light fermion 
contribution both to self-energies, vertex and box diagrams. 

To give explicit expressions for Tj(t) and Hi{h) it is convenient to introduce three 
auxiliary functions Ft{t), Fh{h) and Fl^{h) (see subsection 4.2.5 for their expressions). 



The equations for Tj(t) and Hi{h) have the form |]72| : 



m 



3 9 

2 C2 - S2 



+ 3- 



4 
' 3 ' 

1C2 



32 o 2 , o 



-t 



2 16 2 2 32 2 , X 
- -t s^t Ft t , 

3 9 3 9 ' ^ ' 



/i , , c^h h s2 8 2 
In + m — -n s 

h-1 h-c^ c2 18c2 3 



h 



h — (?• 
4:h 4' 



h 



where = 0.23117 - 



9c^ 9c2 ■ 3^ 

i = A 



9 



Ah 4 



9 3^ 
(1.1205- 2.59(5^2) , 



(183) 



HAih) 



3 9 
'32 



16 _ 
27' 



1 - 2tFt{t) 



At-1 



+ 



4 "2 
— s s 

9 3 2 



-tFt m - - 1 + 2t ^ 

3 ^ ^ 3^ ' 4t- 1 



1,4 



8/i 



l-cV/."-c^ 9(/.-l)^^^^+G-^+^0^^(') 
^ - + ^/^') - + (0.7751 + imSs') . 

i = R 



(184) 



HR{h) 



^lnt + ^-^(l + lW) 

4 /i c2 , /i /4 4, 1,2 
\ In— ' ' u ^ 

3 18 l-c2//i c2 



+ (3-9^+9^ '^'^^^^ 



ln/i + (1.3590 + 0.5155^ 



(185) 
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I = V 



(186) 



The constants Cj are rather comphcated functions of sin^ 9 and we present their numerical 
values near = 0.23117: 



Cm — 


-1.3500 + 4.13 5s^ 


Ca = 


-2.2619 - 2.63 5s^ 


Cr = 


-3.5041 - 5.72 5s^ 


a = 


-1.1638 -4.88 5s^ 



(187) 
(188) 
(189) 
(190) 



4.2.2 Corrections 5Vj 



Functions 5Vi{t, h) in Eq. ( |182|) are small corrections to Vi. They can be presented as a 



sum of five terms 5kVi{k 



5): 



1. Corrections to polarization of electromagnetic vacuum due to W boson loop, S^a, 
and t-quark loop, Sta, are traditionally not included in the running of ^(g^). We 
also prefer to consider them together with electroweak corrections. This is especially 
reasonable because W contribution S^^a is gauge-dependent, while 6at is negligibly 
small. Here and for all other electroweak corrections we use the 't Hooft-Feynman 
gauge. The corrections 5^ a and 6ta were neglected in Ref. [|72| and were introduced 



in Ref. 



16 1 

5iVmit,h) = — —Trs^-{S^a + Sta) = 
16 1 

SiVR{t,h) = — -7is'^c^-{5^a + 5ta) 
o a 

6iVA{t,h) = 6iV,{t,h) = 0, 



-0.055 



-0.181 



(191) 

(192) 
(193) 



where 



a 
a 



1 

2^ 



(3 + 4c' 



' 1 - ^/Ac^ 



4 



;i + 2t)F,(t)-- 



1 arcsin ■ 
4 1 



2c 



457r t 



= 0.0686 
-0.00768 . 



(194) 
(195) 



(Here and in Eqs. (196 - 20£) we use = 175 GeV for numerical estimates.) 



76 



2. Corrections of the order of acts, due to the gluon exchange in the quark electroweak 
loops (see also Ref. ||73|), 62Vi = + d^Vi. For the two generations of light 
quarks (g = u, d, s, c) this gives: 



sivut, h) 

SlVA{t,h) 
S^Vnit, h) 



[4 




13 








l3 


TT 



2 2 , 20 4 
C — S + —s 



IT 



TX 



-0.377) , (196) 



750), (197) 



0, 

5lVA{t,h) . 



(198) 
(199) 



The results of calculations for the third generation are given by rather complicated 
functions of the top mass and s^. Here we present an expansion for large values of 
t and for fixed value of = 0.23117: 



TT 

as{mt) 



-2.86 t + 0.46 Int- 1.540- 



0.68 0.21 



t 



SlVA{t,h) 



TT 
TT 

TT 
TT 

TT 



(-11.67) 



-2.86 t + 0.493 
10.10) , 



0.19 0.05 



t 



-2.86 t + 0.22 Int- 1.513- 



0.42 0.08 



t 



t2 



-11.88) 



SlV,{t,h) = 6iVA{t,h). 



(200) 

(201) 
(202) 

(203) 
(204) 



As these formulas are valid for rrit > M^, in order to go to the region rrit < we 
either put (JjVi = or use a massless limit in which S^Vi = ^^l^i. In any case, this 
region gives a tiny contribution to the global fit. 



3. Corrections of the order of aa^ were calculated for the leading term aa^t only |K 

SsViit, h) ~ -(2.1552 - 0.18094 Nf)al{mt)t ^ -1.250 al{mt)t = -0.06 (205) 

for Nf = 5 light fiavours. [For the numerical estimate we use dg(M^) = 0.125.] 

4. The leading correction of the order a^t^, which originates from the second-order 



Yukawa interaction, was calculated in Refs. [p5[]-||58||: 
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S4Vi{t,h) 



a 



167rs2c2 I t / 



(206) 



where function A{Mjj/mt) can be found in Refs. |^,||5^. For = 175 GeV and 
= 300 GeV one has A = 8.9 and S^Viit, h) = -0.11. 

5. In the second order in electroweak interactions quadratic dependence appears on 



the Higgs mass [|75 



S^Va 
S^Vr 



a 



247r V c2 



0.747 = 0.0011 



a I h\ 

— 1.199 = 0.0057 



247r \s^^ 

a ( h \ — 



247r \ c2 



0.973 



-0.0032 



(207) 
(208) 
(209) 



(Here for numerical estimates we use mn = 300 GeV.) 
4.2.3 Hadronic decays of Z boson 

For the partial width of the Z decay into a pair of quarks qq (g = m, s, c, 6), we use the 
equation 



F, = 12 



{glfRl + {glfW Fo, 



where the final state QED and QCD radiative functions i?^ and i?^ are given by Eqs. 

The electroweak radiative 



(pq) and ( pq ) respectively, and Fq is defined by Eq. 
corrections are included in and g'^ : 



qi = J(3) 



3a 



327rs2c2 



gl/gl = l-A\Q,\s^ + 



47r(c2 - 



-aVRq{t,h) 



(210) 
(211) 



The functions VAqit, h) and Vnqit^ h) in the one- loop electroweak approximation are related 
to the functions VAit, h) and VR(t, h) from leptonic decays ||76| : 



VA«(t, h) 
VAdit, h) 
VruH, h) 



K4c(t, /i) = VA{t, h) + ^-^(F^, + Fau) 

128715^0^ 

V4,(t, /l) = VA{t, h) + — {Fai - FAd) 



VRc{t, h) = VRit,h) 



3a 

16iTsc(c^ — 



3a 



(212) 

(213) 
(214) 



X [Fvi - (1 - As^)Fai + ^(-(1 - ^s')Fau + Fvu)] 
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Vnd{t,h) = Vns{t,h) = Vn{t,h) + ^^^^^^^^—^ (215) 

x[Fvi-{l- As^)Fm + 3((1 - ^s^)FAd - Fvd)] , (216) 

where Fyi = 0.00197 , F^i = 0.00186 , Fy„ = -0.00169 , Fau = -0.00165 , Fyd = 
0.00138 , Faci = 0.00137 ||76[. The difference Vi^g — is due to different electroweak 
corrections to the vertices Zqq and Zll (s^= 0.23117 is assumed: formulas used in the 
code have exphcit dependence). 

The obhque corrections of the order of as, a'^t and at^ to the VAqiVRg) are the same 
as in the case of Va and Vr. But for Z boson decay into pair qq there are additional 
as corrections to the vertices that have not yet been calculated. This brings additional 
uncertainty into the theoretical accuracy. 

4.2.4 Specific features of the decay Z ^ bb 

For Z ^ bb decays we have to take into account corrections to the Z bb vertex which 
depend on t [§3. 

VM{t,h) = VAd{t,h)-^^^^^^m+sm , (217) 

Vm{t,h) = Vn,{t,h)-^^-^^[<p{t) + 6<P{t)] ■ (218) 
For (f){t) we use the following expansion 



0(t) 



3 -2s' 
2s2c2 



t + c^ 



2.88/n ( — 



6.716 + -( 8.368 In f 4 ) - 3.408 



+ ^( 9.126 In (X) + 2.26 c^] + ^ ( 4.043 In f 4 ) + 7.41 c' 
+ ... 



(219) 



and for 6(f)(t) we use the leading approximation calculated in Refs. [0 and |5^, 



5<p{t) 



3 -2s' 
2s'c' 



as[mt) 



TC 



t + 



^\ .2 (2) ( 

IGs'c'KnJ ' [t 



where function r^^^^ has been calculated in Ref. 
For = 175 GeV and = 300 GeV 



(220) 



.(2) 



1.245 



Asymmetries are calculated with the loop corrected values of g^, and 
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4.2.5 Appendix: Auxiliary functions Ft and 



2(1 - V4t - 1 arcsin ^) , 4t > 1 
2(1- v^r^lnii^l) 



4t < 1 



:22ii 



4 J ' 



1 + (^ - f ) In /i - /i J| - 1 arctan J| - 1 , /i< 4 



/i > 4 



(222) 



-1 + ^ln /i+(3-/i) y 
-1 



h-l 
2 



In /i + (3 - /i)^/ arctan J | - 1 



4.3 TOPAZO basics 

The realization given in Ref. describes the couphng of the Z as: 



/i > 4 , 
/K4 . 



j(3) 



(223) 



-2Qfs' + 6gl + {lf +6g{)^,] . (224) 

Before giving the specific expression of the various quantities entering the previous equa- 
tion, we stress that our metric is such that a time-like momentum squared is negative. 
Next we decompose the unrenormalized vector boson self-energies as: 



















'^zz 


(p) 


Z~i 





16 7r2 



16 vr^c; 
g^se 



2„2 zz 



P 



16 TT^Ce 
16^ {p ) ' 



p 



(225) 



where 6 denotes the bare mixing angle. After a re-diagonalization in the neutral sector, 
which makes 5*^^(0) = in the = 1 gauge and replaces ^^'w^ ^33 ^3q with T,^^ + 
4r, S33 + 4r and + 2T with T = M^Bo{0, M^, M^) — M being the bare W mass — 
we consider the ultraviolet and infrared finite corrections to the //-decay 6^ and introduce 



(0) + —^0^ 



(226) 
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which has the virtue of being gauge invariant, a property not satisfied by S^^y(O) alone. 
The extra term induced by the diagonahzation is gauge invariant by construction. From 
now on we will denote the fermionic (bosonic) contributions to a given quantity A with 
the notation A^^^™{A^°^). With these quantities we define 



'^l^ = S° -WeE,3(Af^)+TCeE„(AfJ) + sVM^n^"(A/?), (227) 



where 



2 2 
S C 



(228) 



and n^(M^) = 7^e^^^(M^) - n^^(O). Note that E is ultraviolet finite. Next we can 
introduce as: 



1 

2 



1 - 



4 7ra(M^) 



(229) 



and give the definition of Sg^,Sg^, 



a 
An 

a 
An 



2Ff - '^VfAU, 



2Ff 

A 



ir(3) 



- 2g/As2 



2-'/ 



(230) 



where Vj — if^ — 2Qfs'^, and 



MlAU^ = ^e{(E;j''°^-E^-(M^)-E;-(M^)-M^E;3(M^) 

+ sVU'jMl)} . (231) 

Here /' stands for —df/dp^. The F refers to fiavour-dependent vertex corrections. 

• If no resummation of bosonic self-energies is performed we have 

P:-Pz, and As'^^-^^; (232) 

• otherwise has the same structure of with S^'^"'^ replaced by E^: 
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-L -r 



+ h.o. 

7^eS,^(M 



M2 

z 



(233) 



and 



MS 



(234) 



An^ is the residual Z wave function factor which obtains by writing the Z propagator 

Xz as 



TZex 



-1 

z 



2v^7r2" 



-iferm 



+ An, . 



7^e 



S,,(M^) + S,,(p2 



(235) 



The resummation operates on a gauge invariant quantity since it can be proved that 



Mi 



— s 



2 



bos 



— 



S^^- AX^,,(0)-7^eS,,(M: 



bos 



(236) 

and the sum of the terms in the last parenthesis is automatically gauge invariant. The 
same quantity, however, is not ultraviolet-finite and therefore has to be strictly understood 
in the MS sense at the scale /j, = M^. The higher-order terms in Eq. (|227|) are given by: 



h.o. 



G.M^ 



+ 



27r2 

TT 



1 asjmt) 
2 



TT 



TT 

'^3 



nit 



AS 



Iq 



10.55 



' ai,{mt) 



TC 



mi 



(237) 



where AS2 is the two-loop factor proportional to m'^/M'^, AS^g denotes the 0{aag) light 
quark contribution and the last factor is the three-loop correction computed for six active 
flavours. In the special case of the Z bb width we include the well known correction 
factor |5^, which modifies the vector and axial- vector couplings into 1 — 4/35^ + r and 
1 + r. 

To illustrate the internal structure of the renormalization procedure we use a well 
defined and generalizable example. The three bare parameters of the MSM Lagrangian 
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are related to experimental data by one-loop relations. One example of a solution concerns 
the bare mixing-angle 



sin^^ = 

S2 = 



a 



— 



(c2 - S2)3 



M2 

z 

- n^^(o)si . 



(238) 



For the Z ^ ff amphtude the self-energy corrections give rise to 



r 



if - 2QfVf + if-f' 



(239) 



Thus: 



Vf 



sin^. + ^ 1 ^e^..iMl) 



+ C»(a2) 



q; 



Atx — s' 



+ 



4^/ 



47r — 



+ 



(C2 

a(M^ 
47r 



— 



+ 0(q;=^) 



- 



(240) 



After this result we may perform a partial resummation 



Vf = s' + 



a{M,) 



-ibos 



+ 



(241) 



47r — 

As a final comment let us consider the term proportional to the square of the Z — 7 
transition in the propagators. It would contribute to with 



An. 



27r2 M2(p2 + M2) 



7^eE, fM^) + -^TeeS, 



27r2 M2 



, , 7^eS„ (M2) 



n 2 



+ ... 



, (242) 



where . . . indicates terms of O (p'^ + M^^ . Thus the additional term, being 0(p'^ + M^), 
does not contribute to the Z wave function renormalization factor. It should be noted 
that there is an easy dictionary of translation with other realizations, for instance: 



Pf = 



4p, {lf + 5giy , 

.3) 2Qjs' + 5g[ Sgl 



(243) 
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Once electroweak corrections are included in the formulation, TOPAZO implements initial 
state QED corrections by a convolution on the weakly corrected kernel distributions with 
a radiator function, or with structure functions — depending on the experimental set-up. 
Resummation of soft photon effects and hard photon emission up to is taken into 

account, and final state QED radiation with realistic cuts, QED initial-final interference, 
and initial state leptonic and hadronic pair production are also included. In the next two 
appendices further details are given for pure electroweak corrections. 



4.3.1 Appendix 1: The self-energies 

Starting from the decomposition of Eq. ( 225 ) for the vector boson self-energies we give 



their general expression in terms of two-point scalar form factors ||22[. In the following, 
the first argument, p^, is always left understood: 



n 
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3 

+ 4E 



12 52i(M^,M^) + 75o(M^,M^ 



4 1 

B{mi, nil) + 2 5(m„, m„) + - B{md, rud) 

2 13 
Biggl{- - 10 (M^, M^) + y 5o(M^, M^.) 

+ [B{mi,mi) + 2B{mu,mu) + B{md,md)]Btggr}-2MlBo{M^,M^ 



^33 + ^33 . 



WW WW 



+ 



WW WW 

25 



= -- 952i(M^, M^) + -5o(M^,M^)- 521 (M„M,) 
-5i(M„M,)-^5o(M„M,) 

+ \Y1 l^i^h mi) + B{m^, my) + 3 5(m„, m„) + 3 5(mrf, m^)] 

= -2MlBoiM^,Mj + ^MlB^{M,,M„) + ^MlBo{M,,M,) 
- ^ MlB,{M,,M„) - 1 MlB,{M,,M„) 



\ H [m-lBoimy, m^) + m^Bo{mi,mi) 



+?, mlBo^m^, m„) + 3 m^5o(mrf, m^) 
^ (Ml - Ml)B,iM„M^) + i (13M2 - 21Mj)5o(M„ M^) 

+ l{Ml- Ml)B,{M^^,M,) + i (5M^ - M^)5o(M^, MJ 

+ Y [{m^-ml)B,{my,mi)-mlBo{my,mi) + 3{ml-ml)B,{m^,md) 
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-3mlBo{mu,md)\ , 

<w = 2{Ml-Ml)[2B,{M„Mj + Bo{M„M^)]-2Ml[2B,{0,M^) + Bo{0,M^)], 
Kw = l-9B,,iM„Mj~9B,iM„M^) + lBoiM„Mj 

- 52i(M^,M,) - fii(M,,,M,) - lso(M,,,M,) 

+ 2^ [B2i{mi,,mi) + Bi{mi,,mi) + 3 52i(m„, m<i) + 3 5i(m„,md)] , 

9 

Kw = S[B,,{M^,MJ + B^{M^,MJ-B,,{0,MJ-B,{0,MJ] 

+ 2[Bo{0,Mj-Bo{M^,Mj] , (244) 

where B = 2 B21 — Bq and the sum is over the fermionic generations. The factor A is 
given by A = — 2/(n — 4) + 7^ — ln47r. The functions x{^)y Gn{y) are defined by 



nl 

„n-l 



Gn{y) = / rfxx"-Mn (x - y) . (245) 
Jo 

In terms of x (where — it) we have 

Bq = A — / (ix In X , 

JO 

1 /•! 

i?i = — A+ / dxxlnx, 

2 Jo 

= -A - / rfxxMnx , (246) 
3 JO 



and the corresponding integrals can be written in terms of the G„-functions, for which we 
write recursion formulae to be worked upwards or downwards, according to the magnitude 
of y. 

The 0{aa^) contributions to 11^^ . . . S,^^ are computed according to the formulation 
of Kniehl (see Ref. |5^). For instance, 



n??» = ^ 4 l'--^ + y^ir)] . (247) 

9 TT 



where r = —1/4 p /m^ and the functions X and Vi are explicitly given in Ref. 



4.3.2 Appendix 2: The Zff vertices 

In order to introduce the vertex contributions F^^ we first present the fermion wave 
function renormalization factors (/ 7^ 6): 
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r{3) 



2Qfs'] +1 



1 

32 



F c F 

^ z 8 



(248) 



where = A — In (M^ ) ~ With them we can write (/ 7^ 6): 



V 

F{ 



: t;/ ( t;^ + - j F/ + 4 c'f/ - cHf fI , 



Wflf + WS - 2 Jf ) (3 vj + F/ + 4 Jf c^F/ - c^jf )F3^ , (249) 



where again Vf = I^P — 2 Q/s^. The functions F/ are given in terms of three-point scalar 



formfactors |22]. For instance, 



C2,{mlMl,ml)--p'' 



Cn(m^, Ml,ml) + C2^{ml M^, m^) 
Cn (m^ , , m^) + C23 (m^ , , m^) 



1 



+ 



^3 



-6 C,,{Ml,ml Ml) - [CoiMl,ml + Cn(M^ 
- C23{Ml,ml,Ml)]+l + A-ln{Ml). 



(250) 

In TOPAZO the most general (arbitrary internal and external masses) two-, three- and 
four-point scalar functions are available ||79|. For hb final states the expression for vertices 
contains additional mj-dependent terms which can be found in Ref. . 



4.4 ZFITTER basics 

Here, we introduce explicit expressions for Ar and the weak form factors of Z decay and of 
fermion pair production process e+e" — ff p, k to order (Hoop, a) and order (21oop,««s) 
as used in Eqs.([71| ) -([73| ) and (^)-(0) of the main text. 

4.4.1 Muon life-time 

The virtual, non-photonic one-loop corrections to the muon life-time are: 



^^iioop,. ^ -^j-- l + 2yg^,ln^ + , ^ ^ \W{M^)-Z{M'- 



+ 



W(0) - W(M^ ) - -R(l + R) + — + \nR 



1-R 

Here and in the following sections, we use the abbreviations 



M2 

D W 

Z 



M2 

_ S_ 

w 



M2 

z 



..(251) 



(252) 
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The Ar was introduced in (E.8) and (F.3) of |Q. The t mass dependent terms may be 
found in the appendix of |]8l| . 



4.4.2 Partial widths of the Z boson 



The two form factors for each fermionic partial width of the Z boson are in one loop order 
and in the approximation of vanishing external fermion masses 



Hoop, a 
Pf 



1 + 



a 



47r(l - R) 



Z{Ml) + Z {Ml)-W{Q) 



^ 1 1 

+ -R{l + R)- — ~- 
8 ^ ' 2 A{ 



-—^\nR + Uf + 5p\,j^, (253) 



K 



Hoop, a 
/ 



1 + 



+ 



a 



R 



47r(l - R)\l-R 
l-Rf 



Z{Ml)-W{Ml) + M{Ml) 



R 



f 



(254) 



1 

2R 



l-6\Qf\{l-R) + 12Q}{l-Ry 

3 



-2|Q/|(l-i?) 



+ 2R 



+ 



2R 



(255) 



4.4.3 Auxiliary functions 

All W and Z boson self-energy functions are sums of bosonic and fermionic parts, e.g. 

W{<d) = Wb{0) + Wf{0) . (256) 
The bosonic parts are given in appendix A of I^Op'^ 

3 



bRjl + R) 17 5 

8 4 8i? 8 



9 3 

4 4R ' 1- R 



\nR+ ^""^ 



, Inr,,, , 
4(1 -r^) 

(257) 



157 23 1 



+ 



12R 12R^ 2 

_ _|_ _iL _ w 
'a A 



1 

12 ■ li 
1 



7 7 1 
2 ^ 12^ ^ 24/22 



Ini? 



, lnr„, + , 
24 / ^ \ 2 



17 2 
-2i? - — + — 



1 \^v.z(^, 



6 m 24i?2 



6 24 

W ' 



M2 

w 



-8R'- 



3AR 35 



18 ^ i? 



M2 

2 

2 ^ 12i?^ ^ 



(258) 



4 4 



One should eliminate there the (approximate) fermionic parts proportional to TrQ^ and Nf. 
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Z,{Ml 



51ni? n 



6 24 



+ 



6 3 24 



M2 



M2 



(259) 



-AW 



11 R 23 



+ 



9 18i? 



12i? 2/? V 



+ 



/ 3 3r^ 

L ^ 1 £ 5_ 

6 V 4 8 12 

17i? 1\ L^^(M2) 



Inr, 



6 



12 8 



+ 



1 5r^ r^^ 1 

2 ~ ^ ^ T2" ^ 2(r^ - 4) 



2 ' 



M„ 



(260) 



Function L^^^^(s) = L(-s; Ml^,MlJ is defined in (2.14) of Ref. g and and 
I„(Q^; M^, M|) in appendices C and D of Ref. |8^. The fermionic corrections are: 



Zf{s) 
Wf{s) 



2R 



f 



M2 



(1 + l3(-s; mj, mj) + ^Io(-s; mj, m; 



f={UJd) 
,2 



z 
— S 



m 



l3(-s;mi,mi) + —^Ii(-s;mi,m 



Z;(A^) = -EA'/-|^[l-'-X,^(-^^.r4,m})]+lj(l + „|) 



/ 



6i? 



^ln(r/i?) +rfR + (-- 



(261) 



, r^fR)M^ J^(-M\mlm^f) 



(262) 
(263) 



with Tf = m'j/M'^ and the color factor iV/. The 7Z- mixing function is: 



M{s) = 2j2Ni\Qf\vM-s-m},m}) 



(264) 



The vertex functions are in the hmit of vanishing fermion masses: 



- 2R^ - (3 + 2R^) H-R^) + 2(1 + R^f [Li2(l + Ry) - U^il) 



^^^ + 2R^ {R^+2):F,{s,MI) 



(265) 
(266) 



with 



R. 



7v 



MyV, 



R.. 



(267) 



88 



The additional t mass corrections to the Zbb vertex and the counter term are at the Z 
resonanceP^ 



\ - 32/(1 - V) 



In |ri| + 2ri?ln jrsl - ri? + (2 + i?) 



Fi(r)-Fi(0) 
2rR{\ - R) 



R 



F2{r) - F2{0)\ + -rR{2 + R)F2{r) 



(1 - 4R) 



1 + In IrsI - 4Fi(r) + ^^^(r) 



(268) 



4 ['dyl-{2 + R) fF2(r)-F2(0) 



RJo 



2i?ln IrsI + (In {r^l - 1) 



+ -(r - 2rR - 4i? - 4)Fi(r) + -(1 - r + 2i? + 2rR)F2{r 



with r = Tf and 
ri 



r(l + 2i?) 
" 6(1 -r) 



■(5r- 11) + 



4 

3r(r - 2) 
1 — r 



Inr 



- 1 



y(l -y) 
'^a = ^+(1-2/)^^ 
F,(r) = /,(l,r), 



^2 



(1 -2/)^ 
l-(l-y)^ 



/.(r,l) 



z = 1,2 



(269) 
(270) 



f2ia,b) 



In 



g - yO--y)/R 

a — b — y/R ay + b{l — y) 



(271) 



A completely analytical expression, valid at arbitrary s may be found in [^. In the next 
section, we will need in addition the photonic self energy function A{s), and Dz{s): 



A{s) 



34 



17 



34i? 35 4 InR ^j^2Kwi 



(272) 
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3 18 9R 12R s 
Kw{s)-neL^^{Ml) , l\R 



+ 



17i? 2 

+ 77 + 
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3 24i? 



-S + M2 



(lO - 5r, + r^) i?^ + H-r^fR^ 



~2A 



11 - 4r„ + r' 



+ {'^-r,)R, 



24s 



+ 



+ ^ 

2 6 24 



+ M2 



^Z,/(S) 



M2 



M2 - S 

z 



Zf{s)-neZf{M^ 



(273) 
(274) 



The net one loop finite t mass effect from the two vertices and the counter term is taken into account 
in ZFITTER by a variable VTB. 
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4.4.4 Form factors of the process e+e 



ff 



The virtual, non-photonic corrections to fermion pair production, including Bhabha scat- 
tering, may be described by four form factors per production channel as is indicated in 
Eq. (^). The contributions from the ZZ^ WW box diagrams are vanishingly small at the 
Z peak. Leaving them out^^, it is |p3|: 



Pe/(s) 



lloop,o 



1 + 



+ 



47r(l -R)\ ^ ^' ' V2 8i? 8 4 



T In r 



1 — r„ 



1-R R 



InR 



6s 



-2R + 



+ 4 (^e + Vf) 



Kzis) + -6pi,j\ , (275) 
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I Hoop, a 



1 + 



a 



R 
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+ 
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. Hoop, a 
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47r(l-i?) [1-R 
-2M{s)-^R-2RV,^{s)- {R + R, 

1-R 



9 

A{s) 



6s 



+ 



R 



{Ql + Q] 



RJ. 



m 16R 
1 

2 



(277) 



4.4.5 The 0{aas) corrections to electroweak observables 

Here we follow presentation of Ref. [^, based on the second Ref. ||52 



^Q,21oop,aOa 
^^21oop,aOs 



45Q_2 
4 m? 



QX(MM + Q?nr(M; 



(278) 



The interested reader may find expressions for WW, ZZ boxes in Ref. |Q and in the subroutine 
ROKANCC . .,t - s, -s, -t, . . .) of ZFITTER for the s channel {t > 0). For the t channel (used in Bhabha 
scattering) the call is ROKANCC . . , s, t, u, . . .) . 
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+ 



R 
4^ 



_ 1 (c^ -.Mnr(M2)-< 



(279) 



^^21oop,aas 



acts m| 



7^e< 



4s2 

w 
r2 



(280) 



Ap 



21oop,aa:s 



aas m| 

W 14/ 



+ 



7^e 



^_g!.2^,^^^^4^- 



1 - QVM: 



nr(M^)-nr(Q^) 



+ nr {Ml) - nr (g^) - v^ur {md - k {md 



(281) 



and, when Q"^ = M^, 



w w 



ml 



K + l) + 



,/ , dnr(g^) 

* l ^ — dQ^ — ^ ^ 



+ 4no 



Here, the two-loop functions are: 



1 f 55 

- j — - 26q; + 3 (1 + x„) (1 - Ga) - 2 [a (2x^ - 3^^, + 2) + 2x«] 

+ 2 (4q;^ - l) la + 4q; (2q; - 1) (Aa + 1) | , (283) 

- 1 -T - + 12q;^ + 3(1+ Xa) (l + 12q; + 4q;^) /« + 2 \2xa (Sa'^ - l) 
q;|^42 ^ / I \ / 

+ a (7x1 - 3xa + 7)]/2 - 2 (1 + 2a) (1 + 4a) I« - 4q; (1 + 4a)^ J„| , 



(284) 



nr(g^) 



if 71 

-<! 55 - — a - lOa^ - 8aG{xb) + 2 (g + 9a - 5a^) 



4q; [ 



+ 2 4xb (^a^ - q; - l) + a (5 - 4q;) /, 



+ 4 (a - 2) (a + 1)^ lb + 8a{a- 2) (a + 1) Jb L 



(285) 



rfnr(g^) 



43 



+ 18a + (10a + 3) (1 + x^) fa + 2 (5a - 2) x,,/^ 
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where 

la 

J rv 



- 8 (2a - 1) xlflixl) - 2 Ua^ + 1 ) I« - Sa' (1 + 4a) J 



43 



12a (2a + 3) - (24a2 + 26a - 3 ) (1 + x^) /„ 



- 2(l2a2 + 19a + 2)xafl + 8(1 + 4a) xlfl{xl) + 2 {%a^ - l) I„ 
+ 16a2(l + 4a) J, 



2a 



1 + 2a + VT+Ia ' 
F{l) + F{xl)-2F{x, 



1 



(1 + Xa) a 



[G{xl) - G(xJ] 



a 

lb 
Jb 



0" 



F(l)-F(xf.), 
-^G(xfe) , fb 



a 



TTa ' 

In Xg 

1 -^a 



(286) 



(287) 



F(a;) 



y 



Vl-^/, 



(q2) 

with C(3) = 1.0202... Further, 



an 




r „AFi 

an 




an 






mt=0 




mt=0 




mt=0 



n7 = hm U,{Q' 
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27r^ 



(28^ 



The n, (Q2) 

is a common asymptotic expression for functions ( |283| - |285"D for rut 



55 

— - 12C(3) + 31na, (289) 



(290) 



5 Description of options in different approaches 
5.1 BHM options 

In this subsection we describe the main options implemented in the version of BHM used 
for the present study to estimate the theoretical uncertainties. Out of the actual different 
user-accessible flags and their possible values, only the ones used in the present analysis 
are discussed. The rest of the flags and/or values either have a negligible effect on the 
predictions or exist only for testing purposes. 

• IRES=0,1,2 

The resummation of the one-loop one-particle irreducible contributions in Ar and, in 
general, in the whole set of self-energies is implemented in BHM in three different ways. 
Tsese can be considered representative of at least two rather different philosophies. 
The basic difference comes from the treatment of the (S^ )^ mixing term. If IRES=0 
, then the resummation comes from the resolution of the renormalization equations 
keeping that term (see for instance W.noUik in Ref. 0]). In this case, not only 
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leading top and Higgs terms are re-summed but also non-leading ones. This is 
the default working option. If IRES=1,2 , then the renormalization equations are 
solved keeping strictly one-loop contributions and the inclusion of higher-order terms 
coming from one-loop one-particle irreducible diagrams is explicitly done for Ar and 
for each self-energy. If IRES=1 , the prescription used to treat the top, Higgs and 
remainder terms follows the suggestions of Halzen, Kniehl and Stong , whereas 



if IRES=2 , then it follows the suggestions from S.Fanchiotti and A.Sirlin |^9|. Both 
differ in the detailed treatment of the Higgs and remainder terms. 

• ITW0=1,2 

This flag allows the choice of scale in vertex corrections. The dominant effect hap- 
pens in the b-vertex. If ITW0=1 , then the a(0) scale is used whereas, if ITW0=2, 
then the scale is used. The later is the default working option. 

• IFAC=0,1,2 

This flag allows the choice of different factorization schemes for the final-state cor- 
rections. If IFAC=0 then no factorization at all is applied for weak vertex, QED and 
QCD final-state corrections. That means that these three kinds of corrections are 
independently applied to the vacuum-polarization dressed amplitudes. This is the 
default working option. If IFAC=1, then QED corrections are applied on top of weak 
vertex ones, and if IFAC=2, then QCD corrections are also applied on top of weak 
vertex and QED ones. 

• IQCD=3,4 This flag allows the choice of different treatments of the QCD correc- 
tions to electroweak loops. Of its possible values, only two have been used in the 
present study: if IQCD=3, then the exact AFMT correction is implemented, whereas 
if IQCD=4, then the Sirlin's scale ^ = 0.248 is used. These two approaches have 
already been discussed elsewhere in the text. 

To summarize, BHM runs in 3 x 2x3 electroweak x2 QCD options. We have added 



Table 26, showing the effect of the working options of BHM on theoretical errors. 



5.2 LEPTOP options 

Several options to the preferred formulas used in LEPTOP have been chosen and variations 
of these formulae have been made. Usually each option consists of the addition of an 
extra term corresponding to a rough guess of the value of the uncalculated higher-order 
terms. We then make all possible combinations of these options — 2" in total, where n 
is the number of options. Among all these 2" combinations we locate those yielding the 
minimum and the maximum values of the observables and take as the estimate of the 
theoretical errors their deviations from the central values. 

Theoretical uncertainties come from two different sources: 1) not yet calculated Feyn- 
man diagrams, 2) not yet calculated terms in a given diagram. The first source is rep- 
resented in LEPTOP by the gluonic corrections to the electroweak vertices of light quarks, 
Fiq, where q = u, d, s, c; i = V, A. In this basis for options a crude estimate of 

these corrections is used: 
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= ^F,, . (291) 

71 

The second source is represented by non-leading terms of a) higgs and b) corrections 
to the tt loop in Z boson self energy and t — b loop in W boson self energy denoted by 

2 

2 

5V^ and 6Vi ^ respectively and c) by the tt vertex contribution to Z ^ bb decay denoted 
by 5(f)* . 

In ) and b) the missing non-leading terms are estimated by multiplying the 

corresponding leading terms by a factor 2/t, where t = {rrit/M^)'^. In the case c) the 
leading correction itself is so small that it is taken as a measure of uncertainty. Thus the 
LEPTOP uncertainties are: 

AVf = (mSVf = (2/0 (292) 
in accord with Eq. ( P06| ); 

Avfl = i2/t)6vfl = (2/t)(-1.25a2(mi)t) (293) 



in accord with Eq. ( PO 

3 - 2s2 TT 



A0°^ = -1.37^— ^-al{mt)t (294) 
in accord with the first term in Eq. ( P20| ); 

A*" = = ^^-f ' (295) 

in accord with the second term in Eq. ( P20| ), as A0"» is much larger then 50*^, the latter 
will be neglected; 

Ai^7 = SFt,' = ^F,, (296) 

71 



in accord with Eq. ( ^9 1| ) . 

Explicit calculations [|l], ^ give: 



Vu 



-0.00169 , Fau = -0.00165 , (297) 



Fvd = 0.00138 , FAd = 0.00137 . (298) 

Note that the uncertainties AF^^^ produce corresponding uncertainties in electroweak 
corrections 6Tg to the partial widths Z —> qq: 

6T, = 24scro2jf [(1 - 4\Qg\s^)Fvg + Fa,] . (299) 
With above numbers 

Ar„ = AEc = -1.9 (^^^ MeV ~ -0.08 MeV , (300) 
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Ar^ = Ar, = -2.0 (^^) MeV ~ -0.08 MeV . (301) 
So that the sum over four hght quarks is 

E^Ar, = -0.3 MeV . (302) 
Now the options of LEPTOP can be formulated: 

• Options 1,2 

Uncertainty tsVf given by Eq. (|292| ) is added to (option 1), or subtracted from 



(option 2) the functions Vi, % = A,R. 



Options 3,4 



Uncertainty AV^ ^ given by Eq.( |293| ) is added to (option 3), or subtracted from 
(option 4) the function V^, 2 = m, A, i?. 



Options 5,6 

Uncertainties Arq(g = m, d, s, c) given by Eqs. ( pOOj ), ( P01| ) are added to (option 5) 
or subtracted from (option 6) the partial widths Eg. 

Option 7 

Uncertainty A</)"» is added to the function </)(t) given by Eq. (|219| ). 



Theoretical uncertainties for observables, assuming = 300 GeV, rrit = 175 GeV, 



ac;{M^) = 0.125 are shown in Table |27. 



5.3 TOPAZO options 

In this subsection we describe the options implemented in TOPAZO version 2.0 for studying 
the theoretical uncertainties. A general comment is in order here. Some of the TOPAZO 
electroweak options have been originally designed to produce a conservative estimate of 
the uncertainty. If nothing stands against a certain option then we accept it, even if it 
goes against our own philosophy. 

• OUO.EQ. 'N' CY') 

This is the primordial option since it controls the partial resummation of bosonic 
self-energies. If OUO.EQ. 'N', then in and As^, we use: 



^ ^ ^ferm 

^bos 

— 



As' = (303) 
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i.e., all bosonic self-energies are expanded. Otherwise for OUO.EQ. 'Y' we use: 



M2 

z 

M^As'^ = TleE^^iMl) . (304) 



For OUO.EQ. 'Y' TOPAZO will automatically select 0U2 = 'N' and 0U3 = 'Y'. 

• OUl.EQ. 'N' CY') 

The default choice requires that in 



a 
An 

a 
An 



-2Ff - \vfAIl^ 



1 r(3) , 



2F{ - An, 



(305) 



everything must be expanded in terms of a = a{0) — 1/137.036... while for 
OUl.EQ. 'Y' the parameter expansion is selected as a = a{M^). 

• 0U2.EQ. 'N' CY') 

This option deals with the so-called problem of expansion. The default for TOPAZO is 
the expanded solution where, for instance, a Z partial width is computed according 
to 



— — - N^o 
6V2n " 

If - 2Q,f . 



} + ^ + 2vfSgi + 2lfSgi 



(306) 



where for the moment we assume that there is no final-state QED -|- QCD radiation , 
factor, whose treatment will be explained by one of the following options. If, instead, 
0U2.EQ. 'Y', then the electroweak corrected vector and axial- vector couplings are 
squared numerically. 

• 0U3.EQ. 'Y' ('NO 

Two different procedures are introduced for dealing with the physical Higgs contri- 
bution to the self-energies. If 0U3 . EQ . ' Y ' , then by working in a MS environment 
with a scale /i set to we extract from T,^°^ the physical Higgs contribution, S , 
and redefine 



-iferm 



(307) 
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where T,^ is subject to no additional expansion, such as leading or sub-leading 
behavior with M^. Of course, 0U3.EQ. 'N' leaves the Higgs contribution expanded 
as for any other bosonic contribution. This option reflects and partially illustrates 
one of the defining rules of TOPAZO — a certain reluctance to accept the isolation 
and a different treatment for something which can be considered as the leading part 
of some quantity only for extraordinary values of some of the parameters. 

0U4.EQ. 'N' CY') 

When we consider the mass corrections to the bb partial decay rate there will be 
something like 



-6 



M2 



2lf6g{ + {6g 



(308) 



In these mixed corrections there is an additional uncertainty connected with the 
choice of nib — i-e., pole mass or running mass. 

0U5.EQ. 'N' ('Y') 

According to the strategy that all non-leading and gauge-variant quantities should 
be expanded, in Eq.( p05|) we use as the zero order approximation. 



47ra(M^) 



(309) 



However, since perturbation theory rearranges itself in such a way that the expansion 



of sin^ 6*^3 starts with s^, i.e.. 



sin^^.' 



off 



-6g^ +- {As' 
2 ^^ 2 



I 

A ' 



(310) 



we have envisaged the possibility of reorganizing the structure of the pseudo-obser- 
vables such that to all orders the bare weak mixing angle has an expansion starting 
with s^. Combined with OUl.EQ. 'N' ('Y') this option tells us that the expansion 
parameter, formally a/^Ans'^c^), can be set to 



G„Ml 



X 



{l;p«;l-Aa;p^(l-A«)} . 



(311) 



0U6.EQ. 'N' CY') 

This option deals with factorization of electroweak corrected kernels and QCD radi- 
ation. For instance, if 0U2 . EQ . ' N ' — the expanded solution — we still distinguish 
between a non-factorized solution, 



1-6 



f 

z / 



(312) 
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and a fully factorized solution. A special treatment is of course devoted to bb final 
state in order to reproduce the FTJR correction term. 

• 0U7.EQ. 'N' CY') 

Higher order QCD corrections to are implemented with the exact AFMT term 
or by subtracting from S^*^""™ the leading rrit term and by replacing it with the 
corresponding one evaluated at a scale of ^ = 0.248m(. 

To summarize TOPAZO runs in 2^ or 2^ electroweak x2 QCD options. While we have 
devoted a special section to the effect of DUO . EQ . ' N ' and 0U2 . EQ . ' Y ' , here we present a 
short Table |2^ to indicate the increasing spread in the theoretical errors while the various 
flags are switched on. 



5.4 ZFITTER options 

In this subsection we describe some electroweak and QCD options implemented in ZFITTER 
(version 4.9). Simultaneously, we give a description of the flags, implemented in version 
4.9 for studying the theoretical uncertainties. 



OZl: IAMT4=3,2,1 

The realization of leading and remainder terms, given by formulae ([fl|-|8^), is the 



default, IAMT4=3. If IAMT4=2, then X is not included in the leading terms, it re- 



mains a part of remainders |8j|. If IAMT4=1, then both X and 21oop-aQ;s terms are 
in the remainders. In the last case the Sloop-aa^ term is also placed in remainders. 
In the result of numerical investigations it was revealed that for IAMT4=2 the re- 
mainder terms in and kj are not small. Since this contradicts our philosophy to 
keep remainder terms small, one could exclude it from the set of working options. 
However, it was found that the difference between IAMT4=3 and 1 is rather small, 
therefore this option doesn't sizably influence the theoretical errors. 

0Z2: IHIGS=0,1 

This option governs the resummation of the leading Higgs contribution in Ar 



Ar«'^^^^^^-V^- In-f-- , M^:^M^. (313) 




If IHIGS=0 (the default), then the Higgs contribution is not re-summed. If IHIGS=1 
and if 

M2 5 , , 

In^f-gSO. (314) 

W 

then it is re-summed — i.e., it is extracted from remainders with the scale a/s^r and 
put to the leading terms with the scale G^, as in ( pi3| ). We observed, that 10/12 of 
^^Higgs jg contained in X. Therefore, if IAMT4=3, then only 1/12 of it is additionally 
re-summed. The influence of this option on theoretical errors was found to be tiny. 
For this reason, the Higgs resummation has not been implemented in and Hf. 



98 



• 0Z3: ISCRE=0,1,2 

This option defines the scale of the remainder terms. If ISCRE=0 (the default), 
then the scale is a/s^, if ISCRE=2 it is G^. (More precisely, in pf and Hf and 
G^{1 — Act) in Ar.) For ISCRE=1, the scale of the remainder in Ar is set equal 
to Gf^{l — Ar^), which was not included in the set of working options, since its 
influence is much smaller then the previous one's. The effect of the variation of the 
scale of the remainder term on the theoretical bands was found to be dominating. 
Especially influential is the scale variation in Ar, which introduces terms of the 
order c^/s^ Ap ■ Arj-cm- This illustrates the importance in the calculation of the 
next-to- leading term of the order 0{G'i 

• 0Z4: IFACR=0 , 1 , 2 , 3 

It realizes four subsequent expansions of Ar as they are given in (pQl) . The first, 
fully non-expanded option, is the default. Only the first expansion (the second row) 
was retained in the set of working options. The last two were excluded, since they 
contradict Sirlin's theorem on mass singularities. This leads to a visible spread of 
theoretical bands, although one much smaller than that of the previous option. 

• 0Z5: IFACT=0, 1,2,3,4,5 

It realizes six subsequent expansions for pf and k/. The default, IFACT=0, corre- 
sponds to the non-expanded realizations (f7^([73l). The four first options for p/, 
IFACT=0 ,1,2,3, are exactly the same as are given in (0), while is expanded for 
IFACT=1,2,3: 



(315) 



For IFACT=4,5, we linearize the full expression (p^. Introducing 

1 

Pl 



1 - Ap ' 

see (0) and (|8g), we have for IFACT=4 



W L ' 



24^2 



77 



+ Pr 



W^VL /il^<-™ V 



Realizing, that 



= 1 + ARf 



V,A 



(316) 
(317) 



(318) 
(319) 

for IFACT=5 we finally end up with a fully expanded equation for the partial widths: 



N[ |p. \{9U?K + R{] + Ap;,_ Mj' + 1 



24^2 77 



(320) 
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It was found that the spread of error bands gradually grows while coming from 
IFACT=0 to IFACT=5. For this reason, only these two limiting cases were left among 
the working options. This option was found to be rather influential. 

0Z6: ISCAL=0, 1,2,3,4 

is the only QCD option. At the default, ISCAL=0, we implemented the exact AFMT 
correction. For ISCAL=2,1,3 we implemented the ^-factor as given in Ref. p5[ . 



Finally, for ISCAL=4, Sirlin's scale ^ = 0.248 (see Ref.[^) was implemented. Only 
ISCAL=0,4 were left among working options. 

To summarize, ZFITTER runs in 2^ electroweak x 2 QCD options Q Table ^ has been 
added to show the effect of the working options of ZFITTER on theoretical errors. 
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Figure 11: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for , including an 
estimate of the theoretical error as a function of m^, for = 300 GeV and — 0.125. 
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Figure 12: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for Te, including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 13: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for T^, including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 14: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for Ri, including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 15: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for Rb, including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 16: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for R^, including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 17: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for sin^ 6*^5, including 
an estimate of the theoretical error as a function of mj, for = 300 GeV and = 0.125. 
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Figure 18: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for sin^ 6*^5, including 
an estimate of the theoretical error as a function of mj, for = 300 GeV and = 0.125. 
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Figure 19: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for A^^^ , including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 



114 



_ll 

< 



0.16 



□ BHM 

LEPTOP 
TOPAZO 
A WOH 

* ZFITTER 







0.14 - 



0.12 



100 



150 



as(Mz) = 0.125 
Mh = 300 GeV 



200 



250 

rtit (GeV) 



Figure 20: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for A|°p*, including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 21: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for Aj^^^ , including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 22: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for A^^^ , including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 23: The BHM, LEPTOP, TOPAZO, ZFITTER, WOH predictions for o"^, including an 
estimate of the theoretical error as a function of rrit, for = 300 GeV and dg = 0.125. 
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Figure 24: The BHM (square), TOPAZO (diamond) and ZFITTER (cross) predictions, in- 
cluding an estimate of the theoretical error, for o"^ in a fully extrapolated set-up. Here 
rrit = 175 GeV, = 300 GeV and = 0.125. In the lower part a comparison is also 
shown with the relative deviation of BHM, ZFITTER versus TOPAZO. 
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Figure 25: The same as in Fig. 24 with an s' cut, s' = 0.5s. 
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Figure 26: The same as in Fig. 24 for the hadronic cross-section, cr'*. 
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Figure 27: The same as in Fig. 26 with an s' cut, s' = 0.01s. 
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Figure 28: The same as in Fig. 24 for the leptonic forward-backward asymmetry. In 
the lower part a comparison is also shown with the absolute deviation of BHM , ZFITTER 
versus TOPAZO. 
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Figure 29: The same as in Fig. 28 with an s' cut, s' = 0.5s. 
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Figure 30: The TOPAZO (diamond) and ZFITTER (cross) predictions , including an estimate 
of the theoretical error, for a'^ in the following set-up: 40° < 6- < 140°, ^acoU < 10° and 
Eth = 20GeV. Here rut = 175 GeV, iVf^ = 300 GeV and = 0.125. In the lower part a 
comparison is also shown with the relative deviation of ZFITTER versus TOPAZO. 
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Figure 31: The same as in Fig. 30 for ^acou < 25°. 



126 



88 



90 



92 



94 



v's(GeV) 



o 1.001 

«^ 

b 
b 

1.000 



0.999 



88 



90 



92 



94 



Vs(GeV) 



Figure 32: The TOPAZO (diamond) and ZFITTER (cross) predictions, including an estimate 
of the theoretical error, for cr^ with s-channel electrons, in the following set-up: 40° < 
9^ < 140°, e^oii < 10° and Et^ = IGeV. Here rut = 175 GeV, = 300 GeV and 
dg — 0.125. In the lower part a comparison is also shown with the relative deviation of 
ZFITTER versus TOPAZO. 
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Figure 33: The same as in Fig. 32 for ^acou < 25°. 
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Figure 34: The TOPAZO (diamond) and ZFITTER (cross) predictions, including an estimate 
of the theoretical error, for Ai^^ in the following set-up: 40° < 9- < 140°, 6*30011 < 10° and 
Eti, = 20GeV. Here rut = 175 GeV, = 300 GeV and = 0.125. In the lower part a 
comparison is also shown with the absolute deviation of ZFITTER versus TOPAZO. 
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Figure 35: The same as in Fig. 34 for ^acou < 25°. 
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Figure 36: The TOPAZO (diamond) and ZFITTER (cross) predictions, including an estimate 
of the theoretical error, for A^^^ with s-channel electrons in the following set-up: 40° < 
^_ < 140°, ^acou < 10° and E^], = 1 GeV.Here m* = 175 GeV, = 300 GeV and 
dg — 0.125. In the lower part a comparison is also shown with the absolute deviation of 
ZFITTER versus TOPAZO. 
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5ure 37: The same as in Fig. 36 for ^acou < 25°. 
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Table 1: Maximum Derivatives with respect to a 
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r„(MeV) 


1.3561 


1.3656 


1.3664 


1.3613 


fMeV) 

■■- a ^^■'■v-'-c V j 
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-0.91531x10-^ 


-0.89791x10-^ 


-0.92040x10-^ 




0.14434x10"^ 


0.14411x10"^ 


0.14411x10"! 


0.14421x10"! 


Th (MeV) 


6.3768 


6.4335 


6.4408 


6.4143 


0.16748xl0"2 


0.16718xl0"2 


0.16722xl0"2 


0.16738xl0"2 


Tinv (MeV) 


-0.19877x10-^ 


-0.37484x10-2 


-0.38416x10-2 


-0.56691x10-2 


FB 


0.11116x10-^ 


0.11095x10-1 


0.11098x10-1 


0.11098x10-1 




0.14772x10-3 


0.14744x10-=^ 


0.14730x10-3 


0.14698x10-3 
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Table 2: Maximum Derivatives with respect to mf,. 



Observables 


BHM 


LEPTOP 


TOPAZO 


ZFITTER 


Tb (MeV) 
(MeV) 
Ri 

< (nb) 

Rb 

r, (MeV) 
Rc 


-0.77934 

-0.81949 

-0.92763x10^2 

0.73002x10-2 
-0.34685x10-''^ 

-0.80799 
0.76013x10-^ 


-0.81333 
-0.81344 

-0.96913x10-2 
0.76691x10-2 

-0.36581x10-3 

-0.81344 
0.80403x10-^ 


-0.80247 
-0.84278 
-0.95275x10-2 
0.74958x10-2 
-0.35739x10-3 

-0.83069 
0.78998x10-^ 


-0.79613 
-0.79462 

-0.94665x10-2 
0.74927x10-2 

-0.35819x10-3 

-0.79462 
0.79395x10-^ 



134 



Tabic 3: The experimental data. 



Observables 


Data 




M^(GeV) 


80.22 ±0.18 


(MeV) 


83.96 ±0.18 


(MeV) 


2497.4 ± 3.8 


(nb) 


41.49 ±0.12 


Ri 


20.795 ± 0.040 


Rb 


0.2202 ± 0.0020 


Rc 


0.1583 ±0.0098 




0.2321 ± 0.0004 




0.0170 ±0.0016 


A" 

FB 


0.0967 ± 0.0038 




0.0760 ± 0.0091 


Ah(SLD) 


0.1668 ±0.0077 
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Table 4: The experimental data for M^,r; = Fe,!^ and the theoretical predictions 
corresponding to m* = 175 GeV, = 300 GeV and as(M^) = 0.125 The first entry is 
BHM then LEPTOP, TQPAZO, WOH and ZFITTER. The uncertainties quoted are obtained from 
a variation of program options as described in Section 5. 



Observable 


Exp. 


Theor. Predictions 


Average 


(GeV) 


80.22 ±0.18 


80.3191°:°°? 
80.312l°:°li 

80.3101°:°°? 
80.3191°:°°^ 
80.3171°:°°? 


80.315 


(MeV) 


83.96 ±0.18 


83.9191°:°?° 
83.9301°:°! 
83.93ll°:°i^ 

^ ^ ^ . ^ _i_n noo 

83.943lHi 
83.9411°:°^? 


83.933 


(MeV) 


2497.4 ± 3.8 


2497.41?:^ 
2497.2ll:i 
2497.41^:^ 
2497.4lJ:i 
2497.5lg:^ 


2497.4 
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Table 5: The same as in Table ^ for Ri, Rt, Rc- 



Observable 


Exp. 


Theor. Predictions 


Average 


Ri 


20.795 ± 0.040 


20.788t°:°°^ 
20.7801°:°°^ 
20.7821°:°°^ 
20.780l°:°J;j 
20.78lt°:°°? 


20.782 


Rb 


0.2202 ±0.0020 


n 01 c 77+0.00010 
V.ZLOl < -0.00011 

91 c;fi4+o-oooo9 

U.Zi004_Q 00004 

91 c;r7+0-00003 

U-^-LOO' -0.00012 

21 1^67+0-00018 

U.ZidO^_o.oo006 

n 91 c;7l +0.00001 
U.ZiO 1 i-0.00002 


0.21569 


Rc 


0.1583 ±0.0098 


n 1 7OQC+0.00002 
u.i ^ ZOD_0.00002 

1 7940+0-00002 

u.i ^Z4U_o.00003 

n 1 7OQ7+0.00004 
U'-L'^O 1-0.00000 

1 7240+0-00001 

u.i ^Z4U_o.oo003 

n 1 7OOC+0. 00002 
u.i <^00_o.00000 


0.17238 
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Table 6: The same as in Table ^ for sin^ 6'^^, sin^ 6'^^, 



Observable 


Exp. 


Theor. Predictions 


Average 


sin^ eis 


0.2321 ±0.0004 


n Q7+0. 00004 
U-^-jiy '-0.00007 

n oQonn+o.oooos 

U.ZOZUU_0.00008 

n oo9nn+o.oooo4 

U.ZOZUU_0.00004 
n 9'?1 Q4+0-00003 

n 9Q9nc;+o.oooo4 

U.ZOZUO_o.00014 


0.23199 


sin^ e',s 




n 01.1.'3.^ +0.00004 

U.ZOOOi_o.oo012 
n 9QQ9Q+0.00008 

u.zoozy_o,oooio 
n 9QQQn+o.oooo9 

U.ZOOOU_o. 00001 

n 9QQ9c:+0.0004 
U.ZOOZO_0,00007 

n 9QQQc;+0. 00004 
U.ZOOOO_o.00014 


0.23330 


A' 

FB 


0.0170 ±0.0016 


0.01544l°:S 
0.015391°:™ 

01 c^-ifi+O-OOOOS 
U.UidOO_o.oo007 

n ni c;4Q+o.oooi2 

U.UiO4y_0.00005 

01 ciQI +0.00024 
U.UiO-3i_o.00007 


0.01540 
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Table 7: The same as in Table H for A,„,A'' ,A'^ . 

I_] LR' fb' FB 



Observable 


Exp. 


Theor. Predictions 


Average 




0.1668 ±0.0077 


n 1 4Q4fi+0-00052 
U.i4J40_o.00031 

U.i^OZD_o,oo060 

n 1 /1097+0.00028 
U.i4JZ^_0.oo031 

n 1 /IQ79+0.00057 
U.i'iO 1 ^-0.00024 

n 1 ZL98Q+0-00110 

u.i4/:8y_o.ooo32 


0.14332 


A" 

FB 


0.0967 ±0.0038 


U.iUUdO_o.oo022 
U.iUU^U_o.oo042 
U.iUUOO_o.00023 
U.iUU 1 ^-0.00006 

0.10013i°;S 


0.10042 


A^ 

FB 


0.0760 ±0.0038 


0.07169t°;™ 
0.07158i°:S 
0.071591™ 

n 071 OQ+0.00032 

u.ua«d_o.oooi3 

n nvT Qo+o. 00061 
u.uad»_o.oooi7 


0.07161 
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Table 8: Largest half-differences among central values (dc) and among maximal and 
minimal predictions (dg) for 150 GeV < rrit < 200 GeV, 60GeV < < 1 TeV and 
0.118 < a^iM^) < 0.125. 



Observable 


dcO 


dgO 


M„ (GeV) 


6.5 X 10- 


-3 


1.9 X 10- 


-2 


Te (MeV) 


1.7 X 10- 


-2 


3.7 X 10- 


-2 


r. (MeV) 
sin' 


0.3 




1.6 




6.5 X 10- 


-5 


1.5 X 10" 


-4 


sin' ei^ 


6.0 X 10- 


-5 


1.6 X 10- 


-4 


Ri 


4.0 X 10- 


-3 


1.0 X 10- 


-2 


Rb 


7.0 X 10- 


-5 


2.0 X 10- 


-4 


Rc 


3.0 X 10" 


-5 


5.0 X 10- 


-5 


a', (nb) 


7.5 X 10" 


-3 


8.5 X 10- 


-3 




1.2 X 10- 


-4 


2.5 X 10- 


-4 


A" 

FB 


3.5 X 10" 


-4 


8.2 X 10- 


-4 


A" 

FB 


2.7 X 10- 


-4 


6.3 X 10- 


-4 




5.0 X 10- 


-4 


1.1 X 10- 


-3 
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Table 9: Largest half-differences among central values {dc) and among maximal and 
minimal predictions (dg) for rrit = 175 GeV, 60 GeV < Af ^ < 1 TeV and ^^(M^) = 0.125. 



Observable 


dcO 




(GeV) 


4.5 X 10- 


-3 


1.6 X 10- 


-2 


Te (MeV) 


1.3 X 10" 


-2 


3.1 X 10- 


-2 


(MeV) 
sin' eif, 


0.2 




1.4 




5.5 X 10" 


-5 


1.4 X 10" 


-4 


sin' 


5.0 X 10" 


-5 


1.5 X 10- 


-4 


Ri 


4.0 X 10" 


-3 


9.0 X 10- 


-3 


Rb 


6.5 X 10- 


-5 


1.7 X 10- 


-4 


Rc 


2.0 X 10" 


-5 


4.5 X 10- 


-5 


cTo^ (nb) 


7.0 X 10" 


-3 


8.5 X 10" 


-3 


FB 


9.3 X 10" 


-5 


2.2 X 10" 


-4 


A" 

FB 


3.0 X 10" 


-4 


7.4 X 10" 


-4 


^FB 


2.3 X 10- 


-4 


5.7 X 10- 


-4 




4.2 X 10- 


-4 


8.7 X 10- 


-4 
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Table 10: Largest half-differences among central values (dc) and among maximal and 
minimal predictions (dg) for rrit = 175 GeV, = 300 GeV and 0.118 < d^iM^) < 0.125. 



Observable O 


dc 


O 






dg 







(GeV) 


4.5 X 


10- 


-3 


1.3 


X 


10- 


-2 


Te (MeV) 


1.2 X 


10- 


-2 


3.0 


X 


10- 


-2 


(MeV) 
sin' 


0. 


L5 






1. 


4 




5.5 X 


10- 


-5 


1.1 


X 


10- 


-4 


sin' el^ 


5.0 X 


10" 


-5 


1.1 


X 


10- 


-4 


Ri 


4.0 X 


10- 


-3 


9.0 


X 


10- 


-3 


Rb 


6.5 X 


10- 


-5 


1.6 


X 


10- 


-4 


Rc 


2.5 X 


10- 


-5 


4.5 


X 


10- 


-5 


< (nb) 


6.5 X 


10- 


-3 


8.0 


X 


10- 


-3 




8.9 X 


10- 


-5 


1.9 


X 


10- 


-4 




3.0 X 


10- 


-4 


6.1 


X 


10- 


~4 




2.3 X 


10- 


-4 


4.7 


X 


10- 


-4 




4.2 X 


10- 


-4 


8.6 


X 


10- 


~4 



142 



Table 11: Effects of additional options in TOPAZO for = 175 GeV, 
cf^(M^) = 0.125. 



300 GeV and 



Observables 



M^(GeV) 
(MeV) 
r. (MeV) 
Ri 

Rb 

Rr 



sin^ 



Ai 



FB 



Predictions 



80.3101^°° 



007 

+0.015 



80.310ir°' 



83.931^1j:lji^ 



007 

+0.043 



2497.4lg: 



20.7821°:°°^ 



n 91 c;fi7+0-00003 

U.ZiOU i -0.00012 

A 1 7007+0.00004 
U.l/-iO/_o.oOOOO 

U.Z>3ZUU_o.o0004 



n 01 c;3fi+0-00008 
U.UiOOD_o Q0007 



1 0033+°-°°°23 

U.1UUOO_0.00023 



0.07159l°;ri? 

n 1 4007+0.00028 
U.l^OZ / -0.00031 



83.931^;^:^^^^ 

2497.41JI 



20.7821HJ° 



0.21569l°:S 



A 1 7007+0.00007 
U.l /ZO/ -o.ooOOl 



n o-^oni +0.00004 

U.ZOZU1_0.00026 



n ni c: 00 +0.00075 
U.U1005_o.oo074 



n 1 nn'^7+000137 

U.lUU>3/_o.oo067 



n 071 c;7+0-00ii5 

U.U i iO i -0.00053 



1 4QOn+000193 
U.i4C)ZU_o.oo039 
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Realistic-observables 



Tabic 12: The hadronic cross-section (nb) in two different configurations, NN=NY/ 
YN=YY for inclusion of initial state p'MV ]:)roduction. 



^ (GeV) 


BHM 


TOPAZO 


ZFITTER 


88.45 


5.181 


5.184 


5.185 




5.168 


5.171 


5.172 


89.45 


10.062 


10.068 


10.067 




10.036 


10.042 


10.042 


90.20 


18.033 


18.040 


18.039 




17.983 


17.992 


17.992 


91.1887 


30.446 


30.452 


30.451 




30.366 


30.375 


30.373 


91.30 


30.585 


30.590 


30.590 




30.506 


30.514 


30.514 


91.95 


25.176 


25.176 


25.181 




25.124 


25.126 


25.130 


93.00 


14.131 


14.127 


14.134 




14.117 


14.119 


14.124 


93.70 


10.070 


10.065 


10.072 




10.068 


10.073 


10.074 
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Table 13: The fi forward-backward asymmetry in four different configurations, NN/ 
YN/NY/YY for inclusion of initial state pair production and initial-final QED inter- 
ference. 





BHM 


TOPAZO 


ZFITTER 


88.45 


-0 2552 


-0.2551 


-0.2546 




-0 2552 


-0.2557 


-0.2544 




-0.2546 


-0.2543 


-0.2541 




-0.2546 


-0.2549 


-0 2539 


89.45 


-0.1631 


-0.1632 


-0.1630 




-0.1631 


-0.1636 


-0.1628 




-0.1626 


-0.1624 


-0.1625 




-0.1626 


-0.1628 


-0.1623 


90.20 


-0.0912 


-0.0914 


-0.0913 




-0 091 2 


-0 091 7 


-0 0912 




-0.0907 


-0.0907 


-0 0908 




-0.0907 


-0 0909 


-0.0907 


91.1887 


-0.0012 


-0.0015 


-0.0014 




-0.0012 


-0.0014 


-0.0014 




-0 0008 


-0 0011 

\J .\J\J J. J. 


-0 0011 

\J •\J\J X X 




-0 0008 


-0.0011 


-0.0011 


91.30 


0080 


0.0076 


0.0077 




0.0080 


0.0076 


0.0077 




0.0083 


0.0078 


0.0080 




0.0083 


0.0078 


0.0083 


91.95 


0.0556 


0.0550 


0.0552 




0.0556 


0.0551 


0.0551 




0.0558 


0.0549 


0.0554 




0.0558 


0.0550 


0.0553 


93.00 


0.1114 


0.1105 


0.1108 




0.1114 


0.1106 


0.1106 




0.1116 


0.1104 


0.1109 




0.1116 


0.1104 


0.1106 


93.70 


0.1379 


0.1368 


0.1372 




0.1379 


0.1367 


0.1367 




0.1380 


0.1367 


0.1372 




0.1380 


0.1366 


0.1368 
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Table 14: De-convoluted A^^^ at s = M^, first entry is the ZZ part, then ZZ + 77, 
ZZ + Z7 and total. 



Contribution 


BHM 


TOPAZO 


ZFITTER 


ZZ 
ZZ + 77 
ZZ + Z7 
Total 


0.015443 
0.015351 
0.017026 
0.016925 


0.015358 
0.015267 
0.016735 
0.016636 


0.015279 
0.015189 
0.016725 
0.016627 
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Table 15: The de-convoluted a", and 



y\je\ ) 






7T7TTTT7D 
Zil" 111 tn 


OO.40 


U.o4y i 


n Q/1 oo 


n Q/ioo 
U.o4yz 






( .Uzo 


7 noo 




— U.zooD 


— U.zooo 


— U.ZOo4 


on A K 

oy.4o 


U.Doo4 


n i^QQ7 


A CQQ7 




14. UM 


1 /I HQQ 


1 /I n77 
14. / 




— U. 1404 


n 1 /I 
— U.14D0 


n 1 /I c;q 
— U. 140o 






1 9457 


1 9458 

J- . Z/THJO 




25.650 


25.663 


25.655 




-0.0749 


-0.0751 


-0.0750 


91.1887 


2.0015 


2.0019 


2.0022 




41.400 


41.409 


41.402 




0.0169 


0.0166 


0.0166 


91.30 


1.9812 


1.9816 


1.9820 




40.983 


40.991 


40.987 




0.0271 


0.0268 


0.0268 


91.95 


1.4486 


1.4488 


1.4492 




29.935 


29.936 


29.937 




0.0855 


0.0853 


0.0851 


93.00 


0.6534 


0.6533 


0.6536 




13.411 


13.408 


13.411 




0.1757 


0.1755 


0.1751 


93.70 


0.4105 


0.4104 


0.4106 




8.362 


8.358 


8.361 




0.2323 


0.2320 


0.2316 
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Table 16: 5conv, as defined by Eq.( 



136) for the hadronic cross-section. 



v^(GeV) 


BHM 


TOPAZO 


ZFITTER 


88.450 


-0.2638 


-0.2640 


-0.2634 


89.450 


-0.2869 


-0.2869 


-0.2866 


90.200 


-0.2989 


-0.2989 


-0.2987 


91.189 


-0.2665 


-0.2665 


-0.2664 


91.300 


-0.2556 


-0.2556 


-0.2555 


91.950 


-0.1607 


-0.1607 


-0.1606 


93.000 


+0.0526 


+0.0530 


+0.0532 


93.700 


+0.2040 


+0.2052 


+0.2049 
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Bhabha scattering 



Table 17: The ALIBABA(A) - TOPAZO(T) comparison for the full Bhabha cross-section (in 
pb) for the following set-up: 40° < 9- < 140°, ^acou < 10° and Eth = 1 GeV. I is the 
TOPAZO default for QED final state radiation, II is the ZFITTER-hke default for QED final 
state radiation and III includes initial state pair production. 



V^GoV 


A 


T(I) 


T(II) 


T (III) 


88.45 
89.45 
90.20 
91.19 
91.30 
91.95 
93.00 
93.70 


457.52 ±0.27 
643.95 ±0.31 
908.99 ±0.39 
1183.99 ±0.39 
1163.56 ±0.45 
876.90 ±0.28 
481.35 ±0.14 
355.57 ±0.13 


457.30tg:o6 ± 0.25 
644.37l|i:^^ ± 0.24 
910.461^:?? ± 0.24 
1184.03l[i:l| ±0.24 
1163.5ltg:f|±0.24 
874.021^:?^ ± 0.24 
m.blt'^o^l ± 0.24 
352.52lg:J^ ± 0.25 


457.20 

644.23 

910.26 

1183.78 

1163.26 

873.83 

477.41 

352.45 


456.14 

642.60 

907.86 

1180.79 

1160.38 

872.13 

477.18 

352.72 
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Table 18: The same as in Table 



|T7| 



for the forward backward asymmetry. 



\/ r-J * — ' Vj V 


A 


T(I) 




88.45 


0.44611 ± 1.06 X 


10- 


~3 


O.UbSAtlfs ± 0.79 X 


10- 


-3 


0.44637 


8Q 45 


0.34250 ± 0.90 x 


10- 


-3 


0.34166t[j:|]^ ± 0.51 X 


10- 


-3 


34252 


■>J\J . ^j\J 


0.24956 ±0.81 x 


10" 


-3 


0.24977^{]:|j^ ± 0.33 x 


10" 


-3 


25043 


91.19 


0.13925 ±0.66 x 


10- 


-3 


0.13916i[;:[J^ ± 0.23 X 


10- 


-3 


0.13951 


91.30 


0.13050 ±0.70 X 


10- 


-3 


0.13035t|];|i| ± 0.23 x 


10- 


-3 


0.13067 


91.95 


0.10169 ±0.63 X 


10- 


-3 


0.10139t[]:|]^±0.30 X 


10- 


-3 


0.10158 


93.00 


0.13110 ±0.61 X 


10- 


-3 


0.13055iE]:?^ ± 0.57 X 


10- 


-3 


0.13061 


93.70 


0.18157 ±0.68 X 


10- 


-3 


0.17957;{];?^ ± 0.83 x 


10" 


-3 


0.17944 



150 



Table 19: The same as in Table for ^acou < 25°. 



v/iGeV 


A 


T(l) 


T(ll) 


T (111) 


88.45 


483.27 ±0.25 


A84:.97tom ± 0.22 


484.86 


483.74 


89.45 


672.67 ±0.29 


673.99tl:ll ± 0.22 


673.84 


672.13 


90.20 


942.52 ±0.34 


942.96i[j:?^ ± 0.22 


942.76 


940.27 


91.19 


1218.66 ±0.40 


1219.24l[]:?^ ±0.21 


1218.98 


1215.90 


91.30 


1198.23 ±0.36 


1198.42l[]:?^ ±0.21 


1198.17 


1195.20 


91.95 


908.87 ±0.30 


905.35l[J;?| ± 0.20 


905.16 


903.39 


93.00 


505.38 ±0.15 


504:.2Atom ± 0.20 


504.13 


503.88 


93.70 


378.20 ±0.13 


377Mtom ± 0.20 


377.86 


378.14 
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Table 20: The same as in Table |T8| for ^acou < 25°. 

T(I) 



^/sGeV 



A 



T(lll) 



88.45 
89.45 
90.20 
91.19 
91.30 
91.95 
93.00 
93.70 



0.45843 ±0.94 x 10 
0.35479 ±0.77 x 10 
0.26121 ±0.67 X 10 
0.15114 ±0.70 X 10 
0.14067 ±0.62 X 10 
0.11466 ±0.62 X 10 
0.15628 ±0.59 x 10 
0.21239 ±0.66 x 10 



0.4606lt[j:?^ ± 0.65 X 10 
0.35560tJ]:?7 ± 0.43 x 10 



-3 



± 0.29 X 10 



0.26165tE]:?l 

0.15045i[j:?^ ± 0.20 X 10 
0.14203tl]:?^±0.20 X 10 
0.11773tl]:??±0.25 X 10" 
0.15838i[]:^^±0.45 x 10" 



3 
3 
3 

-3 
"3 
3 







.21360t{]:^^ ± 0.65 X 10-3 



0.46168 
0.67213 
0.26235 
0.15083 
0.14238 
0.11795 
0.15845 
0.21344 



Table 21: The same as in Table [I^ with the percentage relative deviation. 



^GeV 


A 


T(II) 


% 


88.45 


457.52 


457.20 


±0.07 


89.45 


643.95 


644.23 


-0.04 


90.20 


908.99 


910.26 


-0.14 


91.19 


1183.99 


1183.78 


±0.02 


91.30 


1163.56 


1163.26 


±0.03 


91.95 


876.90 


873.83 


±0.35 


93.00 


481.35 


477.41 


±0.82 


93.70 


355.57 


352.45 


±0.88 
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Table 22: The same as in Table ^ for s-channel alone. 



^/sGeV 


A 


T(II) 


% 


o o /I r 

88.45 


173.01 


-1 ^70 '71 

172.71 


+0.17 


89.45 


330.62 


330.70 


-0.02 


90.20 


588.47 


588.81 


-0.06 


91.19 


989.81 


990.90 


-0.11 


91.30 


994.07 


995.38 


-0.13 


91.95 


819.56 


819.97 


-0.05 


93.00 


463.72 


461.68 


+0.44 


93.70 


331.78 


329.83 


+0.59 



Table 23: ALIBABA-TOPAZO comparison for s — t and t — t contributions with the following 
set-up: 40° < ^_ < 140°, ^acoU < 10° and Eth = 1 GeV. 



v/iGeV 


a - cr(s) A 


cr — cr(s) T 


<5(A) 


5(T) 


% 


88.45 


284.51 


284.49 


1.644 


1.647 


-0.18 


89.45 


313.33 


313.53 


0.948 


0.948 


+0.00 


90.20 


320.52 


321.45 


0.545 


0.546 


-0.18 


91.19 


194.18 


192.88 


0.196 


0.195 


+0.51 


91.30 


169.49 


167.88 


0.171 


0.169 


+ 1.18 


91.95 


57.34 


53.86 


0.070 


0.066 


+5.88 


93.00 


17.63 


15.75 


0.038 


0.034 


+11.11 


93.70 


23.79 


22.62 


0.072 


0.069 


+4.26 
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Table 24: The same as in Table 



for the forward-backward asjTiimetry. 



v/iGeV 


A s+t 


T s+t 


A s 


T s 


A t 


T t 


88.45 


0.44611 


0.44534 


-0.22019 


-0.22060 


0.66630 


0.66594 


89.45 


0.34250 


0.34166 


-0.13754 


-0.13847 


0.48004 


0.48013 


90.20 


0.24956 


0.24977 


-0.07774 


-0.07649 


0.32730 


0.32626 


91.19 


0.13925 


0.13916 


-0.00102 


0.00001 


0.14027 


0.13915 


91.30 


0.13050 


0.13035 


0.00745 


0.00779 


0.12305 


0.12256 


91.95 


0.10169 


0.10139 


0.05059 


0.04851 


0.05110 


0.05288 


93.00 


0.13110 


0.13055 


0.09863 


0.09789 


0.03247 


0.03266 


93.70 


0.18157 


0.17957 


0.12237 


0.12248 


0.05920 


0.05709 
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Table 25: Comparison for the full Bhabha cross-section (in pb) for the following set-up: 
40° < 6- < 140°, 6'acoii < 10° and Eth = 1 GeV. First entry is ALIBABA, second entry is 
the TOPAZO default for QED final state radiation and no pair production, 5fsr is the 
uncertainty on final-state QED radiation estimated by TOPAZO, while Ag is the relative 
difference between the maximal and minimal predictions of ALIBABA-TOPAZO. 



^/sGeV 

V/ o vjrc V 


A 


TCT^l 




A % 




4^7 CO _|_ n 97 


A^i '^n+°-^^ -1- n 9'i 


09 


n 1 

yj.i-O 


89.45 


643.95 ±0.31 


644.371^:^^ ± 0.24 


0.02 


0.19 


90.20 


908.99 ±0.39 


910.46tg:?? ± 0.24 


0.02 


0.27 


91.19 


1183.99 ±0.39 


1184.03lg:t| ± 0.24 


0.02 


0.10 


91.30 


1163.56 ±0.45 


1163.5ll[J;g ± 0.24 


0.02 


0.10 


91.95 


876.90 ± 0.28 


874.02^^:?^ ± 0.24 


0.02 


0.37 


93.00 


481.35 ±0.14 


477.511^:^^ ± 0.24 


0.02 


0.86 


93.70 


355.57 ±0.13 


352.52ig:J^ ± 0.25 


0.02 


0.95 
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Effect of working options for different codes 



Table 26: The effect of the working options of BHM on theoretical errors. 



Obsor\'al)lo 


Default 


IRES err 


+IQCD on- 


+IFAC err 


+rTA\'(:) err 


(GeV) 
Te (MeV) 
(MeV) 

Ri 
Rb 

A' 

FB 
FB 

sin' els 


80.319 
83.919 
2497.4 
20.788 
0.21577 
0.015435 
0.10053 
0.23197 
0.23331 


0.004 
0.020 
0.68 
0.001 
0.00000 
0.00011 
0.00037 
0.00007 
0.00007 


0.010 
0.032 
1.07 
0.002 
0.00000 
0.00018 
0.00059 
0.00011 
0.00011 


0.010 
0.033 
1.71 
0.010 
0.00011 
0.00018 
0.00059 
0.00011 
0.00011 


0.010 
0.033 
1.92 
0.012 
0.00020 
0.00018 
0.00059 
0.00011 
0.00016 



Table 27: The effect of the working options of LEPTOP on theoretical errors. The first 
two lines indicate parametric uncertainties caused by Ss' = 0.0003 (which is equivalent 
to Sa~^ = 0.12) and by 6mi, = 0.3GeV. The next six lines refer to the intrinsic theoretical 
uncertainties. This table was calculated with mt=175 GeV, M^=300 GeV and ds=(}.125. 









sin^ ^^ff 


<{nh) 






Ri 


r6 


RbX 






(MeV) 


(MeV) 






(MeV) 


(MeV) 




(MeV) 


10^ 


10^ 




16 


.015 


.00031 


.0014 


.8 


.8 


.0055 


.15 


1.1 


1.8 










.0023 


.2 


.2 


.0029 


.24 


11.0 


2.4 




9 


.015 


.00005 


.0002 


.5 


.4 


.0009 


.08 


.2 


.3 




5 


.008 


.00003 


.0001 


.3 


.2 


.0005 


.04 


.1 


.2 










.0029 


.3 


.3 


.0037 




3.8 


1.5 










.0010 


.1 


.1 


.0012 


.10 


4.6 


1.0 


total 


+13 


+.023 


+.00008 


+.0032 


+ 1.2 


+ 1.0 


+.0063 


+.23 


+8.7 


+ 1.9 




-13 


-.023 


-.00008 


-.0042 


-1.1 


-0.9 


-.0051 


-.13 


-4.1 


-2.9 
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Table 28: The effect of the working options of TOPAZO on theoretical errors. 



Observable 


Default 


OUl err 


+ 0U4 err 


+ 0U5 err 


+ 0U6 err 


+0U7 err 


M,, (GeV) 


80.310 






0.001 


0.001 


0.007 


Te (McV) 


83.931 


0.014 


0.014 


0.015 


0.015 


0.027 


(MeV) 


2497.4 


0.10 


0.10 


0.10 


0.30 


0.70 


Ri 


20.782 


0.004 


0.004 


0.004 


0.006 


0.007 


Rb 


0.21567 


0.00012 


0.00012 


0.00013 


0.00014 


0.00015 


A' 

FB 


0.015362 


0.00006 


0.00006 


0.00008 


0.00008 


0.00015 


FB 


0.10033 


0.00018 


0.00019 


0.00024 


0.00024 


0.00046 


sin' 


0.23200 


0.00004 


0.00004 


0.00004 


0.00004 


0.00008 


sin^ 


0.23330 


0.00005 


0.00005 


0.00006 


0.00006 


0.00010 



Table 29: The effect of the working options of ZFITTER on theoretical errors. 



Observable 


Default 


OZl err 


+0Z2 err 


+0Z3 err 


+0Z4 err 


+0Z5 err 


+0Z6 err 


(GeV) 


80.317 


0.001 


0.002 


0.005 


0.007 


0.007 


0.014 


Te (MeV) 


83.941 


0.004 


0.006 


0.014 


0.014 


0.021 


0.034 


r. (MeV) 


2497.5 


0.1 


0.2 


0.2 


0.2 


0.6 


1.1 


Ri 


20.781 




0.001 


0.002 


0.003 


0.005 


0.007 


Rb 


0.21571 






0.00001 


0.00001 


0.00002 


0.00003 




0.01531 


0.00003 


0.00009 


0.00015 


0.00024 


0.00024 


0.00030 




0.10013 


0.00011 


0.00029 


0.00050 


0.00078 


0.00079 


0.00101 


sin' 


0.23205 


0.00002 


0.00005 


0.00009 


0.00014 


0.00014 


0.00018 


sin' 91^ 


0.23335 


0.00002 


0.00005 


0.00009 


0.00014 


0.00014 


0.00018 
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Note added in proof 



While making proof-reading of this contribution we have been informed on some recent 
development concerning the AFMT term in Ap. A recent calculation by K.G. Chetyrkin, 



J.H. Kiihn and M. Steinhauser ||6^ as well as a revised version of the AFMT calcula- 
tion [53l have shown that the correct coefficient of ({4) in a term proportional to Cp 
of 6^-^ is 4 and not 188/5. This will correspond to some shift in our predictions for 
the central values of the pseudo-observables, as a matter of fact a shift common to all 
codes since the AFMT term is a common external block. We have shown this shift in 
the following Table, where we compare our error bands with the shifted central values at 
the standard reference point. The result can be simply summarized by saying that the 
updated central values remain within our theoretical error bands. As for the theoretical 
bands is concerned, we note first, that they are not dominated by the QCD- uncertainty 
in the calculation of the p-parameter. Second, as it seen from the Table added, the widths 
of the uncertainty bands is only marginally affected by the shift of the AFMT correction. 



Table 30: Change of some of the observables due to the introduction of the revised AFMT 
formulas. Here = 175 GeV, = 300 GeV and = 0.125. In order to estimate the 
size of the non-leading QCD effects, the 6'^'^° correction factor has been implemented 
according to the formulation of Ref . , with a scale which gives the maximum variation 
with respect to the AFMT(revised) term — ^ = 0.248(0.204) — and the difference between 
this and the AFMT (revised) calculation is used as an estimate of the corresponding 
uncertainty. The numbers are calculated by TOPAZO, first row, and ZFITTER, second row. 



Observable 


New AFMT 


AFMT 


Central difference 


(GeV) 


80.307^1{:r 
80.314tHS? 


80.317tnS^ 


-3 MeV 
-3 MeV 


Ti (MeV) 


83.926in|^ 
83.936i°:°^? 


83.93liHI^ 
83.94ltn^? 


-0.005 MeV 
-0.005 MeV 


(MeV) 


2497.2^^:^ 
2497.3+|ji 


2497.4^^:^ 
2497.5+|]i 


-0.2 MeV 
-0.2 MeV 


Ri 


20.782i[!:^^^ 
20.780i°:°°^ 


20.782iH^^ 
20.78linS? 


0.000 
-0.001 


Rb 


n 91 c;71 +0.00001 
V.ZIOI 1-0.00002 


n 91 p;71 +0.00001 
U.ZIOI -L-0.00002 


+1 X 10-^ 
X 10-5 


sin^ ei^ 


n 9Q9ni+o.uooo4 

U.ZOZUl_Q 00003 

n 9Q9nfi+o-oooo4 

U.ZOZUU_o.00014 


n 9Q9nn+'J-U0U04 

U.ZOZUU_o.oo004 

n 9Q9nc;+o-oooo4 

U.ZOZUd_o.00014 


+1 X 10-5 
+1 X 10-5 


sin^ C 


n 9QQQ1 +0.00009 
u.zoooi-o.ooooo 

n 9QQQ7+0.00004 
U.ZOOO/_o.00015 


n 9QQQn+o.ooooy 

U.ZOOOU-o.ooOOl 
n 9QQQc;+0.00004 
U.ZOOOd_o.00014 


+ 1 X 10"5 
+2 X 10-5 


A' 

FB 


n ni c;qq+o.oooo8 

U.U1O-3O-0.00007 

01 c;9S+o-ooo24 

U.U10ZO_o.00007 


U.UlOOO-0.00007 
01 eiS1+°-°°°24 
U.Uld01„o.00007 


-3 X 10-5 
-3 X 10-5 




n lAQI 4+0.00028 
U.14014_o.oo030 

n i497e;+o.ooiii 

U.l^Z/d_o.oo030 


n 1 /1Q97+0.00028 
U.J-40Z ^_0.00031 


-1.3 X 10-^ 
-1.4 X 10-^ 


A'' 

FB 


U.1UUZ^_0.00023 

1 nnn4+o-ooo78 

U.lUUU4_o.oo022 


0.10033tHl]l!i 

o.iooi3t°:S 


-9 X 10-5 
-9 X 10-5 
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